/[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.27 by molod, Tue Sep 20 18:30:41 2005 UTC revision 1.28 by molod, Wed Nov 16 16:55:21 2005 UTC
# Line 270  C---+----1----+----2----+----3----+----4 Line 270  C---+----1----+----2----+----3----+----4
270           CALL DIAGNOSTICS_FILL(tmpMk,'WVELSLT ',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'WVELSLT ',0,Nr,0,1,1,myThid)
271          ENDIF          ENDIF
272                
273            IF ( DIAGNOSTICS_IS_ON('UVELPHI ',myThid) ) THEN
274             DO bj = myByLo(myThid), myByHi(myThid)
275              DO bi = myBxLo(myThid), myBxHi(myThid)
276               DO K=1,Nr
277                DO j = 1,sNy
278                 DO i = 1,sNx+1
279                   tmpMk(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*0.5 _d 0
280         &            *(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i-1,j,K,bi,bj))
281                 ENDDO
282                ENDDO
283               ENDDO
284              ENDDO
285             ENDDO
286             CALL DIAGNOSTICS_FILL(tmpMk,'UVELPHI ',0,Nr,0,1,1,myThid)
287            ENDIF
288          
289            IF ( DIAGNOSTICS_IS_ON('VVELPHI ',myThid) ) THEN
290             DO bj = myByLo(myThid), myByHi(myThid)
291              DO bi = myBxLo(myThid), myBxHi(myThid)
292               DO K=1,Nr
293                DO j = 1,sNy+1
294                 DO i = 1,sNx
295                   tmpMk(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*0.5 _d 0
296         &            *(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i,j-1,K,bi,bj))
297                 ENDDO
298                ENDDO
299               ENDDO
300              ENDDO
301             ENDDO
302             CALL DIAGNOSTICS_FILL(tmpMk,'VVELPHI ',0,Nr,0,1,1,myThid)
303            ENDIF
304          
305          IF ( DIAGNOSTICS_IS_ON('PRESSURE',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('PRESSURE',myThid) ) THEN
306           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
307            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22