pkg/monitor/mon_advcflw2.F

91672e10e3 Alis*0001 #include "MONITOR_OPTIONS.h"
141e421ae7 Jean*0002
7633b97660 Ed H*0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP 0
                0005 C     !ROUTINE: MON_ADVCFLW2
                0006
                0007 C     !INTERFACE:
141e421ae7 Jean*0008       SUBROUTINE MON_ADVCFLW2(
7633b97660 Ed H*0009      I     label, W, rHFac, rDrF, dT,
                0010      I     myThid )
                0011
                0012 C     !DESCRIPTION:
                0013 C     Calculates maximum CFL number in the vertical relevant for tracer
                0014 C     Adv. Pb. with Partial Cell.
141e421ae7 Jean*0015
7633b97660 Ed H*0016 C     !USES:
                0017       IMPLICIT NONE
141e421ae7 Jean*0018 #include "SIZE.h"
                0019 #include "EEPARAMS.h"
                0020 #include "MONITOR.h"
                0021
7633b97660 Ed H*0022 C     !INPUT PARAMETERS:
141e421ae7 Jean*0023       CHARACTER*(*) label
                0024       _RL     W(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0025       _RS rHFac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0026       _RS rDrF(Nr)
                0027       _RL dT
                0028       INTEGER myThid
7633b97660 Ed H*0029 CEOP
141e421ae7 Jean*0030
7633b97660 Ed H*0031 C     !LOCAL VARIABLES:
141e421ae7 Jean*0032       INTEGER bi,bj,I,J,K
                0033       _RL tmpVal,theMax
                0034
                0035       theMax=0.
                0036
                0037       DO bj=myByLo(myThid),myByHi(myThid)
                0038        DO bi=myBxLo(myThid),myBxHi(myThid)
                0039         DO K=2,Nr
                0040          DO J=1,sNy
                0041           DO I=1,sNx
                0042            tmpVal=abs(W(I,J,K,bi,bj))*dT
                0043      &           *max(  rDrF(K)*rHfac(I,J,K,bi,bj),
                0044      &                rDrF(K-1)*rHfac(I,J,K-1,bi,bj) )
                0045            theMax=max(theMax,tmpVal)
                0046           ENDDO
                0047          ENDDO
                0048         ENDDO
                0049        ENDDO
                0050       ENDDO
7163a40534 Jean*0051       _GLOBAL_MAX_RL(theMax,myThid)
141e421ae7 Jean*0052
                0053       CALL MON_OUT_RL( label, theMax, mon_foot_max ,myThid)
                0054
                0055       RETURN
                0056       END
7633b97660 Ed H*0057
                0058 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|