model/src/calc_ivdc.F

d4701cb6da Alis*0001 #include "CPP_OPTIONS.h"
9366854e02 Chri*0002 CBOP
                0003 C     !ROUTINE: CALC_IVDC
                0004 C     !INTERFACE:
fa49f16b06 Jean*0005       SUBROUTINE CALC_IVDC(
                0006      I       bi, bj, iMin, iMax, jMin, jMax, K,
da9b4e6c91 Mart*0007      I       sigmaR,
fa49f16b06 Jean*0008      I       myTime,myIter, myThid )
9366854e02 Chri*0009 C     !DESCRIPTION: \bv
                0010 C     *==========================================================*
                0011 C     | SUBROUTINE CALC_IVDC
                0012 C     | o Calculates Implicit Vertical Diffusivity for Convection
                0013 C     \==========================================================*
                0014 C     \ev
                0015 C     !USES:
d4701cb6da Alis*0016       IMPLICIT NONE
                0017 C     == Global data ==
                0018 #include "SIZE.h"
                0019 #include "EEPARAMS.h"
                0020 #include "PARAMS.h"
                0021 #include "DYNVARS.h"
                0022 #include "GRID.h"
                0023
9366854e02 Chri*0024 C     !INPUT/OUTPUT PARAMETERS:
d4701cb6da Alis*0025 C     == Routine arguments ==
9366854e02 Chri*0026 C     bi,bj,iMin,iMax,jMin,jMax,K :: Loop counters
                0027 C     rhoKm1 :: rho in layer above
                0028 C     rhoKp1 :: rho in layer below
                0029 C     myTime :: Current time in simulation
                0030 C     myIter :: Current iteration in simulation
da9b4e6c91 Mart*0031 C     myThid :: Thread number of this instance of S/R CALC_IVDC
d4701cb6da Alis*0032       INTEGER bi,bj,iMin,iMax,jMin,jMax,K
da9b4e6c91 Mart*0033       _RL sigmaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
d4701cb6da Alis*0034       _RL myTime
                0035       INTEGER myIter
                0036       INTEGER myThid
                0037
9366854e02 Chri*0038 C     !LOCAL VARIABLES:
d4701cb6da Alis*0039 C     == Local variables ==
9366854e02 Chri*0040 C     i,j :: Loop counters
d4701cb6da Alis*0041       INTEGER i,j
9366854e02 Chri*0042 CEOP
d4701cb6da Alis*0043
                0044 C--    Where statically unstable, mix the heat and salt
                0045        DO j=jMin,jmax
                0046         DO i=iMin,imax
da9b4e6c91 Mart*0047          IF ( -sigmaR(i,j,k)*gravitySign.GT.0. ) THEN
fa49f16b06 Jean*0048           IVDConvCount(i,j,k,bi,bj) = 1. _d 0
                0049          ELSE
                0050           IVDConvCount(i,j,k,bi,bj) = 0. _d 0
d4701cb6da Alis*0051          ENDIF
                0052         ENDDO
                0053        ENDDO
                0054
                0055       RETURN
                0056       END