/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_fill_state.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_fill_state.F

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

revision 1.15 by jmc, Fri Jan 28 01:05:02 2005 UTC revision 1.16 by jmc, Thu Feb 10 15:01:34 2005 UTC
# Line 162  C---+----1----+----2----+----3----+----4 Line 162  C---+----1----+----2----+----3----+----4
162           CALL DIAGNOSTICS_FILL(tmpMk,'UV_VEL_Z',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'UV_VEL_Z',0,Nr,0,1,1,myThid)
163          ENDIF          ENDIF
164                
165            IF ( DIAGNOSTICS_IS_ON('WU_VEL  ',myThid) ) THEN
166             DO bj = myByLo(myThid), myByHi(myThid)
167              DO bi = myBxLo(myThid), myBxHi(myThid)
168               DO K=1,Nr
169                km1 = MAX(k-1,1)
170                DO j = 1,sNy
171                 DO i = 1,sNx+1
172                  tmpMk(i,j,K,bi,bj) = 0.25 _d 0
173         &           *(uVel(i,km1,K,bi,bj)+uVel(i,j,K,bi,bj))
174         &           *(wVel(i-1,j,K,bi,bj)*rA(i-1,j,bi,bj)
175         &            +wVel( i ,j,K,bi,bj)*rA( i ,j,bi,bj)
176         &            )*recip_rAw(i,j,bi,bj)
177                 ENDDO
178                ENDDO
179               ENDDO
180              ENDDO
181             ENDDO
182             CALL DIAGNOSTICS_FILL(tmpMk,'WU_VEL  ',0,Nr,0,1,1,myThid)
183            ENDIF
184    
185            IF ( DIAGNOSTICS_IS_ON('WV_VEL  ',myThid) ) THEN
186             DO bj = myByLo(myThid), myByHi(myThid)
187              DO bi = myBxLo(myThid), myBxHi(myThid)
188               DO K=1,Nr
189                km1 = MAX(k-1,1)
190                DO j = 1,sNy+1
191                 DO i = 1,sNx
192                  tmpMk(i,j,K,bi,bj) = 0.25 _d 0
193         &           *(vVel(i,km1,K,bi,bj)+vVel(i,j,K,bi,bj))
194         &           *(wVel(i,j-1,K,bi,bj)*rA(i,j-1,bi,bj)
195         &            +wVel(i, j ,K,bi,bj)*rA(i, j ,bi,bj)
196         &            )*recip_rAs(i,j,bi,bj)
197                 ENDDO
198                ENDDO
199               ENDDO
200              ENDDO
201             ENDDO
202             CALL DIAGNOSTICS_FILL(tmpMk,'WV_VEL  ',0,Nr,0,1,1,myThid)
203            ENDIF
204    
205  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
206    
207          IF ( DIAGNOSTICS_IS_ON('UVELTH  ',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('UVELTH  ',myThid) ) THEN
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22