/[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.7 - (hide 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 cnh 1.7 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 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     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 adcroft 1.1 #ifdef FRUGAL_KPP
25 heimbach 1.5 parameter( imt=(sNx+2)*(sNy+2) )
26     parameter( ibot=0, itop=sNx+1, jbot=0, jtop=sNy+1 )
27 adcroft 1.1 #else
28 heimbach 1.5 parameter( imt=(sNx+2*OLx)*(sNy+2*OLy) )
29     parameter( ibot=1-OLx, itop=sNx+OLx, jbot=1-OLy, jtop=sNy+OLy )
30 adcroft 1.1 #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 adcroft 1.2 C kpp_dumpFreq - KPP dump frequency. (s)
43     C kpp_taveFreq - KPP time-averaging frequency. (s)
44 heimbach 1.3
45 adcroft 1.1
46 heimbach 1.5 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 adcroft 1.1 _RL kpp_freq
51 adcroft 1.2 _RL kpp_dumpFreq
52     _RL kpp_taveFreq
53 adcroft 1.1
54 heimbach 1.5 COMMON /kpp_i/ nzmax
55 adcroft 1.1
56 heimbach 1.5 COMMON /kpp_r1/ pMask, zgrid, hwide
57 adcroft 1.1
58 heimbach 1.5 COMMON /kpp_r2/ kpp_freq, kpp_dumpFreq, kpp_taveFreq
59 adcroft 1.1
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 heimbach 1.5 _KPP_RL minKPPhbl
77 adcroft 1.1
78     COMMON /KPP_PARM_R/
79 heimbach 1.5 & minKPPhbl
80 adcroft 1.1
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 heimbach 1.5 _KPP_RL epsln,phepsi,epsilon,vonk,dB_dz,
102 adcroft 1.1 $ 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 heimbach 1.5 _KPP_RL Ricr,cekman,cmonob,concv,Vtc
131     _RL hbf
132 adcroft 1.1
133 heimbach 1.5 common /kpp_bldepth1/ Ricr,cekman,cmonob,concv,Vtc
134     common /kpp_bldepth2/ hbf
135 adcroft 1.1
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 heimbach 1.5 _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 adcroft 1.1 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 heimbach 1.5 _KPP_RL Riinfty, BVSQcon
188     _KPP_RL difm0 , difs0 , dift0
189     _KPP_RL difmcon, difscon, diftcon
190 adcroft 1.1
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 heimbach 1.5 _KPP_RL Rrho0, dsfmax
209 adcroft 1.1 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 heimbach 1.5 _KPP_RL cstar, cg
222 adcroft 1.1
223     common /kmixcbm/ cstar, cg
224    
225     #endif /* ALLOW_KPP */
  ViewVC Help
Powered by ViewVC 1.1.22