1 |
C $Header: /u/gcmpack/MITgcm/pkg/my82/MY82.h,v 1.1 2004/09/02 09:11:54 mlosch Exp $ |
2 |
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 |
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 |
|
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 */ |