Ocean Modeling Discussion

thank you kate, In the estuary ,i set some sections from head to mouth. And there is no branch, so the tidally averaged water flux in the each section should be equal, at least there is almost the same. But the water flux in the each section is different and from head to mouth the value decrease ...

I made a mistake.
I missed the control line "IF (Zgrd.ge.0.0_r8) THEN" that is set before the early posted "Program Lines" in SUBROUTINE extract_sta.F.
It write the result in a S-coor.
The code is right.

I need to write out the result at a station, i found a bug in the SUBROUTINE extract_sta3d. ! !!!!Interpolate from 3D field at RHO-points.: DO k=N(ng),2,-1 Ztop=GRID(ng)%z_r(i1,j1,k) Zbot=GRID(ng)%z_r(i1,j1,k-1) IF ((Ztop.gt.Zgrd).and.(Zgrd.ge.Zbot)) THEN k1=k-1 k2=k END IF END DO the Zgrd is from ...