/[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.32 by dfer, Tue Aug 29 17:00:30 2006 UTC revision 1.34 by jmc, Sun Dec 24 20:07:02 2006 UTC
# 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 315  C---+----1----+----2----+----3----+----4 Line 315  C---+----1----+----2----+----3----+----4
315              DO j = 1,sNy              DO j = 1,sNy
316               DO i = 1,sNx               DO i = 1,sNx
317                tmpMk(i,j,k,bi,bj) = tmp1k(i,j,bi,bj)                tmpMk(i,j,k,bi,bj) = tmp1k(i,j,bi,bj)
318       &                      + drF(k)*hFacC(i,j,k,bi,bj)*0.5 _d 0       &               + (rF(k)-rC(k))*hFacC(i,j,k,bi,bj)
319    C         above: more general (setInterFDr/setCenterDr) than line below
320    c    &                      + drF(k)*hFacC(i,j,k,bi,bj)*0.5 _d 0
321                tmp1k(i,j,bi,bj) =   tmp1k(i,j,bi,bj)                tmp1k(i,j,bi,bj) =   tmp1k(i,j,bi,bj)
322       &                      + drF(k)*hFacC(i,j,k,bi,bj)       &                      + drF(k)*hFacC(i,j,k,bi,bj)
323               ENDDO               ENDDO
# Line 323  C---+----1----+----2----+----3----+----4 Line 325  C---+----1----+----2----+----3----+----4
325             ENDDO             ENDDO
326            ENDDO            ENDDO
327           ENDDO           ENDDO
328           CALL DIAGNOSTICS_FILL(tmpMk,'PRESSURE',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'RCENTER ',0,Nr,0,1,1,myThid)
329          ENDIF          ENDIF
330    
331  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.32  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22