/[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.34 by adcroft, Mon Sep 23 16:13:31 2002 UTC revision 1.35 by mlosch, Wed Sep 25 19:36:50 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        _RL tmpFac, convertVol2Mass
58        INTEGER numIters        INTEGER numIters
59        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
60  CEOP  CEOP
61    
62          IF ( buoyancyRelation .eq. 'OCEANICP' )  THEN
63           convertVol2Mass = gravity*rhoConstFresh
64          ELSE
65           convertVol2Mass = 1. _d 0
66          ENDIF
67    
68  C--   Save previous solution & Initialise Vector solution and source term :  C--   Save previous solution & Initialise Vector solution and source term :
69        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
70         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
# Line 72  C--   Save previous solution & Initialis Line 78  C--   Save previous solution & Initialis
78           ENDDO           ENDDO
79          ENDDO          ENDDO
80          IF (useRealFreshWaterFlux) THEN          IF (useRealFreshWaterFlux) THEN
81           tmpFac = freeSurfFac           tmpFac = freeSurfFac*convertVol2Mass
82           IF (exactConserv) tmpFac = freeSurfFac*implicDiv2DFlow           IF (exactConserv)
83           IF ( buoyancyRelation .eq. 'OCEANICP' )  THEN       &        tmpFac = freeSurfFac*convertVol2Mass*implicDiv2DFlow
            tmpFac = tmpFac*horiVertRatio  
          ENDIF  
84           DO j=1,sNy           DO j=1,sNy
85            DO i=1,sNx            DO i=1,sNx
86             cg2d_b(i,j,bi,bj) =             cg2d_b(i,j,bi,bj) =
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.22