/[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.19 by dimitri, Fri May 6 03:46:40 2005 UTC revision 1.20 by molod, Fri Jun 17 01:03:56 2005 UTC
# Line 40  C     == Local variables == Line 40  C     == Local variables ==
40        LOGICAL  DIAGNOSTICS_IS_ON        LOGICAL  DIAGNOSTICS_IS_ON
41        EXTERNAL DIAGNOSTICS_IS_ON        EXTERNAL DIAGNOSTICS_IS_ON
42        _RL tmpMk(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL tmpMk(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
43          _RL tmpMk1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+1,nSx,nSy)
44        _RL tmp1k(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL tmp1k(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45        INTEGER i,j,K,bi,bj        INTEGER i,j,K,bi,bj
46        INTEGER km1        INTEGER km1
# Line 328  C---+----1----+----2----+----3----+----4 Line 329  C---+----1----+----2----+----3----+----4
329           CALL DIAGNOSTICS_FILL(tmpMk,'WVELSLT ',0,Nr,0,1,1,myThid)           CALL DIAGNOSTICS_FILL(tmpMk,'WVELSLT ',0,Nr,0,1,1,myThid)
330          ENDIF          ENDIF
331                
332            IF ( DIAGNOSTICS_IS_ON('PRESSURE',myThid) ) THEN
333             DO bj = myByLo(myThid), myByHi(myThid)
334              DO bi = myBxLo(myThid), myBxHi(myThid)
335               do j = 1,sNy
336               do i = 1,sNx
337                do K = 1,Nr
338                 tmpMk1(i,j,K,bi,bj) = 0.
339                enddo
340               enddo
341               enddo
342               do j = 1,sNy
343               do i = 1,sNx
344                if(ksurfC(i,j,bi,bj).ne.0.)
345         .          tmpMk1(i,j,ksurfC(i,j,bi,bj),bi,bj) =
346         .                            (Ro_surf(i,j,bi,bj) + etaH(i,j,bi,bj))
347               enddo
348               enddo
349               do j = 1,sNy
350               do i = 1,sNx
351                do K = ksurfC(i,j,bi,bj)+1,Nr+1
352                 tmpMk1(i,j,K,bi,bj) = tmpMk1(i,j,K-1,bi,bj) -
353         .                        drF(K-1)*hfacC(i,j,K-1,bi,bj)
354                enddo
355                do K = 1,Nr
356                 tmpMk(i,j,K,bi,bj) =
357         .           (tmpMk1(i,j,K,bi,bj) + tmpMk1(i,j,K+1,bi,bj)) /2.
358                enddo
359               enddo
360               enddo
361              ENDDO
362             ENDDO
363             CALL DIAGNOSTICS_FILL(tmpMk,'PRESSURE',0,Nr,0,1,1,myThid)
364            ENDIF
365    
366  C--   fill momentum state-var diagnostics: end  C--   fill momentum state-var diagnostics: end
367        ENDIF        ENDIF
368    
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22