model/src/ini_mixing.F

8fe5594931 Patr*0001 #include "CPP_OPTIONS.h"
4f7b6674b6 Patr*0002
9366854e02 Chri*0003 CBOP
                0004 C     !ROUTINE: INI_MIXING
                0005 C     !INTERFACE:
8fe5594931 Patr*0006       SUBROUTINE INI_MIXING( myThid )
9366854e02 Chri*0007 C     !DESCRIPTION: \bv
                0008 C     *==========================================================*
67be4914b2 Jean*0009 C     | SUBROUTINE INI_MIXING
                0010 C     | o Initialise diffusivity to default constant value.
9366854e02 Chri*0011 C     *==========================================================*
                0012 C     \ev
8fe5594931 Patr*0013
9366854e02 Chri*0014 C     !USES:
                0015       IMPLICIT NONE
8fe5594931 Patr*0016 C     === Global variables ===
                0017 #include "SIZE.h"
                0018 #include "EEPARAMS.h"
                0019 #include "PARAMS.h"
                0020 #include "GRID.h"
                0021 #include "DYNVARS.h"
                0022
9366854e02 Chri*0023 C     !INPUT/OUTPUT PARAMETERS:
4125e77f3d Jean*0024 C     myThid :: my Thread Id number
8fe5594931 Patr*0025       INTEGER myThid
                0026
9366854e02 Chri*0027 C     !LOCAL VARIABLES:
4125e77f3d Jean*0028 #if (defined ALLOW_3D_DIFFKR) || (defined ALLOW_BL79_LAT_VARY)
239f048bdf Jean*0029 C     bi,bj  :: tile indices
                0030 C     i,j,k  :: Loop counters
8fe5594931 Patr*0031       INTEGER bi, bj
4125e77f3d Jean*0032       INTEGER i, j
                0033 #endif
                0034 #ifdef ALLOW_3D_DIFFKR
                0035       INTEGER k
                0036 #endif
9366854e02 Chri*0037 CEOP
8fe5594931 Patr*0038
94a8024bbe Jean*0039 #ifdef ALLOW_3D_DIFFKR
8fe5594931 Patr*0040        DO bj = myByLo(myThid), myByHi(myThid)
                0041         DO bi = myBxLo(myThid), myBxHi(myThid)
8fc5af62b8 Jean*0042          DO k=1,Nr
517dbdc414 Jean*0043           DO j = 1-OLy, sNy+OLy
                0044            DO i = 1-OLx, sNx+OLx
8fc5af62b8 Jean*0045             diffKr(i,j,k,bi,bj) = diffKrNrS(k)
8fe5594931 Patr*0046            ENDDO
                0047           ENDDO
                0048          ENDDO
                0049         ENDDO
                0050        ENDDO
94a8024bbe Jean*0051        IF ( diffKrFile .NE. ' ' ) THEN
                0052          CALL READ_FLD_XYZ_RL( diffKrFile, ' ', diffKr, 0, myThid )
                0053          _EXCH_XYZ_RL( diffKr, myThid )
                0054        ENDIF
                0055 #endif /* ALLOW_3D_DIFFKR */
8fe5594931 Patr*0056
e40c34e398 Dimi*0057 #ifdef ALLOW_BL79_LAT_VARY
                0058        DO bj = myByLo(myThid), myByHi(myThid)
                0059         DO bi = myBxLo(myThid), myBxHi(myThid)
517dbdc414 Jean*0060          DO j = 1-OLy, sNy+OLy
                0061           DO i = 1-OLx, sNx+OLx
e40c34e398 Dimi*0062            BL79LatArray(i,j,bi,bj) = 0.0 _d 0
86e9b00011 Jean*0063            IF ( ABS(YC(i,j,bi,bj)) .LT. BL79LatVary )
4125e77f3d Jean*0064      &       BL79LatArray(i,j,bi,bj) = halfRL
                0065      &         *( 1. _d 0 + COS(YC(i,j,bi,bj)*PI/BL79LatVary) )
e40c34e398 Dimi*0066           ENDDO
                0067          ENDDO
                0068         ENDDO
                0069        ENDDO
                0070 #endif
                0071
8fe5594931 Patr*0072       RETURN
                0073       END