/[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.26 by jmc, Thu Aug 25 20:51:38 2005 UTC revision 1.28 by molod, Wed Nov 16 16:55:21 2005 UTC
# Line 270  C---+----1----+----2----+----3----+----4 Line 270  C---+----1----+----2----+----3----+----4
270           CALL DIAGNOSTICS_FILL(tmpMk,'WVELSLT ',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'WVELSLT ',0,Nr,0,1,1,myThid)
271          ENDIF          ENDIF
272                
273            IF ( DIAGNOSTICS_IS_ON('UVELPHI ',myThid) ) THEN
274             DO bj = myByLo(myThid), myByHi(myThid)
275              DO bi = myBxLo(myThid), myBxHi(myThid)
276               DO K=1,Nr
277                DO j = 1,sNy
278                 DO i = 1,sNx+1
279                   tmpMk(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*0.5 _d 0
280         &            *(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i-1,j,K,bi,bj))
281                 ENDDO
282                ENDDO
283               ENDDO
284              ENDDO
285             ENDDO
286             CALL DIAGNOSTICS_FILL(tmpMk,'UVELPHI ',0,Nr,0,1,1,myThid)
287            ENDIF
288          
289            IF ( DIAGNOSTICS_IS_ON('VVELPHI ',myThid) ) THEN
290             DO bj = myByLo(myThid), myByHi(myThid)
291              DO bi = myBxLo(myThid), myBxHi(myThid)
292               DO K=1,Nr
293                DO j = 1,sNy+1
294                 DO i = 1,sNx
295                   tmpMk(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*0.5 _d 0
296         &            *(totPhiHyd(i,j,K,bi,bj)+totPhiHyd(i,j-1,K,bi,bj))
297                 ENDDO
298                ENDDO
299               ENDDO
300              ENDDO
301             ENDDO
302             CALL DIAGNOSTICS_FILL(tmpMk,'VVELPHI ',0,Nr,0,1,1,myThid)
303            ENDIF
304          
305          IF ( DIAGNOSTICS_IS_ON('PRESSURE',myThid) ) THEN          IF ( DIAGNOSTICS_IS_ON('PRESSURE',myThid) ) THEN
306           DO bj = myByLo(myThid), myByHi(myThid)           DO bj = myByLo(myThid), myByHi(myThid)
307            DO bi = myBxLo(myThid), myBxHi(myThid)            DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 317  C--   fill tracer state-var diagnostics: Line 349  C--   fill tracer state-var diagnostics:
349             call qsat(dummy1,dummy2,dummy3,dummy4,.false.)             call qsat(dummy1,dummy2,dummy3,dummy4,.false.)
350             tmpMk(i,j,K,bi,bj) = hfacC(i,j,K,bi,bj) *             tmpMk(i,j,K,bi,bj) = hfacC(i,j,K,bi,bj) *
351       .                    salt(i,j,k,bi,bj) * 100. / dummy3       .                    salt(i,j,k,bi,bj) * 100. / dummy3
            if(k.lt.4) print *,' Temp ',dummy1,' Pres ',dummy2,' q ',  
      .      salt(i,j,k,bi,bj),' qstar ',dummy3,' rh ',tmpMk(i,j,K,bi,bj)  
352            enddo            enddo
353           enddo           enddo
354          enddo          enddo
Legend:
Removed from v.1.26  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22