/[MITgcm]/MITgcm/model/src/integrate_for_w.F
ViewVC logotype

Diff of /MITgcm/model/src/integrate_for_w.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.8 by jmc, Mon Oct 7 16:24:45 2002 UTC revision 1.9 by jmc, Sun Jan 26 21:06:11 2003 UTC
# Line 25  C     == GLobal variables == Line 25  C     == GLobal variables ==
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
26  #include "PARAMS.h"  #include "PARAMS.h"
27  #include "GRID.h"  #include "GRID.h"
28    #include "SURFACE.h"
29    
30  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
31  C     == Routine arguments ==  C     == Routine arguments ==
# Line 89  C-  o Rigid-Lid case: zero at lower and Line 90  C-  o Rigid-Lid case: zero at lower and
90             ENDDO             ENDDO
91            ENDDO            ENDDO
92          ENDIF          ENDIF
93    #ifdef NONLIN_FRSURF
94          ELSEIF (select_rStar .NE. 0) THEN
95    C-  o rStar case: zero under-ground and at r_lower boundary
96    C     can be non-zero at surface (useRealFreshWaterFlux)
97            IF (k.eq.Nr) THEN
98              DO j=1,sNy
99               DO i=1,sNx
100                 wFld(i,j,k,bi,bj) =
101         &        (       -( uTrans(i+1,j)-uTrans(i,j)
102         &                  +vTrans(i,j+1)-vTrans(i,j)
103         &                 )*recip_rA(i,j,bi,bj)
104         &          - rStarDhCDt(i,j,bi,bj)*drF(k)*h0FacC(i,j,k,bi,bj)
105         &         )*maskC(i,j,k,bi,bj)
106               ENDDO
107              ENDDO
108            ELSE
109              DO j=1,sNy
110               DO i=1,sNx
111                 wFld(i,j,k,bi,bj) =
112         &        ( wFld(i,j,k+1,bi,bj)
113         &                -( uTrans(i+1,j)-uTrans(i,j)
114         &                  +vTrans(i,j+1)-vTrans(i,j)
115         &                 )*recip_rA(i,j,bi,bj)
116         &          - rStarDhCDt(i,j,bi,bj)*drF(k)*h0FacC(i,j,k,bi,bj)
117         &         )*maskC(i,j,k,bi,bj)
118               ENDDO
119              ENDDO
120            ENDIF
121    #endif /* NONLIN_FRSURF */
122        ELSE        ELSE
123  C-  o Free Surface case:  C-  o Free Surface case (r-Coordinate):
124  C      non zero at surface ; zero under-ground and at r_lower boundary  C      non zero at surface ; zero under-ground and at r_lower boundary
125          IF (k.eq.Nr) THEN          IF (k.eq.Nr) THEN
126            DO j=1,sNy            DO j=1,sNy
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22