/[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.61 by jmc, Mon Aug 27 13:18:31 2007 UTC revision 1.62 by jmc, Mon Oct 15 15:28:24 2007 UTC
# Line 105  C--   Save previous solution & Initialis Line 105  C--   Save previous solution & Initialis
105           ENDDO           ENDDO
106          ENDDO          ENDDO
107          IF (useRealFreshWaterFlux) THEN          IF (useRealFreshWaterFlux) THEN
108           tmpFac = freeSurfFac*convertEmP2rUnit           tmpFac = freeSurfFac*mass2rUnit
109           IF (exactConserv)           IF (exactConserv)
110       &        tmpFac = freeSurfFac*convertEmP2rUnit*implicDiv2DFlow       &        tmpFac = freeSurfFac*mass2rUnit*implicDiv2DFlow
111           DO j=1,sNy           DO j=1,sNy
112            DO i=1,sNx            DO i=1,sNx
113             cg2d_b(i,j,bi,bj) =             cg2d_b(i,j,bi,bj) =
# Line 135  C--   Save previous solution & Initialis Line 135  C--   Save previous solution & Initialis
135         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
136          IF ( putPmEinXvector ) THEN          IF ( putPmEinXvector ) THEN
137            tmpFac = 0.            tmpFac = 0.
138            IF (globalArea.GT.0.) tmpFac = freeSurfFac*deltaTfreesurf            IF (globalArea.GT.0.) tmpFac =
139       &                          *convertEmP2rUnit*sumEmP/globalArea       &      freeSurfFac*deltaTfreesurf*mass2rUnit*sumEmP/globalArea
140            DO j=1,sNy            DO j=1,sNy
141             DO i=1,sNx             DO i=1,sNx
142              cg2d_x(i,j,bi,bj) = cg2d_x(i,j,bi,bj)              cg2d_x(i,j,bi,bj) = cg2d_x(i,j,bi,bj)
Legend:
Removed from v.1.61  
changed lines
  Added in v.1.62
  ViewVC Help
Powered by ViewVC 1.1.22