/[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.14 - (show annotations) (download)
Fri Sep 18 11:40:22 2009 UTC (14 years, 7 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint62v, checkpoint62u, checkpoint62t, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62w, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint63g, checkpoint64, checkpoint65, checkpoint62, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint65o, checkpoint62b, checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint61v, checkpoint61w, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.13: +5 -2 lines
File MIME type: text/plain 
add double diffusive contributions as a hack, turned off by default
for now, and the code can be excluded with a CPP-flag
EXCLUDE_KPP_DOUBLEDIFF just as EXCLUDE_KPP_SHEAR_MIX
1 C $Header: /u/gcmpack/MITgcm/pkg/kpp/KPP_PARAMS.h,v 1.13 2009/05/28 22:59:22 dfer 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
17 integer mdiff, Nrm1, Nrp1, Nrp2
18 integer imt
19 parameter( mdiff = 3 )
20 parameter( Nrm1 = Nr-1 )
21 parameter( Nrp1 = Nr+1 )
22 parameter( Nrp2 = Nr+2 )
23 parameter( imt=(sNx+2*OLx)*(sNy+2*OLy) )
24
25 #ifdef ALLOW_KPP
26
27 C Time invariant parameters initialized by subroutine kmixinit
28 C nzmax (nx,ny) - Maximum number of wet levels in each column
29 C pMask - Mask relating to Pressure/Tracer point grid.
30 C 0. if P point is on land.
31 C 1. if P point is in water.
32 C Note: use now maskC since pMask is identical to maskC
33 C zgrid (0:Nr+1) - vertical levels of tracers (<=0) (m)
34 C hwide (0:Nr+1) - layer thicknesses (>=0) (m)
35 C kpp_freq - Re-computation frequency for KPP parameters (s)
36 C kpp_dumpFreq - KPP dump frequency. (s)
37 C kpp_taveFreq - KPP time-averaging frequency. (s)
38
39 INTEGER nzmax ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, nSx, nSy )
40 c _RL pMask ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr, nSx, nSy )
41 _RL zgrid ( 0:Nr+1 )
42 _RL hwide ( 0:Nr+1 )
43 _RL kpp_freq
44 _RL kpp_dumpFreq
45 _RL kpp_taveFreq
46
47 COMMON /kpp_i/ nzmax
48
49 COMMON /kpp_r1/ zgrid, hwide
50
51 COMMON /kpp_r2/ kpp_freq, kpp_dumpFreq, kpp_taveFreq
52
53
54 C-----------------------------------------------------------------------
55 C
56 C KPP flags and min/max permitted values for mixing parameters
57 c
58 C KPPmixingMaps - if true, include KPP diagnostic maps in STDOUT
59 C KPPwriteState - if true, write KPP state to file
60 C minKPPhbl - KPPhbl minimum value (m)
61 C KPP_ghatUseTotalDiffus -
62 C if T : Compute the non-local term using
63 C the total vertical diffusivity ;
64 C if F (=default): use KPP vertical diffusivity
65 C Note: prior to checkpoint55h_post, was using the total Kz
66 C KPPuseDoubleDiff - if TRUE, include double diffusive
67 C contributions
68 C LimitHblStable - if TRUE (the default), limits the depth of the
69 C hbl under stable conditions.
70 C
71 C-----------------------------------------------------------------------
72
73 LOGICAL KPPmixingMaps, KPPwriteState, KPP_ghatUseTotalDiffus
74 LOGICAL KPPuseDoubleDiff
75 LOGICAL LimitHblStable
76 COMMON /KPP_PARM_L/
77 & KPPmixingMaps, KPPwriteState, KPP_ghatUseTotalDiffus,
78 & KPPuseDoubleDiff, LimitHblStable
79
80 _RL minKPPhbl
81 COMMON /KPP_PARM_R/ minKPPhbl
82
83 c====================== file "kmixcom.h" =======================
84 c
85 c-----------------------------------------------------------------------
86 c Define various parameters and common blocks for KPP vertical-
87 c mixing scheme; used in "kppmix.F" subroutines.
88 c Constants are set in subroutine "ini_parms".
89 c-----------------------------------------------------------------------
90 c
91 c-----------------------------------------------------------------------
92 c parameters for several subroutines
93 c
94 c epsln = 1.0e-20
95 c phepsi = 1.0e-10
96 c epsilon = nondimensional extent of the surface layer = 0.1
97 c vonk = von Karmans constant = 0.4
98 c dB_dz = maximum dB/dz in mixed layer hMix = 5.2e-5 s^-2
99 c conc1,conam,concm,conc2,zetam,conas,concs,conc3,zetas
100 c = scalar coefficients
101 c-----------------------------------------------------------------------
102
103 _RL epsln,phepsi,epsilon,vonk,dB_dz,
104 $ conc1,
105 $ conam,concm,conc2,zetam,
106 $ conas,concs,conc3,zetas
107
108 common /kmixcom/ epsln,phepsi,epsilon,vonk,dB_dz,
109 $ conc1,
110 $ conam,concm,conc2,zetam,
111 $ conas,concs,conc3,zetas
112
113 c-----------------------------------------------------------------------
114 c parameters for subroutine "bldepth"
115 c
116 c
117 c to compute depth of boundary layer:
118 c
119 c Ricr = critical bulk Richardson Number = 0.3
120 c cekman = coefficient for ekman depth = 0.7
121 c cmonob = coefficient for Monin-Obukhov depth = 1.0
122 c concv = ratio of interior buoyancy frequency to
123 c buoyancy frequency at entrainment depth = 1.8
124 c hbf = fraction of bounadry layer depth to
125 c which absorbed solar radiation
126 c contributes to surface buoyancy forcing = 1.0
127 c Vtc = non-dimensional coefficient for velocity
128 c scale of turbulant velocity shear
129 c (=function of concv,concs,epsilon,vonk,Ricr)
130 c-----------------------------------------------------------------------
131
132 _RL Ricr,cekman,cmonob,concv,Vtc
133 _RL hbf
134
135 common /kpp_bldepth1/ Ricr,cekman,cmonob,concv,Vtc
136 common /kpp_bldepth2/ hbf
137
138 c-----------------------------------------------------------------------
139 c parameters and common arrays for subroutines "kmixinit"
140 c and "wscale"
141 c
142 c
143 c to compute turbulent velocity scales:
144 c
145 c nni = number of values for zehat in the look up table
146 c nnj = number of values for ustar in the look up table
147 c
148 c wmt = lookup table for wm, the turbulent velocity scale
149 c for momentum
150 c wst = lookup table for ws, the turbulent velocity scale
151 c for scalars
152 c deltaz = delta zehat in table
153 c deltau = delta ustar in table
154 c zmin = minimum limit for zehat in table (m3/s3)
155 c zmax = maximum limit for zehat in table
156 c umin = minimum limit for ustar in table (m/s)
157 c umax = maximum limit for ustar in table
158 c-----------------------------------------------------------------------
159
160 integer nni , nnj
161 parameter (nni = 890, nnj = 480)
162
163 _RL wmt(0:nni+1,0:nnj+1), wst(0:nni+1,0:nnj+1)
164 _RL deltaz,deltau,zmin,zmax,umin,umax
165 common /kmixcws/ wmt, wst
166 $ , deltaz,deltau,zmin,zmax,umin,umax
167
168 c-----------------------------------------------------------------------
169 c parameters for subroutine "ri_iwmix"
170 c
171 c
172 c to compute vertical mixing coefficients below boundary layer:
173 c
174 c num_v_smooth_Ri = number of times Ri is vertically smoothed
175 c num_v_smooth_BV, num_z_smooth_sh, and num_m_smooth_sh are dummy
176 c variables kept for backward compatibility of the data file
177 c Riinfty = local Richardson Number limit for shear instability = 0.7
178 c BVSQcon = Brunt-Vaisala squared (1/s^2)
179 c difm0 = viscosity max due to shear instability (m^2/s)
180 c difs0 = tracer diffusivity .. (m^2/s)
181 c dift0 = heat diffusivity .. (m^2/s)
182 c difmcon = viscosity due to convective instability (m^2/s)
183 c difscon = tracer diffusivity .. (m^2/s)
184 c diftcon = heat diffusivity .. (m^2/s)
185 c-----------------------------------------------------------------------
186
187 INTEGER num_v_smooth_Ri, num_v_smooth_BV
188 INTEGER num_z_smooth_sh, num_m_smooth_sh
189 COMMON /kmixcri_i/ num_v_smooth_Ri, num_v_smooth_BV
190 1 , num_z_smooth_sh, num_m_smooth_sh
191
192 _RL Riinfty, BVSQcon
193 _RL difm0 , difs0 , dift0
194 _RL difmcon, difscon, diftcon
195 COMMON /kmixcri_r/ Riinfty, BVSQcon
196 1 , difm0, difs0, dift0
197 2 , difmcon, difscon, diftcon
198
199 c-----------------------------------------------------------------------
200 c parameters for subroutine "ddmix"
201 c
202 c
203 c to compute additional diffusivity due to double diffusion:
204 c
205 c Rrho0 = limit for double diffusive density ratio
206 c dsfmax = maximum diffusivity in case of salt fingering (m2/s)
207 c-----------------------------------------------------------------------
208
209 _RL Rrho0, dsfmax
210 common /kmixcdd/ Rrho0, dsfmax
211
212 c-----------------------------------------------------------------------
213 c parameters for subroutine "blmix"
214 c
215 c
216 c to compute mixing within boundary layer:
217 c
218 c cstar = proportionality coefficient for nonlocal transport
219 c cg = non-dimensional coefficient for counter-gradient term
220 c-----------------------------------------------------------------------
221
222 _RL cstar, cg
223 common /kmixcbm/ cstar, cg
224
225 #endif /* ALLOW_KPP */
226
227
228 CEH3 ;;; Local Variables: ***
229 CEH3 ;;; mode:fortran ***
230 CEH3 ;;; End: ***
  ViewVC Help
Powered by ViewVC 1.1.22