/[MITgcm]/MITgcm/pkg/my82/MY82.h
ViewVC logotype

Contents of /MITgcm/pkg/my82/MY82.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.1 - (show annotations) (download)
Thu Sep 2 09:11:54 2004 UTC (19 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint55c_post, checkpoint54e_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57e_post, checkpoint55h_post, checkpoint57g_pre, checkpoint55b_post, checkpoint56c_post, checkpoint55, checkpoint57f_pre, checkpoint57a_post, checkpoint54f_post, checkpoint55g_post, checkpoint55f_post, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, eckpoint57e_pre, checkpoint57h_done, checkpoint57f_post, checkpoint57c_post, checkpoint55e_post, checkpoint55a_post, checkpoint57h_pre, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post
File MIME type: text/plain 
o add two new packages
  - pp81 (Packanowski and Philander, 1981), Richardson number and
    stratification dependent mixing
  - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
1 C $Header: $
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.92 )
40 PARAMETER( A2 = 0.74 )
41 PARAMETER( B1 = 16.6 )
42 PARAMETER( B2 = 10.1 )
43 PARAMETER( C1 = 0.08 )
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 */
  ViewVC Help
Powered by ViewVC 1.1.22