/[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.34 by jmc, Sun Dec 24 20:07:02 2006 UTC revision 1.35 by jmc, Mon Feb 5 03:14:26 2007 UTC
# Line 337  C     (next fill for these diagnostics w Line 337  C     (next fill for these diagnostics w
337            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
338             DO K=1,Nr             DO K=1,Nr
339              DO j = 1,sNy              DO j = 1,sNy
340               DO i = 1,sNx               DO i = 1,sNx+1
341                 tmpMk(i,j,K,bi,bj) = -uVel(i,j,K,bi,bj)*                 tmpMk(i,j,K,bi,bj) = -uVel(i,j,K,bi,bj)*
342       .                    *86400./dTtracerLev(1)       .                    *86400./dTtracerLev(1)
343               ENDDO               ENDDO
# Line 352  C     (next fill for these diagnostics w Line 352  C     (next fill for these diagnostics w
352           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
353            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
354             DO K=1,Nr             DO K=1,Nr
355              DO j = 1,sNy              DO j = 1,sNy+1
356               DO i = 1,sNx               DO i = 1,sNx
357                 tmpMk(i,j,K,bi,bj) = -vVel(i,j,K,bi,bj)                 tmpMk(i,j,K,bi,bj) = -vVel(i,j,K,bi,bj)
358       .                    *86400./dTtracerLev(1)       .                    *86400./dTtracerLev(1)
# Line 486  C---+----1----+----2----+----3----+----4 Line 486  C---+----1----+----2----+----3----+----4
486            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
487             DO K=1,Nr             DO K=1,Nr
488              DO j = 1,sNy              DO j = 1,sNy
489                DO i = 1,sNx                DO i = 1,sNx+1
490                  tmpMk(i,j,K,bi,bj)                  tmpMk(i,j,K,bi,bj)
491       &               = uVel(i,j,K,bi,bj)*hFacW(i,j,K,bi,bj)       &               = uVel(i,j,K,bi,bj)*hFacW(i,j,K,bi,bj)
492                ENDDO                ENDDO
# Line 501  C---+----1----+----2----+----3----+----4 Line 501  C---+----1----+----2----+----3----+----4
501           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
502            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
503             DO K=1,Nr             DO K=1,Nr
504              DO j = 1,sNy              DO j = 1,sNy+1
505                DO i = 1,sNx                DO i = 1,sNx
506                  tmpMk(i,j,K,bi,bj)                  tmpMk(i,j,K,bi,bj)
507       &               = vVel(i,j,K,bi,bj)*hFacS(i,j,K,bi,bj)       &               = vVel(i,j,K,bi,bj)*hFacS(i,j,K,bi,bj)
# Line 630  C       this means that diag fill is cal Line 630  C       this means that diag fill is cal
630            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
631             DO K=1,Nr             DO K=1,Nr
632              DO j = 1,sNy              DO j = 1,sNy
633               DO i = 1,sNx               DO i = 1,sNx+1
634                 tmpMk(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)                 tmpMk(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)
635       .                    *86400./dTtracerLev(1)       .                    *86400./dTtracerLev(1)
636               ENDDO               ENDDO
# Line 645  C       this means that diag fill is cal Line 645  C       this means that diag fill is cal
645           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
646            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
647             DO K=1,Nr             DO K=1,Nr
648              DO j = 1,sNy              DO j = 1,sNy+1
649               DO i = 1,sNx               DO i = 1,sNx
650                 tmpMk(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)                 tmpMk(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)
651       .                    *86400./dTtracerLev(1)       .                    *86400./dTtracerLev(1)
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.22