1 |
mlosch |
1.2 |
C $Header: /u/gcmpack/MITgcm/pkg/my82/MY82.h,v 1.1 2004/09/02 09:11:54 mlosch Exp $ |
2 |
mlosch |
1.1 |
C $Name: $ |
3 |
|
|
|
4 |
|
|
#ifdef ALLOW_MY82 |
5 |
|
|
|
6 |
|
|
CBOP |
7 |
|
|
C !ROUTINE: MY82.h |
8 |
|
|
|
9 |
|
|
C !DESCRIPTION: \bv |
10 |
|
|
C /==========================================================\ |
11 |
|
|
C | MY82.h | |
12 |
|
|
C | o Basic header for Pacanowski and Philander (1981) | |
13 |
|
|
C | vertical mixing parameterization. | |
14 |
|
|
C | Contains all MY82 field declarations. | |
15 |
|
|
C \==========================================================/ |
16 |
|
|
|
17 |
|
|
C----------------------------------------------------------------------- |
18 |
|
|
C |
19 |
|
|
C Constants that can be set in data.pp |
20 |
|
|
C MYviscMax, MYdiffMax - maximum allowed viscosity, diffusivity |
21 |
|
|
C MYhblScale - scale the boundary length scale |
22 |
|
|
C RiMax - Maximum of Richardson number |
23 |
|
|
C MYdumpFreq, MYtaveFreq - analogue to dumpFreq and taveFreq (=default) |
24 |
|
|
C MYmixingMaps - if true, include MY diagnostic maps in STDOUT |
25 |
|
|
C MYwriteState - if true, write MY state to file |
26 |
|
|
C |
27 |
|
|
C Time varying parameters computed by subroutine pp_calc |
28 |
|
|
C MYviscAr - Vertical eddy viscosity coefficient (m^2/s) |
29 |
|
|
C MYdiffKr - Vertical diffusion coefficient for heat, |
30 |
|
|
C salt and tracers (m^2/s) |
31 |
|
|
C |
32 |
|
|
C----------------------------------------------------------------------- |
33 |
|
|
C \ev |
34 |
|
|
CEOP |
35 |
|
|
|
36 |
|
|
C Magic parameters of Mellor&Yamada(1982): |
37 |
|
|
C (M. Satoh, p.315) |
38 |
|
|
_RL A1, A2, B1, B2, C1 |
39 |
mlosch |
1.2 |
PARAMETER( A1 = 0.92D0 ) |
40 |
|
|
PARAMETER( A2 = 0.74D0 ) |
41 |
|
|
PARAMETER( B1 = 16.6D0 ) |
42 |
|
|
PARAMETER( B2 = 10.1D0 ) |
43 |
|
|
PARAMETER( C1 = 0.08D0 ) |
44 |
mlosch |
1.1 |
|
45 |
|
|
_RL alpha1, alpha2 |
46 |
|
|
_RL beta1, beta2, beta3, beta4 |
47 |
|
|
_RL RiMax |
48 |
|
|
_RL MYhblScale |
49 |
|
|
_RL MYviscMax, MYdiffMax |
50 |
|
|
_RL MYdumpFreq, MYtaveFreq |
51 |
|
|
COMMON /MY_PARMS_R/ |
52 |
|
|
& alpha1, alpha2, beta1, beta2, beta3, beta4, |
53 |
|
|
& RiMax, MYhblScale, |
54 |
|
|
& MYviscMax, MYdiffMax, |
55 |
|
|
& MYdumpFreq, MYtaveFreq |
56 |
|
|
|
57 |
|
|
_RL MYhbl (1-OLx:sNx+OLx,1-OLy:sNy+OLy, nSx,nSy) |
58 |
|
|
_RL MYviscAr (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
59 |
|
|
_RL MYdiffKr (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) |
60 |
|
|
COMMON /MY_FIELDS/ MYhbl, MYviscAr, MYdiffKr |
61 |
|
|
|
62 |
|
|
LOGICAL MYisOn, MYmixingMaps, MYwriteState |
63 |
|
|
COMMON /MY_PARMS_L/ |
64 |
|
|
& MYisOn, MYmixingMaps, MYwriteState |
65 |
|
|
|
66 |
|
|
#endif /* ALLOW_MY82 */ |