/[MITgcm]/MITgcm/model/src/solve_for_pressure.F
ViewVC logotype

Diff of /MITgcm/model/src/solve_for_pressure.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.69 by mlosch, Mon Nov 23 16:15:54 2009 UTC revision 1.70 by jmc, Wed Nov 25 20:56:15 2009 UTC
# Line 163  C--   Add source term arising from w=d/d Line 163  C--   Add source term arising from w=d/d
163        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
164         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
165  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
166    C--   Add EmPmR contribution to top level cg3d_b:
167    C      (has been done for cg2d_b ; and addMass was added by CALC_DIV_GHAT)
168            IF ( use3Dsolver .AND.
169         &       useRealFreshWaterFlux.AND.fluidIsWater ) THEN
170              tmpFac = freeSurfFac*mass2rUnit
171              IF (exactConserv)
172         &     tmpFac = freeSurfFac*mass2rUnit*implicDiv2DFlow
173              ks = 1
174              IF ( usingPCoords ) ks = Nr
175              DO j=1,sNy
176               DO i=1,sNx
177                cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
178         &        + tmpFac*_rA(i,j,bi,bj)*EmPmR(i,j,bi,bj)/deltaTMom
179               ENDDO
180              ENDDO
181            ENDIF
182          IF ( use3Dsolver .AND. zeroPsNH ) THEN          IF ( use3Dsolver .AND. zeroPsNH ) THEN
183           DO j=1,sNy           DO j=1,sNy
184            DO i=1,sNx            DO i=1,sNx
Legend:
Removed from v.1.69  
changed lines
  Added in v.1.70
  ViewVC Help
Powered by ViewVC 1.1.22