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

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

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

revision 1.27 by cnh, Wed Sep 26 18:09:16 2001 UTC revision 1.28 by jmc, Sun Jan 26 21:06:11 2003 UTC
# Line 72  C-    Compute effective gU term (includi Line 72  C-    Compute effective gU term (includi
72  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
73        IF (.NOT. vectorInvariantMomentum        IF (.NOT. vectorInvariantMomentum
74       &    .AND. nonlinFreeSurf.GT.1) THEN       &    .AND. nonlinFreeSurf.GT.1) THEN
75         DO j=jMin,jMax         IF (select_rStar.GT.0) THEN
76          DO i=iMin,iMax          DO j=jMin,jMax
77           IF ( k.EQ.ksurfW(i,j,bi,bj) ) THEN           DO i=iMin,iMax
78               gUtmp(i,j) = gUtmp(i,j)/rStarExpW(i,j,bi,bj)
79             ENDDO
80            ENDDO
81           ELSE
82            DO j=jMin,jMax
83             DO i=iMin,iMax
84              IF ( k.EQ.ksurfW(i,j,bi,bj) ) THEN
85             gUtmp(i,j) = gUtmp(i,j)             gUtmp(i,j) = gUtmp(i,j)
86       &         *hFacW(i,j,k,bi,bj)/hFac_surfW(i,j,bi,bj)       &         *hFacW(i,j,k,bi,bj)/hFac_surfW(i,j,bi,bj)
87           ENDIF            ENDIF
88             ENDDO
89          ENDDO          ENDDO
90         ENDDO         ENDIF
91        ENDIF        ENDIF
92  #endif  #endif
93    
# Line 123  C-    Compute effective gV term (includi Line 131  C-    Compute effective gV term (includi
131  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
132        IF (.NOT. vectorInvariantMomentum        IF (.NOT. vectorInvariantMomentum
133       &    .AND. nonlinFreeSurf.GT.1) THEN       &    .AND. nonlinFreeSurf.GT.1) THEN
134         DO j=jMin,jMax         IF (select_rStar.GT.0) THEN
135          DO i=iMin,iMax          DO j=jMin,jMax
136           IF ( k.EQ.ksurfS(i,j,bi,bj) ) THEN           DO i=iMin,iMax
137               gVtmp(i,j) = gVtmp(i,j)/rStarExpS(i,j,bi,bj)
138             ENDDO
139            ENDDO
140           ELSE
141            DO j=jMin,jMax
142             DO i=iMin,iMax
143              IF ( k.EQ.ksurfS(i,j,bi,bj) ) THEN
144             gVtmp(i,j) = gVtmp(i,j)             gVtmp(i,j) = gVtmp(i,j)
145       &         *hFacS(i,j,k,bi,bj)/hFac_surfS(i,j,bi,bj)       &         *hFacS(i,j,k,bi,bj)/hFac_surfS(i,j,bi,bj)
146           ENDIF            ENDIF
147             ENDDO
148          ENDDO          ENDDO
149         ENDDO         ENDIF
150        ENDIF        ENDIF
151  #endif  #endif
152    
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22