/[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.17 by jmc, Thu Feb 10 16:30:39 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)
# Line 344  C--   fill tracer state-var diagnostics: Line 396  C--   fill tracer state-var diagnostics:
396           CALL DIAGNOSTICS_FILL(tmpMk,'SALTSQ  ',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'SALTSQ  ',0,Nr,0,1,1,myThid)
397          ENDIF          ENDIF
398                
399            IF ( DIAGNOSTICS_IS_ON('SALTanom',myThid) ) THEN
400             DO bj = myByLo(myThid), myByHi(myThid)
401              DO bi = myBxLo(myThid), myBxHi(myThid)
402               DO K=1,Nr
403                DO j = 1,sNy
404                 DO i = 1,sNx
405                   tmpMk(i,j,K,bi,bj) = salt(i,j,K,bi,bj)-35
406                 ENDDO
407                ENDDO
408               ENDDO
409              ENDDO
410             ENDDO
411             CALL DIAGNOSTICS_FILL(tmpMk,'SALTanom',0,Nr,0,1,1,myThid)
412            ENDIF
413          
414            IF ( DIAGNOSTICS_IS_ON('SALTSQan',myThid) ) THEN
415             DO bj = myByLo(myThid), myByHi(myThid)
416              DO bi = myBxLo(myThid), myBxHi(myThid)
417               DO K=1,Nr
418                DO j = 1,sNy
419                 DO i = 1,sNx
420                   tmpMk(i,j,K,bi,bj) =
421         &               (salt(i,j,K,bi,bj)-35)*(salt(i,j,K,bi,bj)-35)
422                 ENDDO
423                ENDDO
424               ENDDO
425              ENDDO
426             ENDDO
427             CALL DIAGNOSTICS_FILL(tmpMk,'SALTSQan',0,Nr,0,1,1,myThid)
428            ENDIF
429          
430  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
431    
432          IF ( DIAGNOSTICS_IS_ON('UVELMASS',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('UVELMASS',myThid) ) THEN
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22