/[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.63 by jmc, Fri Aug 22 16:04:48 2008 UTC revision 1.64 by jmc, Sun Aug 24 21:40:18 2008 UTC
# Line 93  C--   Initialise the Vector solution wit Line 93  C--   Initialise the Vector solution wit
93  C     instead of simply etaN ; This can speed-up the solver convergence in  C     instead of simply etaN ; This can speed-up the solver convergence in
94  C     the case where |Global_mean_PmE| is large.  C     the case where |Global_mean_PmE| is large.
95        putPmEinXvector = .FALSE.        putPmEinXvector = .FALSE.
96  c     putPmEinXvector = useRealFreshWaterFlux  c     putPmEinXvector = useRealFreshWaterFlux.AND.fluidIsWater
97    
98  C--   Save previous solution & Initialise Vector solution and source term :  C--   Save previous solution & Initialise Vector solution and source term :
99        sumEmP = 0.        sumEmP = 0.
# Line 108  C--   Save previous solution & Initialis Line 108  C--   Save previous solution & Initialis
108            cg2d_b(i,j,bi,bj) = 0.            cg2d_b(i,j,bi,bj) = 0.
109           ENDDO           ENDDO
110          ENDDO          ENDDO
111          IF (useRealFreshWaterFlux) THEN          IF (useRealFreshWaterFlux.AND.fluidIsWater) THEN
112           tmpFac = freeSurfFac*mass2rUnit           tmpFac = freeSurfFac*mass2rUnit
113           IF (exactConserv)           IF (exactConserv)
114       &        tmpFac = freeSurfFac*mass2rUnit*implicDiv2DFlow       &        tmpFac = freeSurfFac*mass2rUnit*implicDiv2DFlow
Legend:
Removed from v.1.63  
changed lines
  Added in v.1.64
  ViewVC Help
Powered by ViewVC 1.1.22