pkg/my82/my82_init_varia.F

08be60903a Mart*0001 #include "MY82_OPTIONS.h"
                0002
60c545bbe1 Mart*0003       SUBROUTINE MY82_INIT_VARIA( myThid )
5e48dccc42 Jean*0004 C     *==========================================================*
60c545bbe1 Mart*0005 C     | SUBROUTINE MY82_INIT_VARIA                               |
08be60903a Mart*0006 C     | o Routine to initialize MY82 parameters and variables.   |
5e48dccc42 Jean*0007 C     *==========================================================*
08be60903a Mart*0008 C     | Initialize MY92 parameters and variables.                |
5e48dccc42 Jean*0009 C     *==========================================================*
08be60903a Mart*0010       IMPLICIT NONE
                0011
                0012 C     === Global variables ===
                0013 #include "SIZE.h"
                0014 #include "EEPARAMS.h"
                0015 #include "PARAMS.h"
                0016 #include "GRID.h"
                0017 #include "MY82.h"
                0018
                0019 C     === Routine arguments ===
                0020 C     myThid -  Number of this instance of MY_INIT
                0021       INTEGER myThid
                0022
                0023 #ifdef ALLOW_MY82
                0024
                0025 C     === Local variables ===
                0026 C     i,j,k,bi,bj - Loop counters
                0027
                0028       INTEGER i, j, k, bi, bj
                0029       _RL gam1, gam2
                0030
                0031 C     magic parameters( M. Satoh, p. 314)
fb62f539dc Jean*0032       gam1 = 1. _d 0/3. _d 0 - 2. _d 0 *A1/B1
                0033       gam2 = (B2+6. _d 0*A1)/B1
                0034       alpha1 = 3. _d 0 *A2*gam1
                0035       alpha2 = 3. _d 0 *A2*(gam1+gam2)
08be60903a Mart*0036       beta1  = A1*B1*(gam1-C1)
fb62f539dc Jean*0037       beta2  = A1*(B1*(gam1-C1) + 6. _d 0 *A1 + 3. _d 0 *A2)
08be60903a Mart*0038       beta3  = A2*B1*gam1
fb62f539dc Jean*0039       beta4  = A2*(B1*(gam1+gam2) - 3. _d 0 *A1)
08be60903a Mart*0040
                0041 C-----------------------------------------------------------------------
                0042 C Initialize MY82 variables MYviscAr, MYdiffKrT
                0043 C-----------------------------------------------------------------------
                0044
5e48dccc42 Jean*0045       DO bj = myByLo(myThid), myByHi(myThid)
                0046        DO bi = myBxLo(myThid), myBxHi(myThid)
                0047         DO k = 1, Nr
                0048          DO j = 1-OLy, sNy+OLy
                0049           DO i = 1-OLx, sNx+OLx
8fef4b1a57 Mart*0050            MYhbl    (i,j  ,bi,bj) = 0. _d 0
5e48dccc42 Jean*0051            MYviscAr (i,j,k,bi,bj) = viscArNr(k)
78524d1402 Jean*0052            MYdiffKr (i,j,k,bi,bj) = diffKrNrS(k)
5e48dccc42 Jean*0053           ENDDO
                0054          ENDDO
                0055         ENDDO
                0056        ENDDO
                0057       ENDDO
08be60903a Mart*0058
                0059 #endif /* ALLOW_MY82 */
                0060
5e48dccc42 Jean*0061       RETURN
                0062       END