/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_calc_rhs.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_calc_rhs.F

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

revision 1.11 by adcroft, Wed Sep 19 20:45:09 2001 UTC revision 1.12 by heimbach, Thu Sep 27 20:12:11 2001 UTC
# Line 104  C--   the kDown is still required Line 104  C--   the kDown is still required
104  #endif  #endif
105        DO j=1-OLy,sNy+OLy        DO j=1-OLy,sNy+OLy
106         DO i=1-OLx,sNx+OLx         DO i=1-OLx,sNx+OLx
107          fZon(i,j)      = 0.0          fZon(i,j)      = 0. _d 0
108          fMer(i,j)      = 0.0          fMer(i,j)      = 0. _d 0
109          fVerT(i,j,kUp) = 0.0          fVerT(i,j,kUp) = 0. _d 0
110         ENDDO         ENDDO
111        ENDDO        ENDDO
112    
# Line 125  C     the tendency to zero. Line 125  C     the tendency to zero.
125       &    advectionScheme.EQ.ENUM_CENTERED_4TH ) THEN       &    advectionScheme.EQ.ENUM_CENTERED_4TH ) THEN
126         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
127          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
128           gTracer(i,j,k,bi,bj)=0.           gTracer(i,j,k,bi,bj)=0. _d 0
129          ENDDO          ENDDO
130         ENDDO         ENDDO
131        ENDIF        ENDIF
# Line 140  C--   Pre-calculate del^2 T if bi-harmon Line 140  C--   Pre-calculate del^2 T if bi-harmon
140  C--   Initialize net flux in X direction  C--   Initialize net flux in X direction
141        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
142         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
143          fZon(i,j) = 0.          fZon(i,j) = 0. _d 0
144         ENDDO         ENDDO
145        ENDDO        ENDDO
146    
# Line 180  C-    Diffusive flux in X Line 180  C-    Diffusive flux in X
180        ELSE        ELSE
181         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
182          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
183           df(i,j) = 0.           df(i,j) = 0. _d 0
184          ENDDO          ENDDO
185         ENDDO         ENDDO
186        ENDIF        ENDIF
# Line 215  C-    Bi-harmonic duffusive flux in X Line 215  C-    Bi-harmonic duffusive flux in X
215  C--   Initialize net flux in Y direction  C--   Initialize net flux in Y direction
216        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
217         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
218          fMer(i,j) = 0.          fMer(i,j) = 0. _d 0
219         ENDDO         ENDDO
220        ENDDO        ENDDO
221    
# Line 255  C-    Diffusive flux in Y Line 255  C-    Diffusive flux in Y
255        ELSE        ELSE
256         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
257          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
258           df(i,j) = 0.           df(i,j) = 0. _d 0
259          ENDDO          ENDDO
260         ENDDO         ENDDO
261        ENDIF        ENDIF
# Line 287  C-    Bi-harmonic flux in Y Line 287  C-    Bi-harmonic flux in Y
287         ENDDO         ENDDO
288        ENDIF        ENDIF
289    
 C--   Initialize net flux in R  
       DO j=1-Oly,sNy+Oly  
        DO i=1-Olx,sNx+Olx  
         fVerT(i,j,kUp) = 0.  
        ENDDO  
       ENDDO  
   
290  C-    Advective flux in R  C-    Advective flux in R
291        IF (.NOT. multiDimAdvection .OR.        IF (.NOT. multiDimAdvection .OR.
292       &    advectionScheme.EQ.ENUM_CENTERED_2ND .OR.       &    advectionScheme.EQ.ENUM_CENTERED_2ND .OR.
# Line 350  C           boundary condition. Line 343  C           boundary condition.
343        IF (implicitDiffusion) THEN        IF (implicitDiffusion) THEN
344         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
345          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
346           df(i,j) = 0.           df(i,j) = 0. _d 0
347          ENDDO          ENDDO
348         ENDDO         ENDDO
349        ELSE        ELSE
# Line 390  C-    Add non local KPP transport term ( Line 383  C-    Add non local KPP transport term (
383        IF (useKPP) THEN        IF (useKPP) THEN
384         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
385          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
386           df(i,j) = 0.           df(i,j) = 0. _d 0
387          ENDDO          ENDDO
388         ENDDO         ENDDO
389         IF (tracerIdentity.EQ.GAD_TEMPERATURE) THEN         IF (tracerIdentity.EQ.GAD_TEMPERATURE) THEN
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22