pkg/my82/my82_init_varia.F

C $Header: $
C $Name:  $

#include "MY82_OPTIONS.h"

      SUBROUTINE MY82_INIT_VARIA( myThid )
C     /==========================================================\
C     | SUBROUTINE MY82_INIT_VARIA                               |
C     | o Routine to initialize MY82 parameters and variables.   |
C     |==========================================================|
C     | Initialize MY92 parameters and variables.                |
C     \==========================================================/
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "MY82.h"

C     === Routine arguments ===
C     myThid -  Number of this instance of MY_INIT
      INTEGER myThid

#ifdef ALLOW_MY82

C     === Local variables ===
C     i,j,k,bi,bj - Loop counters

      INTEGER i, j, k, bi, bj
      _RL gam1, gam2

C     magic parameters( M. Satoh, p. 314)
      gam1 = 1./3.-2.*A1/B1
      gam2 = (B2+6.*A1)/B1
      alpha1 = 3.*A2*gam1
      alpha2 = 3.*A2*(gam1+gam2)
      beta1  = A1*B1*(gam1-C1)
      beta2  = A1*(B1*(gam1-C1) + 6.*A1 + 3.*A2)
      beta3  = A2*B1*gam1
      beta4  = A2*(B1*(gam1+gam2) - 3.*A1)

C-----------------------------------------------------------------------
C Initialize MY82 variables MYviscAr, MYdiffKrT
C-----------------------------------------------------------------------

      do bj = myByLo(myThid), myByHi(myThid)
       do bi = myBxLo(myThid), myBxHi(myThid)
        do k = 1, Nr
         do j = 1-OLy, sNy+OLy
          do i = 1-OLx, sNx+OLx
           MYhbl    (i,j  ,bi,bj) = 0. _d 0
           MYviscAr (i,j,k,bi,bj) = viscAr
           MYdiffKr (i,j,k,bi,bj) = diffKrNrT(k)
          end do
         end do
        end do
       end do
      end do

#endif /* ALLOW_MY82 */

      return
      end
1	mlosch	1.1	C $Header: $
2			C $Name: $
3
4			#include "MY82_OPTIONS.h"
5
6			SUBROUTINE MY82_INIT_VARIA( myThid )
7			C /==========================================================\
8			C \| SUBROUTINE MY82_INIT_VARIA \|
9			C \| o Routine to initialize MY82 parameters and variables. \|
10			C \|==========================================================\|
11			C \| Initialize MY92 parameters and variables. \|
12			C \==========================================================/
13			IMPLICIT NONE
14
15			C === Global variables ===
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19			#include "GRID.h"
20			#include "MY82.h"
21
22			C === Routine arguments ===
23			C myThid - Number of this instance of MY_INIT
24			INTEGER myThid
25
26			#ifdef ALLOW_MY82
27
28			C === Local variables ===
29			C i,j,k,bi,bj - Loop counters
30
31			INTEGER i, j, k, bi, bj
32			_RL gam1, gam2
33
34			C magic parameters( M. Satoh, p. 314)
35			gam1 = 1./3.-2.*A1/B1
36			gam2 = (B2+6.*A1)/B1
37			alpha1 = 3.A2gam1
38			alpha2 = 3.A2(gam1+gam2)
39			beta1 = A1B1(gam1-C1)
40			beta2 = A1(B1(gam1-C1) + 6.A1 + 3.A2)
41			beta3 = A2B1gam1
42			beta4 = A2(B1(gam1+gam2) - 3.*A1)
43
44			C-----------------------------------------------------------------------
45			C Initialize MY82 variables MYviscAr, MYdiffKrT
46			C-----------------------------------------------------------------------
47
48			do bj = myByLo(myThid), myByHi(myThid)
49			do bi = myBxLo(myThid), myBxHi(myThid)
50			do k = 1, Nr
51			do j = 1-OLy, sNy+OLy
52			do i = 1-OLx, sNx+OLx
53			MYhbl (i,j ,bi,bj) = 0. _d 0
54			MYviscAr (i,j,k,bi,bj) = viscAr
55			MYdiffKr (i,j,k,bi,bj) = diffKrNrT(k)
56			end do
57			end do
58			end do
59			end do
60			end do
61
62			#endif /* ALLOW_MY82 */
63
64			return
65			end