/[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.35 by mlosch, Wed Sep 25 19:36:50 2002 UTC revision 1.36 by jmc, Mon Oct 7 16:20:39 2002 UTC
# Line 54  C     == Local variables == Line 54  C     == Local variables ==
54        _RS uf(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RS uf(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
55        _RS vf(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RS vf(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
56        _RL firstResidual,lastResidual        _RL firstResidual,lastResidual
57        _RL tmpFac, convertVol2Mass        _RL tmpFac
58        INTEGER numIters        INTEGER numIters
59        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
60  CEOP  CEOP
61    
       IF ( buoyancyRelation .eq. 'OCEANICP' )  THEN  
        convertVol2Mass = gravity*rhoConstFresh  
       ELSE  
        convertVol2Mass = 1. _d 0  
       ENDIF  
   
62  C--   Save previous solution & Initialise Vector solution and source term :  C--   Save previous solution & Initialise Vector solution and source term :
63        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
64         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 78  C--   Save previous solution & Initialis Line 72  C--   Save previous solution & Initialis
72           ENDDO           ENDDO
73          ENDDO          ENDDO
74          IF (useRealFreshWaterFlux) THEN          IF (useRealFreshWaterFlux) THEN
75           tmpFac = freeSurfFac*convertVol2Mass           tmpFac = freeSurfFac*convertEmP2rUnit
76           IF (exactConserv)           IF (exactConserv)
77       &        tmpFac = freeSurfFac*convertVol2Mass*implicDiv2DFlow       &        tmpFac = freeSurfFac*convertEmP2rUnit*implicDiv2DFlow
78           DO j=1,sNy           DO j=1,sNy
79            DO i=1,sNx            DO i=1,sNx
80             cg2d_b(i,j,bi,bj) =             cg2d_b(i,j,bi,bj) =
Legend:
Removed from v.1.35  
changed lines
  Added in v.1.36
  ViewVC Help
Powered by ViewVC 1.1.22