/[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.31 by jmc, Sun Jun 25 23:06:24 2006 UTC revision 1.33 by jmc, Sun Nov 19 22:03:13 2006 UTC
# Line 277  C---+----1----+----2----+----3----+----4 Line 277  C---+----1----+----2----+----3----+----4
277             DO K=1,Nr             DO K=1,Nr
278              DO j = 1,sNy              DO j = 1,sNy
279               DO i = 1,sNx+1               DO i = 1,sNx+1
280                 tmpMk(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*0.5 _d 0                 tmpMk(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*hFacW(i,j,k,bi,bj)
281       &            *(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i-1,j,K,bi,bj))       &       *0.5 _d 0*(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i-1,j,K,bi,bj))
282               ENDDO               ENDDO
283              ENDDO              ENDDO
284             ENDDO             ENDDO
# Line 293  C---+----1----+----2----+----3----+----4 Line 293  C---+----1----+----2----+----3----+----4
293             DO K=1,Nr             DO K=1,Nr
294              DO j = 1,sNy+1              DO j = 1,sNy+1
295               DO i = 1,sNx               DO i = 1,sNx
296                 tmpMk(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*0.5 _d 0                 tmpMk(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*hFacS(i,j,k,bi,bj)
297       &            *(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i,j-1,K,bi,bj))       &       *0.5 _d 0*(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i,j-1,K,bi,bj))
298               ENDDO               ENDDO
299              ENDDO              ENDDO
300             ENDDO             ENDDO
# Line 303  C---+----1----+----2----+----3----+----4 Line 303  C---+----1----+----2----+----3----+----4
303           CALL DIAGNOSTICS_FILL(tmpMk,'VVELPHI ',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'VVELPHI ',0,Nr,0,1,1,myThid)
304          ENDIF          ENDIF
305    
306          IF ( DIAGNOSTICS_IS_ON('PRESSURE',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('RCENTER ',myThid) ) THEN
307           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
308            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
309             DO j = 1,sNy             DO j = 1,sNy
# Line 323  C---+----1----+----2----+----3----+----4 Line 323  C---+----1----+----2----+----3----+----4
323             ENDDO             ENDDO
324            ENDDO            ENDDO
325           ENDDO           ENDDO
326           CALL DIAGNOSTICS_FILL(tmpMk,'PRESSURE',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'RCENTER ',0,Nr,0,1,1,myThid)
327          ENDIF          ENDIF
328    
329  C First fill sequence for state variable tendency diagnostics: subtract state variable  C First fill sequence for state variable tendency diagnostics: subtract state variable
Legend:
Removed from v.1.31  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22