/[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.15 by jmc, Fri Jan 28 01:05:02 2005 UTC revision 1.18 by dimitri, Thu May 5 16:41:54 2005 UTC
# Line 162  C---+----1----+----2----+----3----+----4 Line 162  C---+----1----+----2----+----3----+----4
162           CALL DIAGNOSTICS_FILL(tmpMk,'UV_VEL_Z',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'UV_VEL_Z',0,Nr,0,1,1,myThid)
163          ENDIF          ENDIF
164                
165            IF ( DIAGNOSTICS_IS_ON('WU_VEL  ',myThid) ) THEN
166             DO bj = myByLo(myThid), myByHi(myThid)
167              DO bi = myBxLo(myThid), myBxHi(myThid)
168               DO K=1,Nr
169                km1 = MAX(k-1,1)
170                DO j = 1,sNy
171                 DO i = 1,sNx+1
172                  tmpMk(i,j,K,bi,bj) = 0.25 _d 0
173         &           *(uVel(i,j,km1,bi,bj)+uVel(i,j,K,bi,bj))
174         &           *(wVel(i-1,j,K,bi,bj)*rA(i-1,j,bi,bj)
175         &            +wVel( i ,j,K,bi,bj)*rA( i ,j,bi,bj)
176         &            )*recip_rAw(i,j,bi,bj)
177                 ENDDO
178                ENDDO
179               ENDDO
180              ENDDO
181             ENDDO
182             CALL DIAGNOSTICS_FILL(tmpMk,'WU_VEL  ',0,Nr,0,1,1,myThid)
183            ENDIF
184    
185            IF ( DIAGNOSTICS_IS_ON('WV_VEL  ',myThid) ) THEN
186             DO bj = myByLo(myThid), myByHi(myThid)
187              DO bi = myBxLo(myThid), myBxHi(myThid)
188               DO K=1,Nr
189                km1 = MAX(k-1,1)
190                DO j = 1,sNy+1
191                 DO i = 1,sNx
192                  tmpMk(i,j,K,bi,bj) = 0.25 _d 0
193         &           *(vVel(i,j,km1,bi,bj)+vVel(i,j,K,bi,bj))
194         &           *(wVel(i,j-1,K,bi,bj)*rA(i,j-1,bi,bj)
195         &            +wVel(i, j ,K,bi,bj)*rA(i, j ,bi,bj)
196         &            )*recip_rAs(i,j,bi,bj)
197                 ENDDO
198                ENDDO
199               ENDDO
200              ENDDO
201             ENDDO
202             CALL DIAGNOSTICS_FILL(tmpMk,'WV_VEL  ',0,Nr,0,1,1,myThid)
203            ENDIF
204    
205  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
206    
207          IF ( DIAGNOSTICS_IS_ON('UVELTH  ',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('UVELTH  ',myThid) ) THEN
# Line 304  C--   fill tracer state-var diagnostics: Line 344  C--   fill tracer state-var diagnostics:
344           CALL DIAGNOSTICS_FILL(tmpMk,'SALTSQ  ',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'SALTSQ  ',0,Nr,0,1,1,myThid)
345          ENDIF          ENDIF
346                
347            IF ( DIAGNOSTICS_IS_ON('SALTanom',myThid) ) THEN
348             DO bj = myByLo(myThid), myByHi(myThid)
349              DO bi = myBxLo(myThid), myBxHi(myThid)
350               DO K=1,Nr
351                DO j = 1,sNy
352                 DO i = 1,sNx
353                   tmpMk(i,j,K,bi,bj) = salt(i,j,K,bi,bj)-35
354                 ENDDO
355                ENDDO
356               ENDDO
357              ENDDO
358             ENDDO
359             CALL DIAGNOSTICS_FILL(tmpMk,'SALTanom',0,Nr,0,1,1,myThid)
360            ENDIF
361          
362            IF ( DIAGNOSTICS_IS_ON('SALTSQan',myThid) ) THEN
363             DO bj = myByLo(myThid), myByHi(myThid)
364              DO bi = myBxLo(myThid), myBxHi(myThid)
365               DO K=1,Nr
366                DO j = 1,sNy
367                 DO i = 1,sNx
368                   tmpMk(i,j,K,bi,bj) =
369         &               (salt(i,j,K,bi,bj)-35)*(salt(i,j,K,bi,bj)-35)
370                 ENDDO
371                ENDDO
372               ENDDO
373              ENDDO
374             ENDDO
375             CALL DIAGNOSTICS_FILL(tmpMk,'SALTSQan',0,Nr,0,1,1,myThid)
376            ENDIF
377          
378  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
379    
380          IF ( DIAGNOSTICS_IS_ON('UVELMASS',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('UVELMASS',myThid) ) THEN
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22