/[MITgcm]/MITgcm/pkg/kpp/KPP_PARAMS.h
ViewVC logotype

Contents of /MITgcm/pkg/kpp/KPP_PARAMS.h

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

Revision 1.7 - (show annotations) (download)
Sun Feb 4 14:38:50 2001 UTC (23 years, 3 months ago) by cnh
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint51k_post, checkpoint47e_post, checkpoint52l_pre, checkpoint44e_post, hrcube4, hrcube5, checkpoint46l_post, checkpoint46g_pre, checkpoint47c_post, release1_p13_pre, checkpoint50c_post, checkpoint46f_post, checkpoint52d_pre, checkpoint48e_post, checkpoint50c_pre, checkpoint44f_post, checkpoint46b_post, checkpoint52j_pre, checkpoint43a-release1mods, checkpoint51o_pre, ecco_c50_e32, ecco_c50_e33, ecco_c50_e30, ecco_c50_e31, release1_p13, checkpoint40pre3, checkpoint40pre1, checkpoint51l_post, checkpoint40pre7, checkpoint40pre6, checkpoint48i_post, checkpoint40pre9, checkpoint40pre8, checkpoint46l_pre, checkpoint52l_post, checkpoint52k_post, chkpt44d_post, checkpoint54, checkpoint51, checkpoint50, checkpoint53, checkpoint52, release1_p8, release1_p9, checkpoint50d_post, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, release1_p7, checkpoint52f_post, checkpoint50b_pre, checkpoint44e_pre, checkpoint51f_post, release1_b1, ecco_c51_e34a, checkpoint48b_post, checkpoint43, checkpoint51d_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint51t_post, checkpoint38, checkpoint51n_post, release1_chkpt44d_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint47a_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint52e_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint47d_post, icebear5, icebear4, icebear3, icebear2, checkpoint53d_post, checkpoint46d_pre, checkpoint40pre2, checkpoint48d_post, release1-branch_tutorials, checkpoint48f_post, checkpoint45d_post, checkpoint52b_pre, checkpoint54b_post, checkpoint46j_pre, checkpoint51l_pre, checkpoint52m_post, chkpt44a_post, checkpoint44h_pre, checkpoint48h_post, checkpoint51q_post, checkpoint40pre4, ecco_c50_e29, checkpoint51b_pre, checkpoint46a_post, checkpoint47g_post, checkpoint52b_post, checkpoint52c_post, checkpoint46j_post, checkpoint51h_pre, checkpoint46k_post, ecco_c50_e28, chkpt44c_pre, checkpoint48a_post, checkpoint45a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint52f_pre, ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, release1_p12, release1_p10, release1_p11, release1_p16, release1_p17, release1_p14, release1_p15, checkpoint47j_post, pre38tag1, checkpoint54a_pre, ecco_c50_e33a, checkpoint53c_post, branch-exfmods-tag, checkpoint44g_post, branchpoint-genmake2, checkpoint54a_post, checkpoint46e_pre, checkpoint51r_post, checkpoint48c_post, checkpoint45b_post, checkpoint46b_pre, checkpoint51i_post, release1-branch-end, c37_adj, release1_final_v1, checkpoint51b_post, checkpoint51c_post, checkpoint46c_pre, checkpoint53a_post, checkpoint46, checkpoint47b_post, checkpoint44b_post, ecco_c51_e34, checkpoint46h_pre, checkpoint52d_post, checkpoint53g_post, checkpoint46m_post, checkpoint46a_pre, checkpoint50g_post, checkpoint45c_post, ecco_ice2, ecco_ice1, checkpoint44h_post, pre38-close, checkpoint46g_post, release1_p12_pre, checkpoint39, checkpoint52a_pre, checkpoint37, checkpoint36, checkpoint35, ecco_c44_e22, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, ecco_c44_e25, checkpoint51i_pre, checkpoint52h_pre, checkpoint53f_post, checkpoint40pre5, checkpoint52j_post, checkpoint47f_post, checkpoint50e_post, chkpt44a_pre, checkpoint46i_post, ecco_c44_e23, ecco_c44_e20, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, checkpoint46c_post, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, branch-netcdf, checkpoint50d_pre, checkpoint52n_post, checkpoint53b_pre, checkpoint46e_post, release1_beta1, checkpoint51e_post, checkpoint44b_pre, checkpoint42, checkpoint40, checkpoint41, checkpoint47, checkpoint44, checkpoint45, checkpoint48, checkpoint49, checkpoint46h_post, checkpoint51o_post, checkpoint51f_pre, chkpt44c_post, checkpoint48g_post, checkpoint53b_post, checkpoint47h_post, checkpoint52a_post, checkpoint44f_pre, checkpoint51g_post, ecco_c52_e35, checkpoint46d_post, checkpoint50b_post, release1-branch_branchpoint, checkpoint51m_post, checkpoint53d_pre, checkpoint51a_post, checkpoint51p_post, checkpoint51u_post
Branch point for: c24_e25_ice, branch-exfmods-curt, release1_final, release1-branch, branch-genmake2, release1, branch-nonh, tg2-branch, ecco-branch, release1_50yr, netcdf-sm0, icebear, checkpoint51n_branch, pre38, release1_coupled
Changes since 1.6: +2 -1 lines
File MIME type: text/plain 
Made sure each .F and .h file had
the CVS keywords Header and Name at its start.
Most had header but very few currently have Name, so
lots of changes!
1 C $Header: /u/gcmpack/models/MITgcmUV/pkg/kpp/KPP_PARAMS.h,v 1.6 2001/01/29 20:09:23 heimbach Exp $
2 C $Name: $
3
4 C /==========================================================\
5 C | KPP_PARAMS.h |
6 C | o Basic parameter header for KPP vertical mixing |
7 C | parameterization. These parameters are initialized by |
8 C | and/or read in from data.kpp file. |
9 C \==========================================================/
10
11 C Parameters used in kpp routine arguments (needed for compilation
12 C of kpp routines even if ALLOW_KPP is not defined)
13 C mdiff - number of diffusivities for local arrays
14 C Nrm1, Nrp1, Nrp2 - number of vertical levels
15 C imt - array dimension for local arrays
16 C ibot, itop, jbot, jtop - array dimension indices
17
18 integer mdiff, Nrm1, Nrp1, Nrp2
19 integer imt, ibot, itop, jbot, jtop
20 parameter( mdiff = 3 )
21 parameter( Nrm1 = Nr-1 )
22 parameter( Nrp1 = Nr+1 )
23 parameter( Nrp2 = Nr+2 )
24 #ifdef FRUGAL_KPP
25 parameter( imt=(sNx+2)*(sNy+2) )
26 parameter( ibot=0, itop=sNx+1, jbot=0, jtop=sNy+1 )
27 #else
28 parameter( imt=(sNx+2*OLx)*(sNy+2*OLy) )
29 parameter( ibot=1-OLx, itop=sNx+OLx, jbot=1-OLy, jtop=sNy+OLy )
30 #endif
31
32 #ifdef ALLOW_KPP
33
34 C Time invariant parameters initialized by subroutine kmixinit
35 C nzmax (nx,ny) - Maximum number of wet levels in each column
36 C pMask - Mask relating to Pressure/Tracer point grid.
37 C 0. if P point is on land.
38 C 1. if P point is in water.
39 C zgrid (0:Nr+1) - vertical levels of tracers (<=0) (m)
40 C hwide (0:Nr+1) - layer thicknesses (>=0) (m)
41 C kpp_freq - Re-computation frequency for KPP parameters (s)
42 C kpp_dumpFreq - KPP dump frequency. (s)
43 C kpp_taveFreq - KPP time-averaging frequency. (s)
44
45
46 INTEGER nzmax ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, nSx, nSy )
47 _KPP_RL pMask ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr, nSx, nSy )
48 _KPP_RL zgrid ( 0:Nr+1 )
49 _KPP_RL hwide ( 0:Nr+1 )
50 _RL kpp_freq
51 _RL kpp_dumpFreq
52 _RL kpp_taveFreq
53
54 COMMON /kpp_i/ nzmax
55
56 COMMON /kpp_r1/ pMask, zgrid, hwide
57
58 COMMON /kpp_r2/ kpp_freq, kpp_dumpFreq, kpp_taveFreq
59
60
61 C-----------------------------------------------------------------------
62 C
63 C KPP flags and min/max permitted values for mixing parameters
64 c
65 C KPPmixingMaps - if true, include KPP diagnostic maps in STDOUT
66 C KPPwriteState - if true, write KPP state to file
67 C minKPPhbl - KPPhbl minimum value (m)
68 C
69 C-----------------------------------------------------------------------
70
71 LOGICAL KPPmixingMaps, KPPwriteState
72
73 COMMON /KPP_PARM_L/
74 & KPPmixingMaps, KPPwriteState
75
76 _KPP_RL minKPPhbl
77
78 COMMON /KPP_PARM_R/
79 & minKPPhbl
80
81 c====================== file "kmixcom.h" =======================
82 c
83 c-----------------------------------------------------------------------
84 c Define various parameters and common blocks for KPP vertical-
85 c mixing scheme; used in "kppmix.F" subroutines.
86 c Constants are set in subroutine "ini_parms".
87 c-----------------------------------------------------------------------
88 c
89 c-----------------------------------------------------------------------
90 c parameters for several subroutines
91 c
92 c epsln = 1.0e-20
93 c phepsi = 1.0e-10
94 c epsilon = nondimensional extent of the surface layer = 0.1
95 c vonk = von Karmans constant = 0.4
96 c dB_dz = maximum dB/dz in mixed layer hMix = 5.2e-5 s^-2
97 c conc1,conam,concm,conc2,zetam,conas,concs,conc3,zetas
98 c = scalar coefficients
99 c-----------------------------------------------------------------------
100
101 _KPP_RL epsln,phepsi,epsilon,vonk,dB_dz,
102 $ conc1,
103 $ conam,concm,conc2,zetam,
104 $ conas,concs,conc3,zetas
105
106 common /kmixcom/ epsln,phepsi,epsilon,vonk,dB_dz,
107 $ conc1,
108 $ conam,concm,conc2,zetam,
109 $ conas,concs,conc3,zetas
110
111 c-----------------------------------------------------------------------
112 c parameters for subroutine "bldepth"
113 c
114 c
115 c to compute depth of boundary layer:
116 c
117 c Ricr = critical bulk Richardson Number = 0.3
118 c cekman = coefficient for ekman depth = 0.7
119 c cmonob = coefficient for Monin-Obukhov depth = 1.0
120 c concv = ratio of interior buoyancy frequency to
121 c buoyancy frequency at entrainment depth = 1.8
122 c hbf = fraction of bounadry layer depth to
123 c which absorbed solar radiation
124 c contributes to surface buoyancy forcing = 1.0
125 c Vtc = non-dimensional coefficient for velocity
126 c scale of turbulant velocity shear
127 c (=function of concv,concs,epsilon,vonk,Ricr)
128 c-----------------------------------------------------------------------
129
130 _KPP_RL Ricr,cekman,cmonob,concv,Vtc
131 _RL hbf
132
133 common /kpp_bldepth1/ Ricr,cekman,cmonob,concv,Vtc
134 common /kpp_bldepth2/ hbf
135
136 c-----------------------------------------------------------------------
137 c parameters and common arrays for subroutines "kmixinit"
138 c and "wscale"
139 c
140 c
141 c to compute turbulent velocity scales:
142 c
143 c nni = number of values for zehat in the look up table
144 c nnj = number of values for ustar in the look up table
145 c
146 c wmt = lookup table for wm, the turbulent velocity scale
147 c for momentum
148 c wst = lookup table for ws, the turbulent velocity scale
149 c for scalars
150 c deltaz = delta zehat in table
151 c deltau = delta ustar in table
152 c zmin = minimum limit for zehat in table (m3/s3)
153 c zmax = maximum limit for zehat in table
154 c umin = minimum limit for ustar in table (m/s)
155 c umax = maximum limit for ustar in table
156 c-----------------------------------------------------------------------
157
158 integer nni , nnj
159 parameter (nni = 890, nnj = 480)
160
161 _KPP_RL wmt(0:nni+1,0:nnj+1), wst(0:nni+1,0:nnj+1)
162 _KPP_RL deltaz,deltau,zmin,zmax,umin,umax
163 common /kmixcws/ wmt, wst
164 $ , deltaz,deltau,zmin,zmax,umin,umax
165
166 c-----------------------------------------------------------------------
167 c parameters for subroutine "ri_iwmix"
168 c
169 c
170 c to compute vertical mixing coefficients below boundary layer:
171 c
172 c num_v_smooth_Ri = number of times Ri is vertically smoothed
173 c num_v_smooth_BV, num_z_smooth_sh, and num_m_smooth_sh are dummy
174 c variables kept for backward compatibility of the data file
175 c Riinfty = local Richardson Number limit for shear instability = 0.7
176 c BVSQcon = Brunt-Vaisala squared (1/s^2)
177 c difm0 = viscosity max due to shear instability (m^2/s)
178 c difs0 = tracer diffusivity .. (m^2/s)
179 c dift0 = heat diffusivity .. (m^2/s)
180 c difmcon = viscosity due to convective instability (m^2/s)
181 c difscon = tracer diffusivity .. (m^2/s)
182 c diftcon = heat diffusivity .. (m^2/s)
183 c-----------------------------------------------------------------------
184
185 INTEGER num_v_smooth_Ri, num_v_smooth_BV
186 INTEGER num_z_smooth_sh, num_m_smooth_sh
187 _KPP_RL Riinfty, BVSQcon
188 _KPP_RL difm0 , difs0 , dift0
189 _KPP_RL difmcon, difscon, diftcon
190
191 COMMON /kmixcri_i/ num_v_smooth_Ri, num_v_smooth_BV
192 1 , num_z_smooth_sh, num_m_smooth_sh
193
194 COMMON /kmixcri_r/ Riinfty, BVSQcon
195 1 , difm0, difs0, dift0
196 2 , difmcon, difscon, diftcon
197
198 c-----------------------------------------------------------------------
199 c parameters for subroutine "ddmix"
200 c
201 c
202 c to compute additional diffusivity due to double diffusion:
203 c
204 c Rrho0 = limit for double diffusive density ratio
205 c dsfmax = maximum diffusivity in case of salt fingering (m2/s)
206 c-----------------------------------------------------------------------
207
208 _KPP_RL Rrho0, dsfmax
209 common /kmixcdd/ Rrho0, dsfmax
210
211 c-----------------------------------------------------------------------
212 c parameters for subroutine "blmix"
213 c
214 c
215 c to compute mixing within boundary layer:
216 c
217 c cstar = proportionality coefficient for nonlocal transport
218 c cg = non-dimensional coefficient for counter-gradient term
219 c-----------------------------------------------------------------------
220
221 _KPP_RL cstar, cg
222
223 common /kmixcbm/ cstar, cg
224
225 #endif /* ALLOW_KPP */
  ViewVC Help
Powered by ViewVC 1.1.22