/[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.28 by jmc, Fri Feb 8 22:13:39 2002 UTC revision 1.29 by jmc, Sun Feb 10 00:39:22 2002 UTC
# Line 6  C $Name$ Line 6  C $Name$
6  CBOP  CBOP
7  C     !ROUTINE: SOLVE_FOR_PRESSURE  C     !ROUTINE: SOLVE_FOR_PRESSURE
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE SOLVE_FOR_PRESSURE(myTime, myIter, myThid)        SUBROUTINE SOLVE_FOR_PRESSURE(myTime, myIter, myThid)
10    
11  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
12  C     *==========================================================*  C     *==========================================================*
# Line 42  C     myIter - Current iteration number Line 42  C     myIter - Current iteration number
42  C     myThid - Thread number for this instance of SOLVE_FOR_PRESSURE  C     myThid - Thread number for this instance of SOLVE_FOR_PRESSURE
43        _RL myTime        _RL myTime
44        INTEGER myIter        INTEGER myIter
45        INTEGER myThid        INTEGER myThid
46    
47  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
48  C     == Local variables ==  C     == Local variables ==
# Line 50  C     == Local variables == Line 50  C     == Local variables ==
50        _RS uf(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RS uf(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
51        _RS vf(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RS vf(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
52        _RL firstResidual,lastResidual        _RL firstResidual,lastResidual
53          _RL tmpFac
54        INTEGER numIters        INTEGER numIters
55        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
56  CEOP  CEOP
# Line 64  C--   Save previous solution & Initialis Line 65  C--   Save previous solution & Initialis
65  #endif  #endif
66            cg2d_x(i,j,bi,bj) = Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)            cg2d_x(i,j,bi,bj) = Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)
67            cg2d_b(i,j,bi,bj) = 0.            cg2d_b(i,j,bi,bj) = 0.
 #ifdef USE_NATURAL_BCS  
      &     + freeSurfFac*_rA(i,j,bi,bj)*  
      &       EmPmR(I,J,bi,bj)/deltaTMom  
 #endif  
68           ENDDO           ENDDO
69          ENDDO          ENDDO
70            IF (useRealFreshWaterFlux) THEN
71             tmpFac = freeSurfFac
72             IF (exactConserv) tmpFac = freeSurfFac*implicDiv2DFlow
73             DO j=1,sNy
74              DO i=1,sNx
75               cg2d_b(i,j,bi,bj) =
76         &       tmpFac*_rA(i,j,bi,bj)*EmPmR(i,j,bi,bj)/deltaTMom
77              ENDDO
78             ENDDO
79            ENDIF
80         ENDDO         ENDDO
81        ENDDO        ENDDO
82    
Legend:
Removed from v.1.28  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.22