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

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

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

Revision 1.14 - (hide annotations) (download)
Fri Sep 18 11:40:22 2009 UTC (14 years, 8 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 mlosch 1.14 C $Header: /u/gcmpack/MITgcm/pkg/kpp/KPP_PARAMS.h,v 1.13 2009/05/28 22:59:22 dfer Exp $
2 jmc 1.8 C $Name: $
3 adcroft 1.1
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 heimbach 1.5 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 dimitri 1.11 integer imt
19 heimbach 1.5 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 adcroft 1.1
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 jmc 1.8 C Note: use now maskC since pMask is identical to maskC
33 adcroft 1.1 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 adcroft 1.2 C kpp_dumpFreq - KPP dump frequency. (s)
37     C kpp_taveFreq - KPP time-averaging frequency. (s)
38 heimbach 1.3
39 heimbach 1.5 INTEGER nzmax ( 1-OLx:sNx+OLx, 1-OLy:sNy+OLy, nSx, nSy )
40 dimitri 1.12 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 adcroft 1.1 _RL kpp_freq
44 adcroft 1.2 _RL kpp_dumpFreq
45     _RL kpp_taveFreq
46 adcroft 1.1
47 heimbach 1.5 COMMON /kpp_i/ nzmax
48 adcroft 1.1
49 jmc 1.8 COMMON /kpp_r1/ zgrid, hwide
50 adcroft 1.1
51 heimbach 1.5 COMMON /kpp_r2/ kpp_freq, kpp_dumpFreq, kpp_taveFreq
52 adcroft 1.1
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 jmc 1.10 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 mlosch 1.14 C KPPuseDoubleDiff - if TRUE, include double diffusive
67     C contributions
68 dfer 1.13 C LimitHblStable - if TRUE (the default), limits the depth of the
69     C hbl under stable conditions.
70 adcroft 1.1 C
71     C-----------------------------------------------------------------------
72    
73 jmc 1.10 LOGICAL KPPmixingMaps, KPPwriteState, KPP_ghatUseTotalDiffus
74 mlosch 1.14 LOGICAL KPPuseDoubleDiff
75 dfer 1.13 LOGICAL LimitHblStable
76 adcroft 1.1 COMMON /KPP_PARM_L/
77 dfer 1.13 & KPPmixingMaps, KPPwriteState, KPP_ghatUseTotalDiffus,
78 mlosch 1.14 & KPPuseDoubleDiff, LimitHblStable
79 adcroft 1.1
80 dimitri 1.12 _RL minKPPhbl
81     COMMON /KPP_PARM_R/ minKPPhbl
82 adcroft 1.1
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 dimitri 1.12 _RL epsln,phepsi,epsilon,vonk,dB_dz,
104 adcroft 1.1 $ 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 dimitri 1.12 _RL Ricr,cekman,cmonob,concv,Vtc
133 heimbach 1.5 _RL hbf
134 adcroft 1.1
135 heimbach 1.5 common /kpp_bldepth1/ Ricr,cekman,cmonob,concv,Vtc
136     common /kpp_bldepth2/ hbf
137 adcroft 1.1
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 dimitri 1.12 _RL wmt(0:nni+1,0:nnj+1), wst(0:nni+1,0:nnj+1)
164     _RL deltaz,deltau,zmin,zmax,umin,umax
165 adcroft 1.1 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 dimitri 1.12 INTEGER num_v_smooth_Ri, num_v_smooth_BV
188     INTEGER num_z_smooth_sh, num_m_smooth_sh
189 adcroft 1.1 COMMON /kmixcri_i/ num_v_smooth_Ri, num_v_smooth_BV
190     1 , num_z_smooth_sh, num_m_smooth_sh
191    
192 dimitri 1.12 _RL Riinfty, BVSQcon
193     _RL difm0 , difs0 , dift0
194     _RL difmcon, difscon, diftcon
195 adcroft 1.1 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 dimitri 1.12 _RL Rrho0, dsfmax
210 adcroft 1.1 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 dimitri 1.12 _RL cstar, cg
223 adcroft 1.1 common /kmixcbm/ cstar, cg
224    
225     #endif /* ALLOW_KPP */
226 edhill 1.9
227    
228     CEH3 ;;; Local Variables: ***
229     CEH3 ;;; mode:fortran ***
230     CEH3 ;;; End: ***
  ViewVC Help
Powered by ViewVC 1.1.22