/[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.18 by dimitri, Thu May 5 16:41:54 2005 UTC revision 1.19 by dimitri, Fri May 6 03:46:40 2005 UTC
# Line 81  C--   fill momentum state-var diagnostic Line 81  C--   fill momentum state-var diagnostic
81          CALL DIAGNOSTICS_FILL(vVel, 'VVEL    ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(vVel, 'VVEL    ',0,Nr,0,1,1,myThid)
82          CALL DIAGNOSTICS_FILL(wVel, 'WVEL    ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(wVel, 'WVEL    ',0,Nr,0,1,1,myThid)
83                
84            IF ( DIAGNOSTICS_IS_ON('UVEL_k2 ',myThid) ) THEN
85             DO bj = myByLo(myThid), myByHi(myThid)
86              DO bi = myBxLo(myThid), myBxHi(myThid)
87               DO j = 1,sNy
88                DO i = 1,sNx
89                  tmp1k(i,j,bi,bj) = UVEL(i,j,2,bi,bj)
90                ENDDO
91               ENDDO
92              ENDDO
93             ENDDO
94             CALL DIAGNOSTICS_FILL(tmp1k,'UVEL_k2 ',0,1,0,1,1,myThid)
95            ENDIF
96          
97            IF ( DIAGNOSTICS_IS_ON('VVEL_k2 ',myThid) ) THEN
98             DO bj = myByLo(myThid), myByHi(myThid)
99              DO bi = myBxLo(myThid), myBxHi(myThid)
100               DO j = 1,sNy
101                DO i = 1,sNx
102                  tmp1k(i,j,bi,bj) = VVEL(i,j,2,bi,bj)
103                ENDDO
104               ENDDO
105              ENDDO
106             ENDDO
107             CALL DIAGNOSTICS_FILL(tmp1k,'VVEL_k2 ',0,1,0,1,1,myThid)
108            ENDIF
109          
110          IF ( DIAGNOSTICS_IS_ON('UVELSQ  ',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('UVELSQ  ',myThid) ) THEN
111           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
112            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 312  C--   fill tracer state-var diagnostics: Line 338  C--   fill tracer state-var diagnostics:
338    
339          CALL DIAGNOSTICS_FILL(theta,'THETA   ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(theta,'THETA   ',0,Nr,0,1,1,myThid)
340          CALL DIAGNOSTICS_FILL(salt, 'SALT    ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(salt, 'SALT    ',0,Nr,0,1,1,myThid)
341          
342            IF ( DIAGNOSTICS_IS_ON('SST     ',myThid) ) THEN
343             DO bj = myByLo(myThid), myByHi(myThid)
344              DO bi = myBxLo(myThid), myBxHi(myThid)
345               DO j = 1,sNy
346                DO i = 1,sNx
347                  tmp1k(i,j,bi,bj) = THETA(i,j,1,bi,bj)
348                ENDDO
349               ENDDO
350              ENDDO
351             ENDDO
352             CALL DIAGNOSTICS_FILL(tmp1k,'SST     ',0,1,0,1,1,myThid)
353            ENDIF
354          
355            IF ( DIAGNOSTICS_IS_ON('SSS     ',myThid) ) THEN
356             DO bj = myByLo(myThid), myByHi(myThid)
357              DO bi = myBxLo(myThid), myBxHi(myThid)
358               DO j = 1,sNy
359                DO i = 1,sNx
360                  tmp1k(i,j,bi,bj) = SALT(i,j,1,bi,bj)
361                ENDDO
362               ENDDO
363              ENDDO
364             ENDDO
365             CALL DIAGNOSTICS_FILL(tmp1k,'SSS     ',0,1,0,1,1,myThid)
366            ENDIF
367    
368          IF ( DIAGNOSTICS_IS_ON('THETASQ ',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('THETASQ ',myThid) ) THEN
369           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22