Context Navigation

Registry.EM.newphys @ 1581

Last change on this file since 1581 was 1236, checked in by aslmd, 12 years ago

MESOSCALE. A necessary complement commit to r1234 where a upgraded interface making use of modules was proposed. Completed the new formulation for module_lmd_driver for newphys with improved interface with both ini/bdy conditions and physical parameterizations. Changed the Registry accordingly. Finished changes about I/O with the LMD physics (see LMDZ.MARS/README). Made all those changes compatible for old interface, and LES runs (checked with test cases), as well as old input files. Changed makemeso to account for full flexibility on changin nx ny ntracers nproc with newphys. Cleaned the now obsolete bits of code used in LMD physics shared with the GCM. ----- Everything is now ready to properly code both restart runs and nesting for mesoscale runs with new physics.

File size: 147.9 KB

Line
1	## WARNING: this file is autogenerated from Registry/Registry.EM. Changes may be lost
2	# Registry file, EM
3	#
4	# At the present time this file is managed manually and edited by hand.
5	#
6	################################################################################
7	# Dimension specifications
8	#
9	# This section of the Registry file is used to specify the dimensions
10	# that will be used to define arrays. Dim is the one-letter name of the
11	# dimension. How defined can either be "standard_domain", which means
12	# that the dimension (1) is one of the three spatial dimensions and (2)
13	# it will be set using the standard namelist mechanism and domain data
14	# structure dimension fields (e.g. sd31,ed31,sd32...).
15	#
16	# Order refers to which of the three sets of just-mentioned internal
17	# dimension variables the dimension is referred to by in the driver.
18	# That is, is it the first, second, or third dimension. The registry
19	# infers the mapping of its internal dimensions according to the
20	# combination of Order and Coord-axis that are specified in this table.
21	# Note that it is all right to more than one dimension name for, say, the
22	# x dimension. However, the Order and Coord-axis relationship must be
23	# consistent throughout.
24	#
25	# Note: these entries do not enforce storage order on a particular field.
26	# That is determined by the dimension strings for each field. But it does
27	# relate the dimspec to the internal data structures that the driver uses
28	# to maintain the three physical domain dimensions.
29	#
30	# "How defined" can also specify the name of a namelist variable from which
31	# the definition for the dimension will come; this is specified as
32	# "namelist=<variable name>". The namelist variable must have been
33	# defined as an integer and with only one entry in the rconfig table. Or
34	# a constant can be specified. The coordinate axis for the dimension is
35	# either X, Y, Z, or C (for "not a spatial dimension"). The Dimname is
36	# the descriptive name of the dimension that will be included in the
37	# metadata in data sets. Note that the b, f, and t modifiers that appear
38	# as the last characters of dimension strings used # in state and # i1
39	# registry definitions are not dimensions and do not need to be declared
40	# here.
41	#
42
43	# Available characters for dimspec: 0123456789@%+=\|?.!&[{}]
44
45	#<Table> <Dim> <Order> <How defined> <Coord-axis> <Dimname in Datasets>
46	dimspec i 1 standard_domain x west_east
47	dimspec j 3 standard_domain y south_north
48	dimspec k 2 standard_domain z bottom_top
49	dimspec l 2 namelist=num_soil_layers z soil_layers
50	dimspec u 2 namelist=num_land_cat z land_cat
51	dimspec s 2 namelist=num_soil_cat z soil_cat
52	dimspec p - constant=7501 c microphysics_rstrt_state
53	dimspec w - namelist=spec_bdy_width c spec_bdy_width
54	dimspec e 3 namelist=ensdim z ensemble dimension
55	dimspec z - namelist=max_obs c max_obs
56	dimspec h - namelist=nobs_err_flds c nobs_err_flds
57	dimspec r - namelist=nobs_ndg_vars c nobs_ndg_vars
58	dimspec g 2 namelist=num_metgrid_levels z num_metgrid_levels
59	dimspec m 2 constant=12 z months_per_year
60	dimspec a - namelist=cam_abs_dim1 c cam_abs_dim1
61	dimspec c - namelist=cam_abs_dim2 z cam_abs_dim2
62	dimspec q 2 namelist=levsiz z levsiz
63	dimspec d 2 namelist=paerlev z paerlev
64	dimspec v - constant=1 z one
65
66
67	################################################################################
68	################################################################################
69	################################################################################
70
71
72	####MARS MARS MARS MARS MARS
73	####MARS MARS MARS MARS MARS
74	####
75	#### variables physiques ecrites dans les fichiers WRF
76	#### - h indique l'ecriture
77	####
78	####
79	#### RADIATIVE TRANSFER
80	state real HR_SW ikj misc 1 - h "HR_SW" "HEATING RATE SW" "K/s"
81	state real HR_LW ikj misc 1 - h "HR_LW" "HEATING RATE LW" "K/s"
82	state real SWDOWNZ ij misc 1 - h "SWDOWNZ" "DOWNWARD SW FLUX AT SURFACE" "W m-2"
83	#### DUST CYCLE
84	state real TAU_DUST ij misc 1 - h "TAU_DUST" "REFERENCE VISIBLE DUST OPACITY" ""
85	state real RDUST ikj misc 1 - h "RDUST" "DUST RADIUS" "m"
86	state real QSURFDUST ij misc 1 - h "QSURFDUST" "DUST MASS AT SURFACE" "kg m-2"
87	#### CLOUDS
88	state real MTOT ij misc 1 - h "MTOT" "TOTAL MASS WATER VAPOR in pmic" "pmic"
89	state real ICETOT ij misc 1 - h "ICETOT" "TOTAL MASS WATER ICE" "kg m-2"
90	state real VMR_ICE ikj misc 1 - h "VMR_ICE" "VOL. MIXING RATIO ICE" "ppm"
91	state real TAU_ICE ij misc 1 - h "TAU_ICE" "CLOUD OD at 825 cm-1 TES" ""
92	state real RICE ikj misc 1 - h "RICE" "ICE RADIUS" "m"
93	#### THERMALS
94	state real ZMAX ij misc 1 - h "ZMAX" "MAXIMUM LEVEL THERMALS" "m"
95	state real HFMAX ij misc 1 - h "HFMAX" "MAXIMUM PBL HEAT FLUX" "m.K/s"
96	####
97	####
98	####MARS MARS MARS MARS MARS
99	####MARS MARS MARS MARS MARS
100
101	####MARS MARS MARS MARS MARS
102	####MARS MARS MARS MARS MARS
103	####
104	####
105	#### TRACEURS: l'ordre dans SCALAR est defini plus bas dans les options mars)
106	#### - effacer le h si l'on veut eviter les sorties
107	####
108	state real - ikjftb scalar 1 - - -
109	state real qh2o ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2O" "Water vapor mixing ratio" "kg kg-1"
110	state real qh2o_ice ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QH2O_ICE" "Water ice mixing ratio" "kg kg-1"
111	state real qdust ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QDUST" "Dust mixing ratio" "kg kg-1"
112	state real qdustn ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QDUSTN" "Dust_number mixing ratio" "kg kg-1"
113	state real qccn ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCCN" "CCN mixing ratio" "kg kg-1"
114	state real qccnn ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCCNN" "CCN number mixing ratio" "kg kg-1"
115	state real qco2 ikjftb scalar 1 - i01rhusdf=(bdy_interp:dt) "QCO2" "CO2 mixing ratio" "kg kg-1"
116	state real chem_co ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_co" "" ""
117	state real chem_o ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_o" "" ""
118	state real chem_o1d ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_o1d" "" ""
119	state real chem_o2 ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_o2" "" ""
120	state real chem_o3 ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_o3" "" ""
121	state real chem_h ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_h" "" ""
122	state real chem_h2 ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_h2" "" ""
123	state real chem_oh ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_oh" "" ""
124	state real chem_ho2 ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_ho2" "" ""
125	state real chem_h2o2 ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_h2o2" "" ""
126	state real chem_ch4 ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_ch4" "" ""
127	state real chem_n2 ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_n2" "" ""
128	state real chem_ar ikjftb scalar 1 - i01rusdf=(bdy_interp:dt) "chem_ar" "" ""
129	####
130	####
131	####MARS MARS MARS MARS MARS
132	####MARS MARS MARS MARS MARS
133
134
135
136	################################################################################
137	################################################################################
138	################################################################################
139
140
141
142	# Lines that start with the word 'state' form a table that is
143	# used by the script use_registry to generate module_state_descript.F
144	# and other files. Also see documentation in use_registry.
145	#
146	# table entries are of the form
147	#<Table> <Type> <Sym> <Dims> <Use> <NumTLev> <Stagger> <IO> <DNAME> <DESCRIP> <UNITS>
148	#
149
150	# It is required that LU_INDEX appears before any variable that is
151	# interpolated with a mask, as lu_index supplies that mask.
152	# this next 1 is for the HFSoLE/PET demo; writing these to auxhist1 output over MCEL for coupling
153	# with wave model, only if compiled with -DMCELIO, JM 2003/05/29
154	state real LU_INDEX ij misc 1 - ir01d=(interp_fcnm)u=(copy_fcnm) "LU_INDEX" "LAND USE CATEGORY" ""
155	state real LU_MASK ij misc 1 - i3h1 "LU_MASK" "0 land 1 water" ""
156
157	#
158	# Variables from WPS
159	#
160	state real u_gc igj dyn_em 1 XZ i1 "UU" "x-wind component" "m s-1"
161	state real v_gc igj dyn_em 1 YZ i1 "VV" "y-wind component" "m s-1"
162	state real t_gc igj dyn_em 1 Z i1 "TT" "temperature" "K"
163	state real rh_gc igj dyn_em 1 Z i1 "RH" "relative humidity" "%"
164	state real ght_gc igj dyn_em 1 Z i1 "GHT" "geopotential height" "m"
165	state real p_gc igj dyn_em 1 Z i1 "PRES" "pressure" "Pa"
166	state real xlat_gc ij dyn_em 1 - i1 "XLAT_M" "latitude, positive north" "degrees"
167	state real xlong_gc ij dyn_em 1 - i1 "XLONG_M" "longitude, positive east" "degrees"
168	state real ht_gc ij dyn_em 1 - i1 "HGT_M" "topography elevation" "m"
169	state real tsk_gc ij dyn_em 1 - i1 "SKINTEMP" "skin temperature" "K"
170	state real tavgsfc ij dyn_em 1 - i1 "TAVGSFC" "daily mean of surface air temperature" "K"
171	state real tmn_gc ij dyn_em 1 - i1 "SOILTEMP" "annual mean deep soil temperature" "K"
172	state real pslv_gc ij dyn_em 1 - i1 "PMSL" "sea level pressure" "Pa"
173	state real greenfrac imj dyn_em 1 Z i1 "GREENFRAC" "monthly greenness fraction" "0 - 1 fraction"
174	state real albedo12m imj dyn_em 1 Z i1 "ALBEDO12M" "background albedo" "0 - 1 fraction"
175	state real pd_gc igj dyn_em 1 Z - "PD" "dry pressure" "Pa"
176	state real psfc_gc ij dyn_em 1 - - "PSFC_GC" "surface pressure" "Pa"
177	state real intq_gc ij dyn_em 1 - - "INTQ" "integrated mixing ratio" "Pa"
178	state real pdhs ij dyn_em 1 - - "PDHS" "hydrostatic dry surface pressure" "Pa"
179	state real qv_gc igj dyn_em 1 Z i1 "QV" "mixing ratio" "kg kg-1"
180	#state real qr_gc igj dyn_em 1 Z i1 "QR" "rain water mixing ratio" "kg kg-1"
181	#state real qc_gc igj dyn_em 1 Z i1 "QC" "cloud water mixing ratio" "kg kg-1"
182	#state real qs_gc igj dyn_em 1 Z i1 "QS" "snow mixing ratio" "kg kg-1"
183	#state real qi_gc igj dyn_em 1 Z i1 "QI" "cloud ice mixing ratio" "kg kg-1"
184	#state real qg_gc igj dyn_em 1 Z i1 "QG" "graupel mixing ratio" "kg kg-1"
185
186	####MARS MARS MARS MARS MARS
187	####MARS MARS MARS MARS MARS
188	####
189	#### INITIALISATION
190	####
191	#### variables dans real
192	#### -- si on change le nom on doit le changer jusque readmeteo !
193	####
194	############# in module_initialize ################# in geo_em.d**.nc ############################################
195	state real albedo_gcm_gc ij dyn_em 1 - i1 "ALBEDO_GCM" "albedo of naked ground" "0 - 1 fraction"
196	state real therm_inert_gc ij dyn_em 1 - i1 "THERMAL_INERTIA" "thermal inertia" "J/m2/K/s0.5"
197	state real z0_gc ij dyn_em 1 - i1 "SURF_ROUGHNESS" "surface roughness" "cm"
198	state real slpx_gc ij dyn_em 1 - i1 "SLPX" "slope x direction" "rad"
199	state real slpy_gc ij dyn_em 1 - i1 "SLPY" "slope y direction" "rad"
200	state real hv_gc igj dyn_em 1 Z i1 "HV" "water vapor" "kg/kg"
201	state real hi_gc igj dyn_em 1 Z i1 "HI" "water ice" "kg/kg"
202	state real tsoil_gc igj dyn_em 1 Z i1 "TSOIL" "soil temperatures" "K"
203	state real isoil_gc igj dyn_em 1 Z i1 "ISOIL" "soil thermal inertia" "tiu"
204	state real dsoil_gc igj dyn_em 1 Z i1 "DSOIL" "soil depths" "m"
205	state real co2_gc igj dyn_em 1 Z i1 "CO2" "co2 mixing ratio" "kg/kg"
206	state real dustq_gc igj dyn_em 1 Z i1 "DUSTQ" "dust mixing ratio" "kg/kg"
207	state real dustn_gc igj dyn_em 1 Z i1 "DUSTN" "dust number density" "part/kg"
208	state real ccnq_gc igj dyn_em 1 Z i1 "CCNQ" "ccn mixing ratio" "kg/kg"
209	state real ccnn_gc igj dyn_em 1 Z i1 "CCNN" "ccn number density" "part/kg"
210	####
211	#### variables dans wrf
212	####
213	#### 'm_' prefix is to avoid redundancies with possible existing terrestrial variables
214	####
215	#### vvvv nom pour module_initialize, solve_em, lmd_driver (interne fortran)
216	#### vvvv nom fichier netCDF (important initialisation
217	################################################################################################################
218	state real m_albedo ij misc 1 - i012rd "MARS_ALB" "albedo of naked ground" "0 - 1 fraction"
219	state real m_ti ij misc 1 - i012rd "MARS_TI" "thermal inertia" "J/m2/K/s0.5"
220	state real m_z0 ij misc 1 - i012rd "MARS_Z0" "surface roughness" "m"
221	state real slpx ij misc 1 - i012rd "SLOPEX" "slope x direction" "rad"
222	state real slpy ij misc 1 - i012rd "SLOPEY" "slope y direction" "rad"
223	state real m_emiss ij misc 1 - i012rd "MARS_EMISS" "emissivity" "0 - 1 fraction"
224	state real m_co2ice ij misc 1 - i012rhd "MARS_CICE" "CO2 GROUND ICE" "kg/m2"
225	state real m_h2oice ij misc 1 - i012rhd "MARS_WICE" "H2O GROUND ICE" "kg/m2"
226	state real m_fluxrad ij misc 1 - i012rhd "MARS_FLUXRAD" "net rad flux at surf" "W.m-2"
227	state real m_wstar ij misc 1 - i012rhd "MARS_WSTAR" "FREE CONVECTION VELOCITY" "m.s-1"
228	state real m_tsurf ij misc 1 - i012rhd "MARS_TSURF" "SURFACE TEMPERATURE" "K"
229	state real m_tsoil ilj misc 1 - i012rhd "MARS_TSOIL" "SOIL TEMPERATURES" "K"
230	state real m_q2 ikj misc 1 Z i012rhd "MARS_Q2" "parameterized TKE" "kg.m-3"
231	state real m_isoil ilj misc 1 - i012rd "MARS_ISOIL" "soil thermal inertia" "tiu"
232	state real m_dsoil ilj misc 1 - i012rd "MARS_DSOIL" "soil depths" "m"
233	state real m_gw ilj misc 1 - i012rd "MARS_GW" "gwparam" " "
234	####
235	####
236	####MARS MARS MARS MARS MARS
237	####MARS MARS MARS MARS MARS
238
239	#-----------------------------------------------------------------------------------------------------------------------------------------------------------------
240	#
241	# Variables for Eulerian mass coordinate dynamics
242	#
243
244	# Velocities
245	#
246	# U Vel
247	state real u ikjb dyn_em 2 X \
248	i01rhusdf=(bdy_interp:dt) "U" "x-wind component" "m s-1"
249	state real ru ikj dyn_em 1 X - "MU_U" "mu-coupled u" "Pa m s-1"
250	state real ru_m ikj dyn_em 1 X - "ru_m" "" ""
251	state real ru_tend ikj dyn_em 1 X - "ru_tend" "" ""
252	i1 real ru_tendf ikj dyn_em 1 X
253	state real u_save ikj dyn_em 1 X - "u_save"
254	#
255	# V Vel
256	state real v ikjb dyn_em 2 Y \
257	i01rhusdf=(bdy_interp:dt) "V" "y-wind component" "m s-1"
258	state real rv ikj dyn_em 1 Y - "MU_V" "mu-coupled v" "Pa m s-1"
259	state real rv_m ikj dyn_em 1 Y - "rv_m"
260	state real rv_tend ikj dyn_em 1 Y - "rv_tend"
261	i1 real rv_tendf ikj dyn_em 1 Y
262	state real v_save ikj dyn_em 1 Y - "v_save"
263	#
264	# Vertical Vel
265	state real w ikjb dyn_em 2 Z \
266	irhusdf=(bdy_interp:dt) "w" "z-wind component" "m s-1"
267	state real ww ikj dyn_em 1 Z r "ww" "mu-coupled eta-dot" "Pa s-1"
268	state real rw ikj dyn_em 1 Z - "rw" "mu-coupled w" "Pa m s-1"
269	i1 real ww1 ikj dyn_em 1 Z
270	state real ww_m ikj dyn_em 1 Z r "ww_m" "time-avg mu-coupled eta-dot" "Pa s-1"
271	i1 real wwp ikj dyn_em 1 Z
272	i1 real rw_tend ikj dyn_em 1 Z
273	i1 real rw_tendf ikj dyn_em 1 Z
274	i1 real w_save ikj dyn_em 1 Z
275
276	# Geopotential
277	state real ph ikjb dyn_em 2 Z \
278	irusdf=(bdy_interp:dt) "ph" "perturbation geopotential" "m2 s-2"
279	state real phb ikj dyn_em 1 Z irdus "phb" "base-state geopotential" "m2 s-2"
280	state real phb_fine ikj dyn_em 1 Z - "phb_fine" "for nesting, temp holding interpolated coarse grid phb" "m2 s-2"
281	state real ph0 ikj dyn_em 1 Z r "ph0" "initial geopotential" "m2 s-2"
282	state real php ikj dyn_em 1 Z r "php" "geopotential" "m2 s-2"
283	i1 real ph_tend ikj dyn_em 1 Z
284	i1 real ph_tendf ikj dyn_em 1 Z
285	i1 real ph_save ikj dyn_em 1 Z
286	##
287	## MARS: TO GET LIGHTER HISTORY FILES
288	##
289	state real phtot ikj dyn_em 1 Z irh "PHTOT" "total geopotential" "m2 s-2"
290
291	# Potential Temperature
292	state real t ikjb dyn_em 2 - \
293	i01rhusdf=(bdy_interp:dt) "t" "perturbation potential temperature (theta-t0)" "K"
294
295	state real t_init ikj dyn_em 1 - ir "t_init" "initial potential temperature" "K"
296	i1 real t_tend ikj dyn_em 1 -
297	i1 real t_tendf ikj dyn_em 1 -
298	state real tp ikj dyn_em 2 -
299	i1 real t_2save ikj dyn_em 1 -
300	state real t_save ikj dyn_em 1 - "t_save"
301	#
302
303
304	# Mass
305	state real mu ijb dyn_em 2 - \
306	irusdf=(bdy_interp:dt) "mu" "perturbation dry air mass in column" "Pa"
307	state real mub ij dyn_em 1 - irdus "mub" "base state dry air mass in column" "Pa"
308	state real mub_fine ij dyn_em 1 - - "mub_fine" "nest temp, holds interpolated coarse grid mub" "Pa"
309	state real mu0 ij dyn_em 1 - i01rdu "mu0" "initial dry mass in column" "Pa"
310	state real mudf ij dyn_em 1 - - "mudf" "" ""
311	state real muu ij dyn_em 1 - "muu"
312	i1 real muus ij dyn_em 1 -
313	state real muv ij dyn_em 1 - "muv"
314	i1 real muvs ij dyn_em 1 -
315	state real mut ij dyn_em 1 - "mut"
316	state real muts ij dyn_em 1 - "muts"
317	i1 real muave ij dyn_em 1 -
318	i1 real mu_save ij dyn_em 1 -
319	i1 real mu_tend ij dyn_em 1 -
320	i1 real mu_tendf ij dyn_em 1 -
321
322	#diagnostic for looking at nest position in output. A mungy version of terrain height.
323	state real nest_pos ij misc 1 - rhu=(mark_domain) "NEST_POS"
324	state real nest_mask ij misc 1 - ru=(mark_domain) "NEST_MASK" "LOCATION OF NEST IF ANY"
325	state real ht_coarse ij misc 1 - r - "STORAGE FOR LOW-RES TERRAIN"
326
327
328	# TKE
329	state real tke ikj dyn_em 2 - r "tke" "TURBULENCE KINETIC ENERGY" "m2 s-2"
330	i1 real tke_tend ikj dyn_em 1 -
331
332	# Pressure and Density
333	state real p ikj dyn_em 1 - r "p" "perturbation pressure" "Pa"
334	state real al ikj dyn_em 1 - r "al" "inverse perturbation density" "m3 kg-1"
335	state real alt ikj dyn_em 1 - r "alt" "inverse density" "m3 kg-1"
336	state real alb ikj dyn_em 1 - rdus "alb" "inverse base density" "m3 kg-1"
337	state real zx ikj dyn_em 1 X - " " " " " "
338	state real zy ikj dyn_em 1 Y - " " " " " "
339	state real rdz ikj dyn_em 1 Z - " " " " " "
340	state real rdzw ikj dyn_em 1 Z - " " " " " "
341	state real pb ikj dyn_em 1 - rdus "pb" "BASE STATE PRESSURE " "Pa"
342	##
343	## MARS: TO GET LIGHTER HISTORY FILES
344	##
345	state real ptot ikj dyn_em 1 - irh "PTOT" "total pressure" "Pa"
346
347	#
348	# Other dyn
349	#
350	i1 real advect_tend ikj dyn_em 1 -
351	i1 real alpha ikj dyn_em 1 -
352	i1 real a ikj dyn_em 1 -
353	i1 real gamma ikj dyn_em 1 -
354	i1 real c2a ikj dyn_em 1 - -
355	i1 real rho ikj dyn_em 1 - -
356	i1 real phm ikj dyn_em 1 - -
357	i1 real cqu ikj dyn_em 1 - -
358	i1 real cqv ikj dyn_em 1 - -
359	i1 real cqw ikj dyn_em 1 - -
360	i1 real pm1 ikj dyn_em 1 - -
361	#
362	# sortie desactivée
363	#
364	state real sr ij dyn_em 1 - ir "sr" "fraction of frozen precipitation"
365	state real potevp ij dyn_em 1 - "potevp" "whatever"
366	state real snopcx ij dyn_em 1 - "snopcx" "whatever"
367	state real soiltb ij dyn_em 1 - ir "soiltb" "whatever"
368	state real fnm k dyn_em 1 - irh "fnm" "upper weight for vertical stretching" ""
369	state real fnp k dyn_em 1 - irh "fnp" "lower weight for vertical stretching" ""
370	state real rdnw k dyn_em 1 - irh "rdnw" "inverse d(eta) values between full (w) levels" ""
371	state real rdn k dyn_em 1 - irh "rdn" "inverse d(eta) values between half (mass) levels" ""
372	state real dnw k dyn_em 1 - irh "dnw" "d(eta) values between full (w) levels" ""
373	state real dn k dyn_em 1 - irh "dn " "d(eta) values between half (mass) levels" ""
374	state real znu k dyn_em 1 - irh "znu" "eta values on half (mass) levels" ""
375	state real znw k dyn_em 1 Z i01rh "znw" "eta values on full (w) levels" ""
376	state real t_base k dyn_em 1 - ir "t_base" "BASE STATET T IN IDEALIZED CASES" "K"
377	state real z ikj dyn_em 1 - - " " " " " "
378	i1 real mu_3d ikj dyn_em 1 -
379	i1 real z_at_w ikj dyn_em 1 Z
380	state real cfn - misc - - irh "cfn" "extrapolation constant" ""
381	state real cfn1 - misc - - irh "cfn1" "extrapolation constant" ""
382	state integer step_number - misc - - ir "step_number" ""
383
384	#
385	# UN BON EXEMPLE: irhd et i01rh
386	#
387	# 2m and 10m output diagnostics
388	state real Q2 ij misc 1 - ird "Q2" "QV at 2 M" "kg kg-1"
389	state real T2 ij misc 1 - ird "T2" "TEMP at 2 M" "K"
390	state real TH2 ij misc 1 - ird "TH2" "POT TEMP at 2 M" "K"
391	state real PSFC ij misc 1 - i01rh "PSFC" "SFC PRESSURE" "Pa"
392
393	#
394	# ICI irh01d
395	#
396	# these next 2 are for the HFSoLE/PET demo; writing these to auxhist1 output over MCEL for coupling
397	# with wave model, only if compiled with -DMCELIO, JM 2003/05/29
398	state real U10 ij misc 1 - ir01d "U10" "U at 10 M" "m s-1"
399	state real V10 ij misc 1 - ir01d "V10" "V at 10 M" "m s-1"
400
401	# these next 4 are for observational nudging
402	state real uratx ij misc 1 - ir "URATX" "Ratio of U over U10 on mass points " "dimensionless"
403	state real vratx ij misc 1 - ir "VRATX" "Ratio of V over V10 on mass points " "dimensionless"
404	state real tratx ij misc 1 - ir "TRATX" "Ratio of T over TH2 on mass points " "dimensionless"
405	state real obs_savwt hikj dyn_em 1 X - "OBS_SAVWT"
406
407	# Other
408	state real rdx - misc - - irh "rdx" "INVERSE X GRID LENGTH" ""
409	state real rdy - misc - - irh "rdy" "INVERSE Y GRID LENGTH" ""
410	state real dts - misc - - ir "dts" "SMALL TIMESTEP" ""
411	state real dtseps - misc - - ir "dtseps" "TIME WEIGHT CONSTANT FOR SMALL STEPS" ""
412	state real resm - misc - - irh "resm" "TIME WEIGHT CONSTANT FOR SMALL STEPS" ""
413	state real zetatop - misc - - irh "zetatop" "ZETA AT MODEL TOP" ""
414	state real cf1 - misc - - irh "cf1" "2nd order extrapolation constant" ""
415	state real cf2 - misc - - irh "cf2" "2nd order extrapolation constant" ""
416	state real cf3 - misc - - irh "cf3" "2nd order extrapolation constant" ""
417	state integer number_at_same_level - - - - - "number_at_same_level" "" ""
418
419	# State for derived time quantities.
420	state integer itimestep - - - - rh "itimestep" "" ""
421	state real xtime - - - - rh "xtime" "minutes since simulation start" ""
422	state real julian - - - - - "julian" "day of year, 0.0 at 0Z on 1 Jan." "days"
423
424
425	# input file descriptor for lbcs on parent domain
426	state integer lbc_fid - - - - - "lbc_fid" "" ""
427	# indicates if tiling has been computed
428	state logical tiled - - - - - "tiled" "" ""
429	# indicates if patches have been computed
430	state logical patched - - - - - "patched" "" ""
431	# indicates whether to read input from file or generate
432	#state logical input_from_file - - - - - "input_from_file" "" ""
433
434	# Mask for moving nest interpolations
435	state integer imask_nostag ij misc 1 -
436	state integer imask_xstag ij misc 1 X
437	state integer imask_ystag ij misc 1 Y
438	state integer imask_xystag ij misc 1 XY
439	# vortex center indices; need for restarts of moving nests
440	state real xi - misc - - r
441	state real xj - misc - - r
442	state real vc_i - misc - - r
443	state real vc_j - misc - - r
444
445	#-----------------------------------------------------------------------------------------------------------------------------------------------------------------
446	# Scalar (4D) arrays
447
448	# Moist Scalars
449	#
450	# The first line ensures that there will be identifiers named moist and
451	# moist_tend even if there are not any moist scalars (so the essentially
452	# dry code will will still link properly)
453	#
454	#-----------------------------
455	# afin d eviter les mic macs terrestres, ceux là ne sont pas a utiliser
456	# reglés a zero
457	#-----------------------------
458	state real - ikjftb moist 1 - - -
459	state real qv ikjftb moist 1 - \
460	i01rusdf=(bdy_interp:dt) "QVAPOR" "Water vapor mixing ratio" "kg kg-1"
461	state real qc ikjftb moist 1 - \
462	i01rusdf=(bdy_interp:dt) "QCLOUD" "Cloud water mixing ratio" "kg kg-1"
463	state real qr ikjftb moist 1 - \
464	i01rusdf=(bdy_interp:dt) "QRAIN" "Rain water mixing ratio" "kg kg-1"
465	state real qi ikjftb moist 1 - \
466	i01rusdf=(bdy_interp:dt) "QICE" "Ice mixing ratio" "kg kg-1"
467	state real qs ikjftb moist 1 - \
468	i01rusdf=(bdy_interp:dt) "QSNOW" "Snow mixing ratio" "kg kg-1"
469	state real qg ikjftb moist 1 - \
470	i01rusdf=(bdy_interp:dt) "QGRAUP" "Graupel mixing ratio" "kg kg-1"
471
472	# Chem Scalars
473	state real - ikjftb chem 1 - - -
474
475	# Other Scalars: see also above
476	state real qni ikjftb scalar 1 - \
477	i01rusdf=(bdy_interp:dt) "QNICE" "Ice Number concentration" "# kg(-1)"
478	state real qt ikjftb scalar 1 - \
479	i01rusdf=(bdy_interp:dt) "CWM" "Total condensate mixing ratio" "kg kg-1"
480
481	#-----------------------------------------------------------------------------------------------------------------------------------------------------------------
482
483	# Arrays for Specified LBCs (lbc arrays REMOVED; Boundary arrays are now specified with the state array; see above, 20050413 JM )
484
485	state real fcx w misc - - ir "fcx" "RELAXATION TERM FOR BOUNDARY ZONE" ""
486	state real gcx w misc - - ir "gcx" "2ND RELAXATION TERM FOR BOUNDARY ZONE" ""
487	state real dtbc - misc - - ir "dtbc" "TIME SINCE BOUNDARY READ" ""
488
489	#-------------------------------------------------------------------------------------------------------------------------------------------
490	#-------------------------------------------------------------------------------------------------------------------------------------------
491	# Physics Related State Varibles
492
493	#-------------------------------------------------------------------------------------------------------------------------------------------
494	# SI - start variables from netCDF format from Standard Initialization, most eventually for use in LSM schemes
495	#-------------------------------------------------------------------------------------------------------------------------------------------
496
497	state real sm000007 ij misc 1 - i1 "SM000007" "LAYER SOIL MOISTURE" "m3 m-3"
498	state real sm007028 ij misc 1 - i1 "SM007028" "LAYER SOIL MOISTURE" "m3 m-3"
499	state real sm028100 ij misc 1 - i1 "SM028100" "LAYER SOIL MOISTURE" "m3 m-3"
500	state real sm100255 ij misc 1 - i1 "SM100255" "LAYER SOIL MOISTURE" "m3 m-3"
501	state real st000007 ij misc 1 - i1 "ST000007" "LAYER SOIL TEMPERATURE" "K"
502	state real st007028 ij misc 1 - i1 "ST007028" "LAYER SOIL TEMPERATURE" "K"
503	state real st028100 ij misc 1 - i1 "ST028100" "LAYER SOIL TEMPERATURE" "K"
504	state real st100255 ij misc 1 - i1 "ST100255" "LAYER SOIL TEMPERATURE" "K"
505	state real sm000010 ij misc 1 - i1 "SM000010" "LAYER SOIL MOISTURE" "m3 m-3"
506	state real sm010040 ij misc 1 - i1 "SM010040 " "LAYER SOIL MOISTURE" "m3 m-3"
507	state real sm040100 ij misc 1 - i1 "SM040100 " "LAYER SOIL MOISTURE" "m3 m-3"
508	state real sm100200 ij misc 1 - i1 "SM100200 " "LAYER SOIL MOISTURE" "m3 m-3"
509	state real sm010200 ij misc 1 - i1 "SM010200" "LAYER SOIL MOISTURE" "m3 m-3"
510	state real soilm000 ij misc 1 - i1 "SOILM000" "LAYER SOIL MOISTURE" "m3 m-3"
511	state real soilm005 ij misc 1 - i1 "SOILM005" "LAYER SOIL MOISTURE" "m3 m-3"
512	state real soilm020 ij misc 1 - i1 "SOILM020" "LAYER SOIL MOISTURE" "m3 m-3"
513	state real soilm040 ij misc 1 - i1 "SOILM040" "LAYER SOIL MOISTURE" "m3 m-3"
514	state real soilm160 ij misc 1 - i1 "SOILM160" "LAYER SOIL MOISTURE" "m3 m-3"
515	state real soilm300 ij misc 1 - i1 "SOILM300" "LAYER SOIL MOISTURE" "m3 m-3"
516	state real sw000010 ij misc 1 - i1 "SW000010" "LAYER SOIL LIQUID" "m3 m-3"
517	state real sw010040 ij misc 1 - i1 "SW010040" "LAYER SOIL LIQUID" "m3 m-3"
518	state real sw040100 ij misc 1 - i1 "SW040100" "LAYER SOIL LIQUID" "m3 m-3"
519	state real sw100200 ij misc 1 - i1 "SW100200" "LAYER SOIL LIQUID" "m3 m-3"
520	state real sw010200 ij misc 1 - i1 "SW010200" "LAYER SOIL LIQUID" "m3 m-3"
521	state real soilw000 ij misc 1 - i1 "SOILW000" "LAYER SOIL LIQUID" "m3 m-3"
522	state real soilw005 ij misc 1 - i1 "SOILW005" "LAYER SOIL LIQUID" "m3 m-3"
523	state real soilw020 ij misc 1 - i1 "SOILW020" "LAYER SOIL LIQUID" "m3 m-3"
524	state real soilw040 ij misc 1 - i1 "SOILW040" "LAYER SOIL LIQUID" "m3 m-3"
525	state real soilw160 ij misc 1 - i1 "SOILW160" "LAYER SOIL LIQUID" "m3 m-3"
526	state real soilw300 ij misc 1 - i1 "SOILW300" "LAYER SOIL LIQUID" "m3 m-3"
527	state real st000010 ij misc 1 - i1 "ST000010" "LAYER SOIL TEMPERATURE" "K"
528	state real st010040 ij misc 1 - i1 "ST010040" "LAYER SOIL TEMPERATURE" "K"
529	state real st040100 ij misc 1 - i1 "ST040100" "LAYER SOIL TEMPERATURE" "K"
530	state real st100200 ij misc 1 - i1 "ST100200" "LAYER SOIL TEMPERATURE" "K"
531	state real st010200 ij misc 1 - i1 "ST010200" "LAYER SOIL TEMPERATURE" "K"
532	state real soilt000 ij misc 1 - i1 "SOILT000" "LAYER SOIL TEMPERATURE" "K"
533	state real soilt005 ij misc 1 - i1 "SOILT005" "LAYER SOIL TEMPERATURE" "K"
534	state real soilt020 ij misc 1 - i1 "SOILT020" "LAYER SOIL TEMPERATURE" "K"
535	state real soilt040 ij misc 1 - i1 "SOILT040" "LAYER SOIL TEMPERATURE" "K"
536	state real soilt160 ij misc 1 - i1 "SOILT160" "LAYER SOIL TEMPERATURE" "K"
537	state real soilt300 ij misc 1 - i1 "SOILT300" "LAYER SOIL TEMPERATURE" "K"
538	state real landmask ij misc 1 - i012rd=(interp_fcnm)u=(copy_fcnm) "LANDMASK" "LAND MASK (1 FOR LAND, 0 FOR WATER)" ""
539	state real topostdv ij misc 1 - i12 "TOPOSTDV" "ELEVATION STD DEV" "m"
540	state real toposlpx ij misc 1 - i12 "TOPOSLPX" "ELEVATION X SLOPE" ""
541	state real toposlpy ij misc 1 - i12 "TOPOSLPY" "ELEVATION Y SLOPE" ""
542	state real shdmax ij misc 1 - i012r "SHDMAX" "ANNUAL MAX VEG FRACTION" ""
543	state real shdmin ij misc 1 - i012r "SHDMIN" "ANNUAL MIN VEG FRACTION" ""
544	state real snoalb ij misc 1 - i012r "SNOALB" "ANNUAL MAX SNOW ALBEDO IN FRACTION" ""
545	state real slopecat ij misc 1 - i12 "SLOPECAT" "SLOPE CATEGORY" ""
546	state real toposoil ij misc 1 - i12 "SOILHGT" "ELEVATION OF LSM DATA" "m"
547	state real landusef iuj misc 1 Z i12 "LANDUSEF" "LANDUSE FRACTION BY CATEGORY" ""
548	state real soilctop isj misc 1 Z i12 "SOILCTOP" "SOIL CAT FRACTION (TOP)" ""
549	state real soilcbot isj misc 1 Z i1 "SOILCBOT" "SOIL CAT FRACTION (BOTTOM)" ""
550	state real soilcat ij misc 1 - i12 "SOILCAT" "SOIL CAT DOMINANT TYPE" ""
551	state real vegcat ij misc 1 - i12 "VEGCAT" "VEGETATION CAT DOMINANT TYPE" ""
552
553	#---------------------------------------------------------------------------------------------------------------------------------------
554	# SI - end variables from netCDF format from Standard Initialization
555	#---------------------------------------------------------------------------------------------------------------------------------------
556
557	# soil model variables (Note that they are marked as staggered in the vertical dimension
558	# because they are "fully dimensioned" -- they use every element in that dim
559	state real TSLB ilj misc 1 Z i02rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TSLB" "SOIL TEMPERATURE" "K"
560	state real ZS l misc - Z irh "ZS" "DEPTHS OF CENTERS OF SOIL LAYERS" "m"
561	state real DZS l misc - Z irh "DZS" "THICKNESSES OF SOIL LAYERS" "m"
562
563	# urban model variables
564	state real DZR l em - Z r "DZR" "THICKNESSES OF ROOF LAYERS" "m"
565	state real DZB l em - Z r "DZB" "THICKNESSES OF WALL LAYERS" "m"
566	state real DZG l em - Z r "DZG" "THICKNESSES OF ROAD LAYERS" "m"
567
568	# lsm State Variables
569
570	#
571	# sortie desactivé, mais peut servir pour le modèle de sol
572	#
573	state real SMOIS ilj - 1 Z i02rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "SMOIS" "SOIL MOISTURE" "m3 m-3"
574	state real SH2O ilj - 1 Z i02rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "SH2O" "SOIL LIQUID WATER" "m3 m-3"
575	state real XICE ij misc 1 - i012rd=(interp_mask_water_field:lu_index)u=(copy_fcnm) "XICE" "SEA ICE FLAG" ""
576	state real SMSTAV ij misc 1 - rd=(interp_mask_land_field:lu_index) "SMSTAV" "MOISTURE AVAILABILITY" ""
577	state real SMSTOT ij misc 1 - r "SMSTOT" "TOTAL SOIL MOISTURE" "m3 m-3"
578
579	state real SFCRUNOFF ij misc 1 - rd=(interp_mask_land_field:lu_index) "SFROFF" "SURFACE RUNOFF" "mm"
580	state real UDRUNOFF ij misc 1 - rd=(interp_mask_land_field:lu_index) "UDROFF" "UNDERGROUND RUNOFF" "mm"
581	state integer IVGTYP ij misc 1 - i02rd=(interp_fcni)u=(copy_fcni) "IVGTYP" "DOMINANT VEGETATION CATEGORY" ""
582	state integer ISLTYP ij misc 1 - i02rd=(interp_fcni)u=(copy_fcni) "ISLTYP" "DOMINANT SOIL CATEGORY" ""
583	state real VEGFRA ij misc 1 - i0125rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "VEGFRA" "VEGETATION FRACTION" ""
584	state real SFCEVP ij misc 1 - r "SFCEVP" "SURFACE EVAPORATION" "kg m-2"
585	state real GRDFLX ij misc 1 - r "GRDFLX" "GROUND HEAT FLUX" "W m-2"
586	state real SFCEXC ij misc 1 - r "SFCEXC " "SURFACE EXCHANGE COEFFICIENT" "m s-1"
587
588	state real ACSNOW ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "ACSNOW" "ACCUMULATED SNOW" "kg m-2"
589	state real ACSNOM ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "ACSNOM" "ACCUMULATED MELTED SNOW" "kg m-2"
590	#
591	# ici sortie désactivée
592	#
593	state real SNOW ij misc 1 - i01rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "SNOW" "SNOW WATER EQUIVALENT" "kg m-2"
594	state real SNOWH ij misc 1 - i01rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "SNOWH" "PHYSICAL SNOW DEPTH" "m"
595	state real RHOSN ij misc 1 - i01rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "RHOSN" " SNOW DENSITY" "kg m-3"
596	state real CANWAT ij misc 1 - i01rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "CANWAT" "CANOPY WATER" "kg m-2"
597	state real SST ij misc 1 - i0125rd=(interp_mask_water_field:lu_index) "SST" "SEA SURFACE TEMPERATURE" "K"
598	state integer IFNDSNOWH - misc 1 - i "FNDSNOWH" "SNOWH_LOGICAL"
599	state integer IFNDSOILW - misc 1 - i "FNDSOILW" "SOILW_LOGICAL"
600
601	# urban state variables
602	state real TR_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TR_URB" "URBAN ROOF SKIN TEMPERATURE" "K"
603	state real TB_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TB_URB" "URBAN WALL SKIN TEMPERATURE" "K"
604	state real TG_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TG_URB" "URBAN ROAD SKIN TEMPERATURE" "K"
605	state real TC_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TC_URB" "URBAN CANOPY TEMPERATURE" "K"
606	state real QC_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "QC_URB" "URBAN CANOPY HUMIDITY" "kg kg{-1}"
607	state real UC_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "UC_URB" "URBAN CANOPY WIND" "m s{-1}"
608	state real XXXR_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "XXXR_URB" "M-O LENGTH ABOVE URBAN ROOF" "dimensionless"
609	state real XXXB_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "XXXB_URB" "M-O LENGTH ABOVE URBAN WALL" "dimensionless"
610	state real XXXG_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "XXXG_URB" "M-O LENGTH ABOVE URBAN ROAD" "dimensionless"
611	state real XXXC_URB2D ij misc 1 - rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "XXXC_URB" "M-O LENGTH ABOVE URBAN CANOPY" "dimensionless"
612	state real TRL_URB3D ilj misc 1 Z rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TRL_URB" "ROOF LAYER TEMPERATURE" "K"
613	state real TBL_URB3D ilj misc 1 Z rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TBL_URB" "WALL LAYER TEMPERATURE" "K"
614	state real TGL_URB3D ilj misc 1 Z rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TGL_URB" "ROAD LAYER TEMPERATURE" "K"
615	state real SH_URB2D ij misc 1 - r "SH_URB" "SENSIBLE HEAT FLUX FROM URBAN SFC" "W m{-2}"
616	state real LH_URB2D ij misc 1 - r "LH_URB" "LATENT HEAT FLUX FROM URBAN SFC" "W m{-2}"
617	state real G_URB2D ij misc 1 - r "G_URB" "GROUND HEAT FLUX INTO URBAN" "W m{-2}"
618	state real RN_URB2D ij misc 1 - r "RN_URB" "NET RADIATION ON URBAN SFC" "W m{-2}"
619	state real TS_URB2D ij misc 1 - r "TS_URB" "SKIN TEMPERATURE" "K"
620	state real FRC_URB2D ij misc 1 - rd=(interp_fcnm)u=(copy_fcnm) "FRC_URB" "URBAN FRACTION" "dimensionless"
621	state integer UTYPE_URB2D ij misc 1 - rd=(interp_fcnm)u=(copy_fcnm) "UTYPE_URB" "URBAN TYPE" "dimensionless"
622
623
624	# urban variables from radiation model
625	state real COSZ_URB2D ij misc 1 - r "COSZ_URB" "COS of SOLAR ZENITH ANGLE" "dimensionless"
626	state real OMG_URB2D ij misc 1 - r "OMG_URB" "SOLAR HOUR ANGLE" "dimensionless"
627	state real DECLIN_URB - misc 1 - r "DECLIN_URB" "SOLAR DECLINATION" "dimensionless"
628
629
630	# RUC LSM
631	state real SMFR3D ilj misc 1 Z r "SMFR3D" "SOIL ICE" ""
632	state real KEEPFR3DFLAG ilj misc 1 Z r "KEEPFR3DFLAG" "FLAG - 1. FROZEN SOIL YES, 0 - NO" ""
633
634	# MRF PBL variables
635	i1 real PSIM ij misc 1 - - "PSIM" "SIMILARITY FUNCTION FOR MOMENTUM" ""
636	i1 real PSIH ij misc 1 - - "PSIH" "SIMILARITY FUNCTION FOR HEAT" ""
637	i1 real WSPD ij misc 1 - - "WSPD" "Wind speed" "m s-1"
638	i1 real GZ1OZ0 ij misc 1 - - "GZ1OZ0" "LOG OF Z1 over Z0" ""
639	i1 real BR ij misc 1 - - "BR" "Bulk Richardson" ""
640
641	# MYJ PBL variables
642	state real tke_myj ikj misc 1 - r "tke_myj" "TKE FROM MELLOR-YAMADA-JANJIC" "m2 s-2"
643	state real EL_MYJ ikj misc 1 - - "el_myj" "MIXING LENGTH FROM MELLOR-YAMADA-JANJIC" "m"
644	state real EXCH_H ikj misc 1 - r "EXCH_H" "EXCHANGE COEFFICIENTS "
645	state real CT ij misc 1 - r "CT" "COUNTERGRADIENT TERM" "K"
646	state real THZ0 ij misc 1 - r "THZ0" "POTENTIAL TEMPERATURE AT ZNT" "K"
647	state real Z0 ij misc 1 - r "Z0" "Background ROUGHNESS LENGTH" "m"
648	state real QZ0 ij misc 1 - r "QZ0" "SPECIFIC HUMIDITY AT ZNT" "kg kg-1"
649	state real UZ0 ij misc 1 - r "UZ0" "U WIND COMPONENT AT ZNT" "m s-1"
650	state real VZ0 ij misc 1 - r "VZ0" "V WIND COMPONENT AT ZNT" "m s-1"
651	state real QSFC ij misc 1 - r "QSFC" "SPECIFIC HUMIDITY AT LOWER BOUNDARY" "kg kg-1"
652	state real AKHS ij misc 1 - r "AKHS" "SFC EXCH COEFF FOR HEAT" "m s-1"
653	state real AKMS ij misc 1 - r "AKMS" "SFC EXCH COEFF FOR MOMENTUM" "m s-1"
654	state integer KPBL ij misc 1 - r "KPBL" "LEVEL OF PBL TOP" ""
655	i1 real TSHLTR ij misc 1 - - "TSHLTR" "SHELTER THETA FROM MYJ" "K"
656	i1 real QSHLTR ij misc 1 - - "QSHLTR" "SHELTER SPECIFIC HUMIDITY FROM MYJ" "kg kg-1"
657	i1 real PSHLTR ij misc 1 - - "PSHLTR" "SHELTER PRESSURE FROM MYJ" "Pa"
658	i1 real TH10 ij misc 1 - - "TH10" "10-M THETA FROM MYJ" "K"
659	i1 real Q10 ij misc 1 - - "Q10" "10-M SPECIFIC HUMIDITY FROM MYJ" "kg kg-1"
660	i1 real CHKLOWQ ij misc 1 - - "CHKLOWQ" "SURFACE SATURATION FLAG" ""
661
662	# gfdl (eta) radiation State Variables
663	state real HTOP ij misc 1 - r "HTOP" "TOP OF CONVECTION LEVEL" ""
664	state real HBOT ij misc 1 - r "HBOT" "BOT OF CONVECTION LEVEL" ""
665	state real HTOPR ij misc 1 - r "HTOPR" "TOP OF CONVECTION LEVEL FOR RADIATION" ""
666	state real HBOTR ij misc 1 - r "HBOTR" "BOT OF CONVECTION LEVEL FOR RADIATION" ""
667	state real CUTOP ij misc 1 - r "CUTOP" "TOP OF CONVECTION LEVEL FROM CUMULUS PAR" ""
668	state real CUBOT ij misc 1 - r "CUBOT" "BOT OF CONVECTION LEVEL FROM CUMULUS PAR" ""
669	state REAL CUPPT ij misc 1 - r "CUPPT" "ACCUMULATED CONVECTIVE RAIN SINC LAST CALL TO THE RADIATION" ""
670	state real rswtoa ij misc 1 - i
671	state real rlwtoa ij misc 1 - i
672	state real czmean ij misc 1 - i
673	state real cfracl ij misc 1 - i
674	state real cfracm ij misc 1 - i
675	state real cfrach ij misc 1 - i
676	state real acfrst ij misc 1 - i
677	state integer ncfrst ij misc 1 - i
678	state real acfrcv ij misc 1 - i
679	state integer ncfrcv ij misc 1 - i
680
681	# cam radiation variables
682	state real - iqjf ozmixm 1 - - -
683	state real mth01 iqjf ozmixm 1 - - -
684	state real mth02 iqjf ozmixm 1 - - -
685	state real mth03 iqjf ozmixm 1 - - -
686	state real mth04 iqjf ozmixm 1 - - -
687	state real mth05 iqjf ozmixm 1 - - -
688	state real mth06 iqjf ozmixm 1 - - -
689	state real mth07 iqjf ozmixm 1 - - -
690	state real mth08 iqjf ozmixm 1 - - -
691	state real mth09 iqjf ozmixm 1 - - -
692	state real mth10 iqjf ozmixm 1 - - -
693	state real mth11 iqjf ozmixm 1 - - -
694	state real mth12 iqjf ozmixm 1 - - -
695	state real pin q misc 1 - - "PIN" "PRESSURE LEVEL OF OZONE MIXING RATIO" "millibar"
696	state real m_ps ij misc 2 - - "m_ps" "PS from MATCH on WRF grids"
697	state real - idjf aerosolc 2 - - -
698	state real SUL idjf aerosolc 2 - - "SUL" "SUL aerosol concentration"
699	state real SSLT idjf aerosolc 2 - - "SSLT" "SSLT aerosol concentration"
700	state real DUST1 idjf aerosolc 2 - - "DUST1" "DUST1 aerosol concentration"
701	state real DUST2 idjf aerosolc 2 - - "DUST2" "DUST2 aerosol concentration"
702	state real DUST3 idjf aerosolc 2 - - "DUST3" "DUST3 aerosol concentration"
703	state real DUST4 idjf aerosolc 2 - - "DUST4" "DUST4 aerosol concentration"
704	state real OCPHO idjf aerosolc 2 - - "OCPHO" "OCPHO aerosol concentration"
705	state real BCPHO idjf aerosolc 2 - - "BCPHO" "BCPHO aerosol concentration"
706	state real OCPHI idjf aerosolc 2 - - "OCPHI" "OCPHI aerosol concentration"
707	state real BCPHI idjf aerosolc 2 - - "BCPHI" "BCPHI aerosol concentration"
708	state real BG idjf aerosolc 2 - - "BG" "BG aerosol concentration"
709	state real VOLC idjf aerosolc 2 - - "VOLC" "VOLC aerosol concentration"
710	state real m_hybi d misc 1 - - "m_hybi" "MATCH hybi"
711
712	# new eta microphpysics State Variables
713	state real F_ICE_PHY ikj misc 1 - rdu "F_ICE_PHY" "FRACTION OF ICE" ""
714	state real F_RAIN_PHY ikj misc 1 - rdu "F_RAIN_PHY" "FRACTION OF RAIN " ""
715	state real F_RIMEF_PHY ikj misc 1 - rdu "F_RIMEF_PHY" "MASS RATIO OF RIMED ICE " ""
716	# Other Misc State Variables
717
718
719	state real ht ij misc 1 - i012rhdus "HGT" "Terrain Height" "m"
720	state real ht_fine ij misc 1 - - "HGT_FINE" "Fine Terrain Height" "m"
721	state real ht_int ij misc 1 - - "HGT_INT" "Terrain Height Horizontally Interpolated" "m"
722	state real ht_input ij misc 1 - - "HGT_INPUT" "Terrain Height from FG Input File" "m"
723
724	#
725	# h is removed from TSK, parce que la temperature de surface dans TSK n'est pas modifiee
726	#
727	state real TSK ij misc 1 - i012rdu=(copy_fcnm) "TSK" "SURFACE SKIN TEMPERATURE" "K"
728	state real TSK_SAVE ij misc 1 - - "TSK_SAVE" "SURFACE SKIN TEMPERATURE, EXTRA COPY FOR SEA ICE TESTS in REAL" "K"
729
730	state real h_diabatic ikj misc 1 - r "h_diabatic" "PREVIOUS TIMESTEP CONDENSATIONAL HEATING" ""
731	state real msft ij misc 1 - i012rhdu=(copy_fcnm) "MAPFAC_M" "Map scale factor on mass grid" ""
732	state real msfu ij misc 1 X i012rhdu=(copy_fcnm) "MAPFAC_U" "Map scale factor on u-grid" ""
733	state real msfv ij misc 1 Y i012rhdu=(copy_fcnm) "MAPFAC_V" "Map scale factor on v-grid" ""
734	state real f ij misc 1 - i012rhdu=(copy_fcnm) "f" "Coriolis sine latitude term" "s-1"
735	state real e ij misc 1 - i012rhdu=(copy_fcnm) "e" "Coriolis cosine latitude term" "s-1"
736	state real sina ij misc 1 - i012rhdu=(copy_fcnm) "SINALPHA" "Local sine of map rotation" ""
737	state real cosa ij misc 1 - i012rhdu=(copy_fcnm) "COSALPHA" "Local cosine of map rotation" ""
738
739	state real u_base k misc 1 - ir "u_base" "BASE STATE X WIND IN IDEALIZED CASES" ""
740	state real v_base k misc 1 - ir "v_base" "BASE STATE Y WIND IN IDEALIZED CASES" ""
741	state real qv_base k misc 1 - ir "qv_base" "BASE STATE QV IN IDEALIZED CASES" ""
742	state real z_base k misc 1 - ir "z_base" "BASE STATE HEIGHT IN IDEALIZED CASES" ""
743	state real u_frame - misc 1 - ir "u_frame" "FRAME X WIND" "m s-1"
744	state real v_frame - misc 1 - ir "v_frame" "FRAME Y WIND" "m s-1"
745	# p_top appears as metadata between SI and real but as a state variable in real and WRF
746	# since it is a scalar and a constant, it makes sense to have it as metadata -- there
747	# are, however, probably post-processing programs that expect to see it as an I/O record
748	# another problem: share/input_wrf tries to read this as metadata (fine for real reading
749	# SI, but with model reading real output, it generates a warning when debug is > 0 in
750	# namelist and causes repeated questions from users. A third problem is the potential
751	# collision between a metadata name and a field record in the I/O data
752	# resolve this how? Have the real program throw a switch to tell the code to get it
753	# from the metadata? Otherwise it's a field?
754	state real p_top - misc - - irh "p_top" "PRESSURE TOP OF THE MODEL" "Pa"
755	state real lat_ll_t - dyn_em - - irh "lat_ll_t" "latitude lower left, temp point" "degrees"
756	state real lat_ul_t - dyn_em - - irh "lat_ul_t" "latitude up left, temp point" "degrees"
757	state real lat_ur_t - dyn_em - - irh "lat_ur_t" "latitude up right, temp point" "degrees"
758	state real lat_lr_t - dyn_em - - irh "lat_lr_t" "latitude lower right, temp point" "degrees"
759	state real lat_ll_u - dyn_em - - irh "lat_ll_u" "latitude lower left, u point" "degrees"
760	state real lat_ul_u - dyn_em - - irh "lat_ul_u" "latitude up left, u point" "degrees"
761	state real lat_ur_u - dyn_em - - irh "lat_ur_u" "latitude up right, u point" "degrees"
762	state real lat_lr_u - dyn_em - - irh "lat_lr_u" "latitude lower right, u point" "degrees"
763	state real lat_ll_v - dyn_em - - irh "lat_ll_v" "latitude lower left, v point" "degrees"
764	state real lat_ul_v - dyn_em - - irh "lat_ul_v" "latitude up left, v point" "degrees"
765	state real lat_ur_v - dyn_em - - irh "lat_ur_v" "latitude up right, v point" "degrees"
766	state real lat_lr_v - dyn_em - - irh "lat_lr_v" "latitude lower right, v point" "degrees"
767	state real lat_ll_d - dyn_em - - irh "lat_ll_d" "latitude lower left, massless point" "degrees"
768	state real lat_ul_d - dyn_em - - irh "lat_ul_d" "latitude up left, massless point" "degrees"
769	state real lat_ur_d - dyn_em - - irh "lat_ur_d" "latitude up right, massless point" "degrees"
770	state real lat_lr_d - dyn_em - - irh "lat_lr_d" "latitude lower right, massless point" "degrees"
771	state real lon_ll_t - dyn_em - - irh "lon_ll_t" "longitude lower left, temp point" "degrees"
772	state real lon_ul_t - dyn_em - - irh "lon_ul_t" "longitude up left, temp point" "degrees"
773	state real lon_ur_t - dyn_em - - irh "lon_ur_t" "longitude up right, temp point" "degrees"
774	state real lon_lr_t - dyn_em - - irh "lon_lr_t" "longitude lower right, temp point" "degrees"
775	state real lon_ll_u - dyn_em - - irh "lon_ll_u" "longitude lower left, u point" "degrees"
776	state real lon_ul_u - dyn_em - - irh "lon_ul_u" "longitude up left, u point" "degrees"
777	state real lon_ur_u - dyn_em - - irh "lon_ur_u" "longitude up right, u point" "degrees"
778	state real lon_lr_u - dyn_em - - irh "lon_lr_u" "longitude lower right, u point" "degrees"
779	state real lon_ll_v - dyn_em - - irh "lon_ll_v" "longitude lower left, v point" "degrees"
780	state real lon_ul_v - dyn_em - - irh "lon_ul_v" "longitude up left, v point" "degrees"
781	state real lon_ur_v - dyn_em - - irh "lon_ur_v" "longitude up right, v point" "degrees"
782	state real lon_lr_v - dyn_em - - irh "lon_lr_v" "longitude lower right, v point" "degrees"
783	state real lon_ll_d - dyn_em - - irh "lon_ll_d" "longitude lower left, massless point" "degrees"
784	state real lon_ul_d - dyn_em - - irh "lon_ul_d" "longitude up left, massless point" "degrees"
785	state real lon_ur_d - dyn_em - - irh "lon_ur_d" "longitude up right, massless point" "degrees"
786	state real lon_lr_d - dyn_em - - irh "lon_lr_d" "longitude lower right, massless point" "degrees"
787
788	# Other physics variables
789
790	state real RTHCUTEN ikj misc 1 - r "RTHCUTEN" "COUPLED THETA TENDENCY DUE TO CUMULUS SCHEME" "Pa K s-1"
791	state real RQVCUTEN ikj misc 1 - r "RQVCUTEN" "COUPLED Q_V TENDENCY DUE TO CUMULUS SCHEME" "Pa kg kg-1 s-1"
792	state real RQRCUTEN ikj misc 1 - r "RQRCUTEN" "COUPLED Q_R TENDENCY DUE TO CUMULUS SCHEME" "Pa kg kg-1 s-1"
793	state real RQCCUTEN ikj misc 1 - r "RQCCUTEN" "COUPLED Q_C TENDENCY DUE TO CUMULUS SCHEME" "Pa kg kg-1 s-1"
794	state real RQSCUTEN ikj misc 1 - r "RQSCUTEN" "COUPLED Q_S TENDENCY DUE TO CUMULUS SCHEME" "Pa kg kg-1 s-1"
795	state real RQICUTEN ikj misc 1 - r "RQICUTEN" "COUPLED Q_I TENDENCY DUE TO CUMULUS SCHEME" "Pa kg kg-1 s-1"
796	state real W0AVG ikj misc 1 - r "W0AVG" "AVERAGE VERTICAL VELOCITY FOR KF CUMULUS SCHEME" "m s-1"
797
798	#
799	# Variables 2D ne servent pas - peuvent servir pour la microphysique Mars
800	# rajouter 'h' si l'on veut des sorties
801	#
802	state real RAINC ij misc 1 - rdu "RAINC" "ACCUMULATED TOTAL CUMULUS PRECIPITATION" "mm"
803	state real RAINNC ij misc 1 - rdu "RAINNC" "ACCUMULATED TOTAL GRID SCALE PRECIPITATION" "mm"
804	state real RAINCV ij misc 1 - r "RAINCV" "TIME-STEP CUMULUS PRECIPITATION" "mm"
805	state real RAINNCV ij misc 1 - r "RAINNCV" "TIME-STEP NONCONVECTIVE PRECIPITATION" "mm"
806	state real RAINBL ij misc 1 - r "RAINBL" "PBL TIME-STEP TOTAL PRECIPITATION" "mm"
807	state real SNOWNC ij misc 1 - rdu "SNOWNC" "ACCUMULATED TOTAL GRID SCALE SNOW AND ICE" "mm"
808	state real GRAUPELNC ij misc 1 - rdu "GRAUPELNC" "ACCUMULATED TOTAL GRID SCALE GRAUPEL" "mm"
809	state real SNOWNCV ij misc 1 - r "SNOWNCV" "TIME-STEP NONCONVECTIVE SNOW AND ICE" "mm"
810	state real GRAUPELNCV ij misc 1 - r "GRAUPELNCV" "TIME-STEP NONCONVECTIVE GRAUPEL" "mm"
811	state real NCA ij misc 1 - r "NCA" "COUNTER OF THE CLOUD RELAXATION TIME IN KF CUMULUS SCHEME" ""
812	state integer LOWLYR ij misc 1 - - "LOWLYR" "INDEX OF LOWEST MODEL LAYER ABOVE THE GROUND IN BMJ SCHEME" ""
813	state real MASS_FLUX ij misc 1 - r "MASS_FLUX" "DOWNDRAFT MASS FLUX FOR IN GRELL CUMULUS SCHEME" "mb hour-1"
814	state real apr_gr ij misc 1 - r "APR_GR" "PRECIP FROM CLOSURE OLD_GRELL" "mm hour-1"
815	state real apr_w ij misc 1 - r "APR_W" "PRECIP FROM CLOSURE W" "mm hour-1"
816	state real apr_mc ij misc 1 - r "APR_MC" "PRECIP FROM CLOSURE KRISH MV" "mm hour-1"
817	state real apr_st ij misc 1 - r "APR_ST" "PRECIP FROM CLOSURE STABILITY" "mm hour-1"
818	state real apr_as ij misc 1 - r "APR_AS" "PRECIP FROM CLOSURE AS-TYPE" "mm hour-1"
819	state real apr_capma ij misc 1 - r "APR_CAPMA" "PRECIP FROM MAX CAP" "mm hour-1"
820	state real apr_capme ij misc 1 - r "APR_CAPME" "PRECIP FROM MEAN CAP" "mm hour-1"
821	state real apr_capmi ij misc 1 - r "APR_CAPMI" "PRECIP FROM MIN CAP" "mm hour-1"
822	state real xf_ens ije misc 1 Z r "XF_ENS" "MASS FLUX PDF IN GRELL CUMULUS SCHEME" "mb hour-1"
823	state real pr_ens ije misc 1 Z r "PR_ENS" "PRECIP RATE PDF IN GRELL CUMULUS SCHEME" "mb hour-1"
824	state real RTHFTEN ikj misc 1 - r "RTHFTEN" "TEMPERATURE TENDENCY USED IN GRELL CUMULUS SCHEME" "K s-1"
825	state real RQVFTEN ikj misc 1 - r "RQVFTEN" "MOISTURE TENDENCY USED IN GRELL CUMULUS SCHEME" "kg s-1"
826
827	state integer STEPCU - misc 1 - r "STEPCU" "NUMBER OF FUNDAMENTAL TIMESTEPS BETWEEN CONVECTION CALLS" ""
828
829	state real RTHRATEN ikj misc 1 - rd "RTHRATEN" "COUPLED THETA TENDENCY DUE TO RADIATION" "Pa K s-1"
830	state real RTHRATENLW ikj misc 1 - r "RTHRATLW" "COUPLED THETA TENDENCY DUE TO LONG WAVE RADIATION" "Pa K s-1"
831	state real RTHRATENSW ikj misc 1 - r "RTHRATSW" "COUPLED THETA TENDENCY DUE TO SHORT WAVE RADIATION" "Pa K s-1"
832	state real CLDFRA ikj misc 1 - r "CLDFRA" "CLOUD FRACTION" ""
833	state real GLW ij misc 1 - rd "GLW" "DOWNWARD LONG WAVE FLUX AT GROUND SURFACE" "W m-2"
834	state real SWDOWN ij misc 1 - rd "SWDOWN" "DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE" "W m-2"
835	state real SWDOWNC ij misc 1 - - "SWDOWNC" "DOWNWARD CLEAR-SKY SHORT WAVE FLUX AT GROUND SURFACE" "W m-2"
836	state real GSW ij misc 1 - rd "GSW" "NET SHORT WAVE FLUX AT GROUND SURFACE" "W m-2"
837
838
839	# upward and downward clearsky and total diagnostic fluxes for CAM radiation
840	#state real ACSWUPT ij misc 1 - rhdu "ACSWUPT" "ACCUMULATED UPWELLING SHORTWAVE FLUX AT TOP" "J m-2"
841	#state real ACSWUPTC ij misc 1 - rhdu "ACSWUPTC" "ACCUMULATED UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP" "J m-2"
842	#state real ACSWDNT ij misc 1 - rhdu "ACSWDNT" "ACCUMULATED DOWNWELLING SHORTWAVE FLUX AT TOP" "J m-2"
843	#state real ACSWDNTC ij misc 1 - rhdu "ACSWDNTC" "ACCUMULATED DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP" "J m-2"
844	#state real ACSWUPB ij misc 1 - rhdu "ACSWUPB" "ACCUMULATED UPWELLING SHORTWAVE FLUX AT BOTTOM" "J m-2"
845	#state real ACSWUPBC ij misc 1 - rhdu "ACSWUPBC" "ACCUMULATED UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM" "J m-2"
846	#state real ACSWDNB ij misc 1 - rhdu "ACSWDNB" "ACCUMULATED DOWNWELLING SHORTWAVE FLUX AT BOTTOM" "J m-2"
847	#state real ACSWDNBC ij misc 1 - rhdu "ACSWDNBC" "ACCUMULATED DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM" "J m-2"
848	#state real ACLWUPT ij misc 1 - rhdu "ACLWUPT" "ACCUMULATED UPWELLING LONGWAVE FLUX AT TOP" "J m-2"
849	#state real ACLWUPTC ij misc 1 - rhdu "ACLWUPTC" "ACCUMULATED UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP" "J m-2"
850	#state real ACLWDNT ij misc 1 - rhdu "ACLWDNT" "ACCUMULATED DOWNWELLING LONGWAVE FLUX AT TOP" "J m-2"
851	#state real ACLWDNTC ij misc 1 - rhdu "ACLWDNTC" "ACCUMULATED DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP" "J m-2"
852	#state real ACLWUPB ij misc 1 - rhdu "ACLWUPB" "ACCUMULATED UPWELLING LONGWAVE FLUX AT BOTTOM" "J m-2"
853	#state real ACLWUPBC ij misc 1 - rhdu "ACLWUPBC" "ACCUMULATED UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM" "J m-2"
854	#state real ACLWDNB ij misc 1 - rhdu "ACLWDNB" "ACCUMULATED DOWNWELLING LONGWAVE FLUX AT BOTTOM" "J m-2"
855	#state real ACLWDNBC ij misc 1 - rhdu "ACLWDNBC" "ACCUMULATED DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM" "J m-2"
856	#state real SWUPT ij misc 1 - rhdu "SWUPT" "INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT TOP" "J m-2"
857	#state real SWUPTC ij misc 1 - rhdu "SWUPTC" "INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP" "J m-2"
858	#state real SWDNT ij misc 1 - rhdu "SWDNT" "INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT TOP" "J m-2"
859	#state real SWDNTC ij misc 1 - rhdu "SWDNTC" "INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP" "J m-2"
860	#state real SWUPB ij misc 1 - rhdu "SWUPB" "INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT BOTTOM" "J m-2"
861	#state real SWUPBC ij misc 1 - rhdu "SWUPBC" "INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM" "J m-2"
862	#state real SWDNB ij misc 1 - rhdu "SWDNB" "INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT BOTTOM" "J m-2"
863	#state real SWDNBC ij misc 1 - rhdu "SWDNBC" "INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM" "J m-2"
864	#state real LWUPT ij misc 1 - rhdu "LWUPT" "INSTANTANEOUS UPWELLING LONGWAVE FLUX AT TOP" "J m-2"
865	#state real LWUPTC ij misc 1 - rhdu "LWUPTC" "INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP" "J m-2"
866	#state real LWDNT ij misc 1 - rhdu "LWDNT" "INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT TOP" "J m-2"
867	#state real LWDNTC ij misc 1 - rhdu "LWDNTC" "INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP" "J m-2"
868	#state real LWUPB ij misc 1 - rhdu "LWUPB" "INSTANTANEOUS UPWELLING LONGWAVE FLUX AT BOTTOM" "J m-2"
869	#state real LWUPBC ij misc 1 - rhdu "LWUPBC" "INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM" "J m-2"
870	#state real LWDNB ij misc 1 - rhdu "LWDNB" "INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT BOTTOM" "J m-2"
871	#state real LWDNBC ij misc 1 - rhdu "LWDNBC" "INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM" "J m-2"
872
873	state real SWCF ij misc 1 - r "SWCF" "SHORT WAVE CLOUD FORCING AT TOA" "W m-2"
874	state real LWCF ij misc 1 - r "LWCF" "LONG WAVE CLOUD FORCING AT TOA" "W m-2"
875	state real OLR ij misc 1 - r "OLR" "TOA OUTGOING LONG WAVE" "W m-2"
876
877	# these next 2 are for the HFSoLE/PET demo; writing these to auxhist1 output over MCEL for coupling
878	# with wave model, only if compiled with -DMCELIO, JM 2003/05/29
879	state real XLAT ij misc 1 - i0123rh01du=(copy_fcnm) "XLAT" "LATITUDE, SOUTH IS NEGATIVE" "degree_north"
880	state real XLONG ij misc 1 - i0123rh01du=(copy_fcnm) "XLONG" "LONGITUDE, WEST IS NEGATIVE" "degree_east"
881	state real XLAT_U ij dyn_em 1 X i012rh01du=(copy_fcnm) "XLAT_U" "LATITUDE, SOUTH IS NEGATIVE" "degree_north"
882	state real XLONG_U ij dyn_em 1 X i012rh01du=(copy_fcnm) "XLONG_U" "LONGITUDE, WEST IS NEGATIVE" "degree_east"
883	state real XLAT_V ij dyn_em 1 Y i012rh01du=(copy_fcnm) "XLAT_V" "LATITUDE, SOUTH IS NEGATIVE" "degree_north"
884	state real XLONG_V ij dyn_em 1 Y i012rh01du=(copy_fcnm) "XLONG_V" "LONGITUDE, WEST IS NEGATIVE" "degree_east"
885	state real ALBEDO ij misc 1 - r "ALBEDO" "ALBEDO"
886	state real ALBBCK ij misc 1 - i012r "ALBBCK" "BACKGROUND ALBEDO" ""
887	state real EMISS ij misc 1 - r "EMISS" "SURFACE EMISSIVITY" ""
888
889	state real CLDEFI ij misc 1 - r "CLDEFI" "precipitation efficiency in BMJ SCHEME" ""
890	state integer STEPRA - misc 1 - r "STEPRA" "NUMBER OF FUNDAMENTAL TIMESTEPS BETWEEN RADIATION CALLS" ""
891
892	state real RUBLTEN ikj misc 1 - r "RUBLTEN" "COUPLED X WIND TENDENCY DUE TO PBL PARAMETERIZATION" "Pa m s-2"
893	state real RVBLTEN ikj misc 1 - r "RVBLTEN" "COUPLED Y WIND TENDENCY DUE TO PBL PARAMETERIZATION" "Pa m s-2"
894	state real RTHBLTEN ikj misc 1 - r "RTHBLTEN" "COUPLED THETA TENDENCY DUE TO PBL PARAMETERIZATION" "Pa K s-1"
895	state real RQVBLTEN ikj misc 1 - r "RQVBLTEN" "COUPLED Q_V TENDENCY DUE TO PBL PARAMETERIZATION" "Pa kg kg-1 s-1"
896	state real RQCBLTEN ikj misc 1 - r "RQCBLTEN" "COUPLED Q_C TENDENCY DUE TO PBL PARAMETERIZATION" "Pa kg kg-1 s-1"
897	state real RQIBLTEN ikj misc 1 - r "RQIBLTEN" "COUPLED Q_I TENDENCY DUE TO PBL PARAMETERIZATION" "Pa kg kg-1 s-1"
898
899	# State vector for etampnew microphysics. Must be declared state because it is not read-once and is needed for restarting.
900	state real mp_restart_state p misc 1 - r "MP_RESTART_STATE" "STATE VECTOR FOR MICROPHYSICS RESTARTS"
901	state real tbpvs_state p misc 1 - r "TBPVS_STATE" "STATE FOR ETAMPNEW MICROPHYSICS"
902	state real tbpvs0_state p misc 1 - r "TBPVS0_STATE" "STATE FOR ETAMPNEW MICROPHYSICS"
903
904	# State variables for landuse_init, Must be declared state because they are read in and needed for restarts. Had been SAVE vars in
905	# landuse_init (phys/module_physics_init.F)
906	state integer landuse_isice - misc - - r
907	state integer landuse_lucats - misc - - r
908	state integer landuse_luseas - misc - - r
909	state integer landuse_isn - misc - - r
910	state real lu_state p misc - - r
911
912	i1 real th_phy ikj misc 1 -
913	i1 real pi_phy ikj misc 1 -
914	i1 real p_phy ikj misc 1 -
915	i1 real t_phy ikj misc 1 -
916	i1 real u_phy ikj misc 1 -
917	i1 real v_phy ikj misc 1 -
918	i1 real w_phy ikj misc 1 -
919	i1 real dz8w ikj misc 1 Z
920	i1 real p8w ikj misc 1 Z
921	i1 real t8w ikj misc 1 Z
922	i1 real rho_phy ikj misc 1 -
923	i1 logical CU_ACT_FLAG ij misc 1 -
924
925
926	state real TMN ij misc 1 - i012rd=(interp_mask_land_field:lu_index)u=(copy_fcnm) "TMN" "SOIL TEMPERATURE AT LOWER BOUNDARY" "K"
927	state real XLAND ij misc 1 - i02rd=(interp_fcnm)u=(copy_fcnm) "XLAND" "LAND MASK (1 FOR LAND, 2 FOR WATER)" ""
928	state real ZNT ij misc 1 - i3r "ZNT" "TIME-VARYING ROUGHNESS LENGTH" "m"
929	##### LES mettre h pour les deux variables ci-dessous
930	state real UST ij misc 1 - r "UST" "U* IN SIMILARITY THEORY" "m s-1"
931	state real USTM ij misc 1 - rh "USTM" "U* IN SIMILARITY THEORY - LES MODE" "m s-1" #SAVEMARS2 ustar
932	i1 real HOL ij misc 1 - - "HOL" "PBL HEIGHT OVER MONIN-OBUKHOV LENGTH" ""
933	state real RMOL ij misc 1 - r "RMOL" "1./Monin Ob. Length" ""
934	state real MOL ij misc 1 - r "MOL" "T* IN SIMILARITY THEORY" "K"
935	state real PBLH ij misc 1 - r "PBLH" "PBL HEIGHT" "m"
936	state real CAPG ij misc 1 - r "CAPG" "HEAT CAPACITY FOR SOIL" "J K-1 m-3"
937	state real THC ij misc 1 - r "THC" "THERMAL INERTIA" "Cal cm-1 K-1 s-0.5"
938	state real HFX ij misc 1 - rh "HFX" "UPWARD HEAT FLUX AT THE SURFACE" "W m-2" #SAVEMARS2 sensheat
939	state real QFX ij misc 1 - r "QFX" "UPWARD MOISTURE FLUX AT THE SURFACE" "kg m-2 s-1"
940	state real LH ij misc 1 - r "LH" "LATENT HEAT FLUX AT THE SURFACE" "W m-2"
941	state real FLHC ij misc 1 - r "FLHC" "SURFACE EXCHANGE COEFFICIENT FOR HEAT" ""
942	state real FLQC ij misc 1 - r "FLQC" "SURFACE EXCHANGE COEFFICIENT FOR MOISTURE" ""
943	state real QSG ij misc 1 - r "QSG" "SURFACE SATURATION WATER VAPOR MIXING RATIO" "kg kg-1"
944	state real QVG ij misc 1 - r "QVG" "WATER VAPOR MIXING RATIO AT THE SURFACE" "kg kg-1"
945	state real QCG ij misc 1 - r "QCG" "CLOUD WATER MIXING RATIO AT THE SURFACE" "kg kg-1"
946	state real SOILT1 ij misc 1 - r "SOILT1" "TEMPERATURE INSIDE SNOW " "K"
947	state real TSNAV ij misc 1 - r "TSNAV" "AVERAGE SNOW TEMPERATURE " "C"
948	i1 real REGIME ij misc 1 -
949	state real SNOWC ij misc 1 - ird=(interp_mask_land_field:lu_index)u=(copy_fcnm) "SNOWC" "FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER)" ""
950	state real MAVAIL ij misc 1 - r "MAVAIL" "SURFACE MOISTURE AVAILABILITY" ""
951
952	state real tkesfcf ij misc 1 - r "tkesfcf" "TKE AT THE SURFACE" "m2 s-2"
953
954	state integer STEPBL - misc 1 - r "STEPBL" "NUMBER OF FUNDAMENTAL TIMESTEPS BETWEEN PBL CALLS" ""
955	state real taucldi ikj misc 1 - r "TAUCLDI" "CLOUD OPTICAL THICKNESS FOR ICE" ""
956	state real taucldc ikj misc 1 - r "TAUCLDC" "CLOUD OPTICAL THICKNESS FOR WATER" ""
957
958	state real defor11 ikj misc 1 - r "defor11" "DEFORMATION 11" "s-1"
959	state real defor22 ikj misc 1 - r "defor22" "DEFORMATION 22" "s-1"
960	state real defor12 ikj misc 1 - r "defor12" "DEFORMATION 12" "s-1"
961	state real defor33 ikj misc 1 z r "defor33" "DEFORMATION 33" "s-1"
962	state real defor13 ikj misc 1 z r "defor13" "DEFORMATION 13" "s-1"
963	state real defor23 ikj misc 1 z r "defor23" "DEFORMATION 23" "s-1"
964	state real xkmv ikj misc 1 - r "xkmv" "VERTICAL EDDY VISCOSITY" "m2 s-1"
965	state real xkmh ikj misc 1 - r "xkmh" "HORIZONTAL EDDY VISCOSITY" "m2 s-1"
966	state real xkmhd ikj misc 1 - r "xkmhd" "HORIZONTAL EDDY DIFFUSIVITY" "m2 s-1"
967	state real xkhv ikj misc 1 - r "xkhv" "VERTICAL EDDY DIFFUSIVITY OF HEAT" "m2 s-1"
968	state real xkhh ikj misc 1 - r "xkhh" "HORIZONTAL EDDY DIFFUSIVITY OF HEAT" "m2 s-1"
969	state real div ikj misc 1 - r "div" "DIVERGENCE" "s-1"
970	state real BN2 ikj misc 1 - r "BN2" "BRUNT-VAISALA FREQUENCY" "s-2"
971	state logical warm_rain - misc 1 - - "warm_rain" "WARM_RAIN_LOGICAL"
972	state logical adv_moist_cond - misc 1 - - "adv_moist_cond" "ADVECT MOIST CONDENSATES LOGICAL"
973
974	## FDDA variables
975
976	state integer STEPFG - misc 1 - r "STEPFG" "NUMBER OF FUNDAMENTAL TIMESTEPS BETWEEN FDDA GRID CALLS" ""
977	state real RUNDGDTEN ikj misc 1 X r "RUNDGDTEN" "COUPLED X WIND TENDENCY DUE TO FDDA GRID NUDGING" "Pa m s-2"
978	state real RVNDGDTEN ikj misc 1 Y r "RVNDGDTEN" "COUPLED Y WIND TENDENCY DUE TO FDDA GRID NUDGING" "Pa m s-2"
979	state real RTHNDGDTEN ikj misc 1 - r "RTHNDGDTEN" "COUPLED THETA TENDENCY DUE TO FDDA GRID NUDGING" "Pa K s-1"
980	state real RQVNDGDTEN ikj misc 1 - r "RQVNDGDTEN" "COUPLED Q_V TENDENCY DUE TO FDDA GRID NUDGING" "Pa kg kg-1 s-1"
981	state real RMUNDGDTEN ij misc 1 - r "RMUNDGDTEN" "MU TENDENCY DUE TO FDDA GRID NUDGING" "Pa s-1"
982	state real - ikjf fdda3d 1 - - -
983	state real U_NDG_NEW ikjf fdda3d 1 X igr "U_NDG_NEW" "NEW X WIND FOR FDDA GRID NUDGING" "m s-1"
984	state real V_NDG_NEW ikjf fdda3d 1 Y igr "V_NDG_NEW" "NEW Y WIND FOR FDDA GRID NUDGING" "m s-1"
985	state real T_NDG_NEW ikjf fdda3d 1 - igr "T_NDG_NEW" "NEW PERT POT TEMP FOR FDDA GRID NUDGING" "K"
986	state real Q_NDG_NEW ikjf fdda3d 1 - igr "Q_NDG_NEW" "NEW WATER VAPOR MIX RATIO FOR FDDA GRID NUDGING" "kg/kg"
987	state real PH_NDG_NEW ikjf fdda3d 1 Z igr "PH_NDG_NEW" "NEW PERT GEOPOTENTIAL FOR FDDA GRID NUDGING" "kg/kg"
988	state real U_NDG_OLD ikjf fdda3d 1 X igr "U_NDG_OLD" "OLD X WIND FOR FDDA GRID NUDGING" "m s-1"
989	state real V_NDG_OLD ikjf fdda3d 1 Y igr "V_NDG_OLD" "OLD Y WIND FOR FDDA GRID NUDGING" "m s-1"
990	state real T_NDG_OLD ikjf fdda3d 1 - igr "T_NDG_OLD" "OLD PERT POT TEMP FOR FDDA GRID NUDGING" "K"
991	state real Q_NDG_OLD ikjf fdda3d 1 - igr "Q_NDG_OLD" "OLD WATER VAPOR MIX RATIO FOR FDDA GRID NUDGING" "kg/kg"
992	state real PH_NDG_OLD ikjf fdda3d 1 Z igr "PH_NDG_OLD" "OLD PERT GEOPOTENTIAL FOR FDDA GRID NUDGING" "kg/kg"
993	state real - ivjf fdda2d 1 Z - -
994	state real MU_NDG_NEW ivjf fdda2d 1 Z igr "MU_NDG_NEW" "NEW PERT COLUMN DRY MASS FOR FDDA GRID NUDGING" "Pa"
995	state real MU_NDG_OLD ivjf fdda2d 1 Z igr "MU_NDG_OLD" "OLD PERT COLUMN DRY MASS FOR FDDA GRID NUDGING" "Pa"
996
997
998	# flag for nest movement
999	state logical moved - misc 1 - -
1000
1001	# special cam radiation restart arrays
1002	state real abstot ikcj misc 1 Z - "" "" " "
1003	state real absnxt ikaj misc 1 - - "" "" " "
1004	state real emstot ikj misc 1 Z - "" "" " "
1005
1006	# model diagnostics
1007	state real dpsdt ij misc 1 - - "dpsdt" "surface pressure tendency" "Pa/sec"
1008	state real dmudt ij misc 1 - - "dmudt" "mu tendency" "Pa/sec"
1009	state real pk1m ij misc 1 - - "pk1m" "surface pressure at previous step" "Pa"
1010	state real mu_2m ij misc 1 - - "mu_2m" "mu_2 at previous step" "Pa"
1011
1012	#
1013	#---------------------------------------------------------------------------------------------------------------------------------------
1014	#
1015
1016	######
1017	#
1018	# Variables that are set at run-time to control configuration (namelist-settable)
1019	#
1020	#<Table> <Type> <Sym> <How set> <Nentries> <Default>
1021
1022
1023	# Time Control
1024	rconfig integer run_days namelist,time_control 1 0 irh "run_days" "NUMBER OF DAYS TO RUN"
1025	rconfig integer run_hours namelist,time_control 1 0 irh "run_hours" "NUMBER OF HOURS TO RUN"
1026	rconfig integer run_minutes namelist,time_control 1 0 irh "run_minutes" "NUMBER OF MINUTES TO RUN"
1027	rconfig integer run_seconds namelist,time_control 1 0 irh "run_seconds" "NUMBER OF SECONDS TO RUN"
1028	rconfig integer start_year namelist,time_control max_domains 1999 irh "start_year" "4 DIGIT YEAR OF START OF MODEL" "YEARS"
1029	rconfig integer start_month namelist,time_control max_domains 09 irh "start_month" "2 DIGIT MONTH OF THE YEAR OF START OF MODEL, 1-12" "MONTHS"
1030	rconfig integer start_day namelist,time_control max_domains 09 irh "start_day" "2 DIGIT DAY OF THE MONTH OF START OF MODEL, 1-31" "DAYS"
1031	rconfig integer start_hour namelist,time_control max_domains 09 irh "start_hour" "2 DIGIT HOUR OF THE DAY OF START OF MODEL, 0-23" "HOURS"
1032	rconfig integer start_minute namelist,time_control max_domains 00 irh "start_minute" "2 DIGIT MINUTE OF THE HOUR OF START OF MODEL, 0-59" "MINUTES"
1033	rconfig integer start_second namelist,time_control max_domains 00 irh "start_second" "2 DIGIT SECOND OF THE MINUTE OF START OF MODEL, 0-59" "SECONDS"
1034	rconfig integer end_year namelist,time_control max_domains 1999 irh "end_year" "4 DIGIT YEAR OF END OF MODEL" "YEARS"
1035	rconfig integer end_month namelist,time_control max_domains 09 irh "end_month" "2 DIGIT MONTH OF THE YEAR OF END OF MODEL, 1-12" "MONTHS"
1036	rconfig integer end_day namelist,time_control max_domains 09 irh "end_day" "2 DIGIT DAY OF THE MONTH OF END OF MODEL, 1-31" "DAYS"
1037	rconfig integer end_hour namelist,time_control max_domains 19 irh "end_hour" "2 DIGIT HOUR OF THE DAY OF END OF MODEL, 0-23" "HOURS"
1038	rconfig integer end_minute namelist,time_control max_domains 00 irh "end_minute" "2 DIGIT MINUTE OF THE HOUR OF END OF MODEL, 0-59" "MINUTES"
1039	rconfig integer end_second namelist,time_control max_domains 00 irh "end_second" "2 DIGIT SECOND OF THE MINUTE OF END OF MODEL, 0-59" "SECONDS"
1040	rconfig integer interval_seconds namelist,time_control 1 3700 irh "interval_seconds" "SECONDS BETWEEN ANALYSIS AND BOUNDARY PERIODS" "SECONDS"
1041	rconfig logical input_from_file namelist,time_control max_domains .true. irh "input_from_file" "T/F INPUT FOR THIS DOMAIN FROM A SEPARATE INPUT FILE" ""
1042	rconfig integer fine_input_stream namelist,time_control max_domains 0 irh "fine_input_stream" "0 THROUGH 5, WHAT INPUT STREAM IS FINE GRID IC FROM" ""
1043	rconfig logical input_from_hires namelist,time_control max_domains .false. irh "input_from_hires" "T/F INPUT FOR THIS DOMAIN FROM USGS HI RES TERRAIN" ""
1044	rconfig character rsmas_data_path namelist,time_control 1 "." - "rsmas_data_path" "" ""
1045
1046	include registry.io_boilerplate
1047
1048	rconfig integer JULYR namelist,time_control max_domains 0 h "JULYR" "" ""
1049	rconfig integer JULDAY namelist,time_control max_domains 1 h "JULDAY" "" ""
1050	rconfig real GMT namelist,time_control max_domains 0. h "GMT" "" ""
1051	rconfig character input_inname namelist,time_control 1 "wrfinput_d<domain>" - "name of input infile" "" ""
1052	rconfig character input_outname namelist,time_control 1 "wrfinput_d<domain>" - "name of input outfile" "" ""
1053	rconfig character bdy_inname namelist,time_control 1 "wrfbdy_d<domain>" - "name of boundary infile" "" ""
1054	rconfig character bdy_outname namelist,time_control 1 "wrfbdy_d<domain>" - "name of boundary outfile" "" ""
1055	rconfig character rst_inname namelist,time_control 1 "wrfrst_d<domain>_<date>" - "name of restrt infile" "" ""
1056	rconfig character rst_outname namelist,time_control 1 "wrfrst_d<domain>_<date>" - "name of restrt outfile" "" ""
1057	rconfig logical write_input namelist,time_control 1 .false. - "write input data for 3dvar etc." "" ""
1058	rconfig logical write_restart_at_0h namelist,time_control 1 .false. h "write_restart_at_0h" "" ""
1059	rconfig logical adjust_output_times namelist,time_control 1 .false. - "adjust_output_times"
1060	rconfig logical adjust_input_times namelist,time_control 1 .false. - "adjust_input_times"
1061
1062	rconfig integer diag_print namelist,time_control 1 0 - "print out time series of model diagnostics"
1063	rconfig logical nocolons namelist,time_control 1 .false. - "nocolons"
1064
1065	# Domains
1066	rconfig integer time_step namelist,domains 1 - ih "time_step"
1067	rconfig integer time_step_fract_num namelist,domains 1 0 ih "time_step_fract_num"
1068	rconfig integer time_step_fract_den namelist,domains 1 1 ih "time_step_fract_den"
1069	rconfig integer max_dom namelist,domains 1 1 irh "max_dom" "" ""
1070	rconfig integer s_we namelist,domains max_domains 1 irh "s_we" "" ""
1071	rconfig integer e_we namelist,domains max_domains 101 irh "e_we" "" ""
1072	rconfig integer s_sn namelist,domains max_domains 1 irh "s_sn" "" ""
1073	rconfig integer e_sn namelist,domains max_domains 101 irh "e_sn" "" ""
1074	rconfig integer s_vert namelist,domains max_domains 1 irh "s_vert" "" ""
1075	rconfig integer e_vert namelist,domains max_domains 51 irh "e_vert" "" ""
1076	rconfig integer num_metgrid_levels namelist,domains 1 26 irh "num_metgrid_levels" "" ""
1077	rconfig real p_top_requested namelist,domains 1 2 irh "p_top_requested" "Pa" ""
1078	rconfig integer interp_type namelist,domains 1 2 irh "interp_type" "1=interp in pressure, 2=interp in LOG pressure" ""
1079	rconfig logical lowest_lev_from_sfc namelist,domains 1 .false. irh "lowest_lev_from_sfc" "" ""
1080	rconfig integer lagrange_order namelist,domains 1 1 irh "lagrange_order" "1=linear, 2=quadratic vertical interpolation" ""
1081	rconfig integer force_sfc_in_vinterp namelist,domains 1 8 irh "force_sfc_in_vinterp" "number of eta levels forced to use sfc in vert interp" ""
1082	rconfig real zap_close_levels namelist,domains 1 0 irh "zap_close_levels" "delta p where level is removed in vert interp" "Pa"
1083	rconfig logical sfcp_to_sfcp namelist,domains 1 .true. irh "sfcp_to_sfcp" "T/F use incoming sfc pres to compute new sfc pres" "flag"
1084	rconfig logical adjust_heights namelist,domains 1 .false. irh "adjust_heights" "T/F adjust pressure level input to match 500 mb height" "flag"
1085	rconfig real dx namelist,domains max_domains 5000 h "dx" "X HORIZONTAL RESOLUTION" "METERS"
1086	rconfig real dy namelist,domains max_domains 5000 h "dy" "Y HORIZONTAL RESOLUTION" "METERS"
1087	rconfig integer grid_id namelist,domains max_domains 1 irh "id" "" ""
1088	rconfig integer parent_id namelist,domains max_domains 0 h "parent_id" "" ""
1089	rconfig integer i_parent_start namelist,domains max_domains 0 rh "i_parent_start" "" ""
1090	rconfig integer j_parent_start namelist,domains max_domains 0 rh "j_parent_start" "" ""
1091	rconfig integer parent_grid_ratio namelist,domains max_domains 1 h "parent_grid_ratio" "" ""
1092	rconfig integer parent_time_step_ratio namelist,domains max_domains 1 h "parent_time_step_ratio" "" ""
1093	rconfig integer feedback namelist,domains 1 0 h "feedback" "" ""
1094	rconfig integer smooth_option namelist,domains 1 2 h "smooth_option" "" ""
1095	rconfig integer blend_width namelist,domains 1 5 h "blend_width" "width of cg fg terrain blended zone" ""
1096	rconfig real ztop namelist,domains max_domains 35000. h "ztop" "" ""
1097	rconfig integer moad_grid_ratio namelist,domains max_domains 1 h "moad_grid_ratio" "" ""
1098	rconfig integer moad_time_step_ratio namelist,domains max_domains 1 h "moad_time_step_ratio" "" ""
1099	rconfig integer shw namelist,domains max_domains 2 h "stencil_half_width" "HORIZONTAL INTERPOLATION STENCIL HALF-WIDTH" "GRID POINTS"
1100	rconfig integer tile_sz_x namelist,domains 1 0 - "tile_sz_x" "" ""
1101	rconfig integer tile_sz_y namelist,domains 1 0 - "tile_sz_y" "" ""
1102	rconfig integer numtiles namelist,domains 1 1 - "numtiles" "" ""
1103	rconfig integer nproc_x namelist,domains 1 -1 - "nproc_x" "-1 means not set" ""
1104	rconfig integer nproc_y namelist,domains 1 -1 - "nproc_y" "-1 means not set" ""
1105	rconfig integer irand namelist,domains 1 0 - "irand" "" ""
1106	rconfig real dt derived max_domains 2. h "dt" "TEMPORAL RESOLUTION" "SECONDS"
1107	rconfig integer num_moves namelist,domains 1 0
1108	rconfig integer vortex_interval namelist,domains max_domains 15 - "" "" "minutes"
1109	rconfig integer max_vortex_speed namelist,domains max_domains 40 - "" "" "meters per second"
1110	rconfig integer corral_dist namelist,domains max_domains 8
1111	rconfig integer move_id namelist,domains max_moves 0
1112	rconfig integer move_interval namelist,domains max_moves 999999999
1113	rconfig integer move_cd_x namelist,domains max_moves 0
1114	rconfig integer move_cd_y namelist,domains max_moves 0
1115	rconfig logical swap_x namelist,domains max_domains .false. rh "swap_x" "" ""
1116	rconfig logical swap_y namelist,domains max_domains .false. rh "swap_y" "" ""
1117	rconfig logical cycle_x namelist,domains max_domains .false. rh "cycle_x" "" ""
1118	rconfig logical cycle_y namelist,domains max_domains .false. rh "cycle_y" "" ""
1119	rconfig logical reorder_mesh namelist,domains 1 .false. rh "reorder_mesh" "" ""
1120	rconfig logical perturb_input namelist,domains 1 .false. h "" "" ""
1121	rconfig real eta_levels namelist,domains max_eta -1.
1122	rconfig real max_dz namelist,domains 1 1500.
1123
1124	# Physics
1125	#######
1126	####### MARS SPECIFIC NAMELIST ENTRIES
1127	#######
1128	rconfig integer mars namelist,physics max_domains 0 rh "mars" "model mode" ""
1129	rconfig real init_TI namelist,physics max_domains 0 irh "init_TI" "thermal inertia" ""
1130	rconfig real init_AL namelist,physics max_domains 0 irh "init_AL" "albedo" ""
1131	rconfig real init_Z0 namelist,physics max_domains 0 irh "init_Z0" "surface roughness" ""
1132	rconfig real init_U namelist,physics max_domains 0 irh "init_U" "zonal wind" ""
1133	rconfig real init_V namelist,physics max_domains 0 irh "init_V" "meridional wind" ""
1134	rconfig integer init_WX namelist,physics max_domains 0 irh "init_WX" "x coord" ""
1135	rconfig integer init_WY namelist,physics max_domains 0 irh "init_WY" "y coord" ""
1136	rconfig real init_MU namelist,physics max_domains 0 irh "init_MU" "zonal wind mult." ""
1137	rconfig real init_MV namelist,physics max_domains 0 irh "init_MV" "meridional wind mult." ""
1138	rconfig logical init_LES namelist,physics max_domains .false. irh "init_LES" "LES mode" ""
1139	rconfig logical modif_wrf namelist,physics max_domains .true. ir "modif_wrf" "true if tweaked WRF" ""
1140	#######
1141	####### MARS SPECIFIC NAMELIST ENTRIES
1142	#######
1143	rconfig integer mp_physics namelist,physics max_domains 0 rh "mp_physics" "" ""
1144	rconfig integer ra_lw_physics namelist,physics max_domains 0 rh "ra_lw_physics" "" ""
1145	rconfig integer ra_sw_physics namelist,physics max_domains 0 rh "ra_sw_physics" "" ""
1146	rconfig real radt namelist,physics max_domains 1 h "RADT" "" ""
1147	rconfig integer sf_sfclay_physics namelist,physics max_domains 0 rh "sf_sfclay_physics" "" ""
1148	rconfig integer sf_surface_physics namelist,physics max_domains 0 rh "sf_surface_physics" "" ""
1149	rconfig integer bl_pbl_physics namelist,physics max_domains 0 rh "bl_pbl_physics" "" ""
1150	rconfig real BLDT namelist,physics max_domains 0 h "BLDT" "" ""
1151	rconfig integer cu_physics namelist,physics max_domains 0 rh "cu_physics" "" ""
1152	rconfig real CUDT namelist,physics max_domains 100000 h "CUDT" "" ""
1153	rconfig real GSMDT namelist,physics max_domains 0 h "GSMDT" "" ""
1154	rconfig integer ISFFLX namelist,physics 1 1 irh "ISFFLX" "" ""
1155	rconfig integer IFSNOW namelist,physics 1 0 irh "IFSNOW" "" ""
1156	rconfig integer ICLOUD namelist,physics 1 1 irh "ICLOUD" "" ""
1157	rconfig real swrad_scat namelist,physics 1 1 irh "SWRAD_SCAT" "SCATTERING FACTOR IN SWRAD" ""
1158	rconfig integer surface_input_source namelist,physics 1 1 irh "surface_input_source" "1=static (fractional), 2=time dependent (dominant), 3=hybrid (not yet implemented)" ""
1159	rconfig integer num_soil_layers namelist,physics 1 10 irh "num_soil_layers" "" ""
1160	rconfig integer num_months namelist,physics 1 12 irh "num_months" "" ""
1161	rconfig integer maxiens namelist,physics 1 1 irh "maxiens" "" ""
1162	rconfig integer maxens namelist,physics 1 3 irh "maxens" "" ""
1163	rconfig integer maxens2 namelist,physics 1 3 irh "maxens2" "" ""
1164	rconfig integer maxens3 namelist,physics 1 16 irh "maxens3" "" ""
1165	rconfig integer ensdim namelist,physics 1 144 irh "ensdim" "" ""
1166	rconfig integer chem_opt namelist,physics max_domains 0 rh "chem_opt" "" ""
1167	rconfig integer num_land_cat namelist,physics 1 24 - "num_land_cat" "" ""
1168	rconfig integer num_soil_cat namelist,physics 1 16 - "num_soil_cat" "" ""
1169	rconfig integer mp_zero_out namelist,physics 1 0 - "mp_zero_out" "microphysics fields set to zero 0=no action taken, 1=all fields but Qv, 2=all fields including Qv" "flag"
1170	rconfig real mp_zero_out_thresh namelist,physics 1 1.e-8 - "mp_zero_out_thresh" "minimum threshold for non-Qv moist fields, below are set to zero" "kg/kg"
1171	rconfig real seaice_threshold namelist,physics 1 271 h "seaice_threshold" "tsk below which which water points are set to sea ice for slab scheme" "K"
1172	rconfig integer sst_update namelist,physics 1 0 h "sst_update" "update sst from wrflowinp file 0=no, 1=yes" ""
1173	rconfig integer ucmcall namelist,physics 1 0 h "ucmcall" "activate urban model 0=no, 1=yes" ""
1174	rconfig integer co2tf namelist,physics 1 0 - "co2tf" "GFDL radiation co2 flag" ""
1175	rconfig integer ra_call_offset namelist,physics 1 0 - "ra_call_offset" "radiation call offset in timesteps (-1=old, 0=new offset)" ""
1176	rconfig real cam_abs_freq_s namelist,physics 1 21600. - "cam_abs_freq_s" "CAM radiation frequency for clear-sky longwave calculations" "s"
1177	rconfig integer levsiz namelist,physics 1 1 - "levsiz" "Number of ozone data levels for CAM radiation (59)" ""
1178	rconfig integer paerlev namelist,physics 1 1 - "paerlev" "Number of aerosol data levels for CAM radiation (29)" ""
1179	rconfig integer cam_abs_dim1 namelist,physics 1 1 - "cam_abs_dim1" "dimension for absnxt in CAM radiation" ""
1180	rconfig integer cam_abs_dim2 namelist,physics 1 1 - "cam_abs_dim2" "dimension for abstot in CAM radiation" ""
1181	rconfig logical cu_rad_feedback namelist.physics max_domains .false. - "feedback cumulus to radiation"
1182
1183	#FDDA namelist parameters
1184	rconfig real FGDT namelist,fdda max_domains 0 h "FGDT" "" ""
1185	rconfig integer grid_fdda namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1186	rconfig integer if_no_pbl_nudging_uv namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1187	rconfig integer if_no_pbl_nudging_t namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1188	rconfig integer if_no_pbl_nudging_q namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1189	rconfig integer if_zfac_uv namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1190	rconfig integer k_zfac_uv namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1191	rconfig integer if_zfac_t namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1192	rconfig integer k_zfac_t namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1193	rconfig integer if_zfac_q namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1194	rconfig integer k_zfac_q namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1195	rconfig real guv namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1196	rconfig real gt namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1197	rconfig real gq namelist,fdda max_domains 0 rh "grid_fdda" "" ""
1198	rconfig real dtramp_min namelist,fdda 1 0 h "grid_fdda" "" ""
1199	rconfig integer if_ramping namelist,fdda 1 0 h "grid_fdda" "" ""
1200
1201	#Observational Nudging
1202	rconfig integer obs_nudge_opt namelist,fdda max_domains 0 rh "obs_nudge_opt" "Obs-nudging flag for domain" ""
1203	rconfig integer max_obs namelist,fdda 1 0 h "max_obs" "Maximum number of observations" ""
1204	rconfig integer nobs_ndg_vars namelist,fdda 1 0 h "num_ndg_vars" "Number of nudging variables" ""
1205	rconfig integer nobs_err_flds namelist,fdda 1 0 h "num_err_flds" "Number of error fields" ""
1206	rconfig real fdda_start namelist,fdda max_domains 0 rh "fdda_start" "Nudging start time for domain" "min"
1207	rconfig real fdda_end namelist,fdda max_domains 0 rh "fdda_end" "Nudging end time for domain" "min"
1208	rconfig integer obs_nudge_wind namelist,fdda max_domains 0 rh "obs_nudge_wind" "Wind-nudging flag for domain" ""
1209	rconfig real obs_coef_wind namelist,fdda max_domains 0 rh "obs_coef_wind" "Wind-nudging coeficient for domain" "s-1"
1210	rconfig integer obs_nudge_temp namelist,fdda max_domains 0 rh "obs_nudge_temp" "Temperature-nudging flag for domain" ""
1211	rconfig real obs_coef_temp namelist,fdda max_domains 0 rh "obs_coef_temp" "Temperature-nudging coef for domain" "s-1"
1212	rconfig integer obs_nudge_mois namelist,fdda max_domains 0 rh "obs_nudge_mois" "Moisture-nudging flag for domain" ""
1213	rconfig real obs_coef_mois namelist,fdda max_domains 0 rh "obs_coef_mois" "Moisture-nudging coef for domain" "s-1"
1214	rconfig integer obs_nudge_pstr namelist,fdda max_domains 0 rh "obs_nudge_pstr" "Not used" ""
1215	rconfig real obs_coef_pstr namelist,fdda max_domains 0 rh "obs_coef_pstr" "Not used" ""
1216	rconfig real obs_rinxy namelist,fdda max_domains 0 rh "obs_rinxy" "Horizontal radius of influence" "km"
1217	rconfig real obs_rinsig namelist,fdda 1 0 h "obs_rinsig" "Vertical radius of influence" "sigma"
1218	rconfig real obs_twindo namelist,fdda 1 0 h "obs_twindo" "Half-period time window for nudging" "min"
1219	rconfig integer obs_npfi namelist,fdda 1 0 h "obs_npfi" "Freq in cg timesteps for diag print" ""
1220	rconfig integer obs_ionf namelist,fdda 1 0 h "obs_ionf" "Freq in cg timesteps for obs input and error calc" ""
1221	rconfig integer obs_idynin namelist,fdda 1 0 h "obs_idynin" "Flag for dynamic initialization" ""
1222	rconfig real obs_dtramp namelist,fdda 1 0 h "obs_dtramp" "Time period for ramping (idynin)" "min"
1223	rconfig logical obs_ipf_in4dob namelist,fdda 1 .false. h "obs_ipf_in4dob" "Print obs input diagnostics" "min"
1224	rconfig logical obs_ipf_errob namelist,fdda 1 .false. h "obs_ipf_errob" "Print obs error diagnostics" "min"
1225	rconfig logical obs_ipf_nudob namelist,fdda 1 .false. h "obs_ipf_nudob" "Print obs nudge diagnostics" "min"
1226
1227	# Dynamics
1228	# dynamics option (see package definitions, below)
1229	rconfig integer dyn_opt namelist,dynamics 1 2 irh "dyn_opt" "" ""
1230	rconfig integer rk_ord namelist,dynamics 1 3 irh "rk_order" "" ""
1231	rconfig integer w_damping namelist,dynamics 1 0 irh "w_damping" "" ""
1232	# diff_opt 1=old diffusion, 2=new
1233	rconfig integer diff_opt namelist,dynamics 1 1 irh "diff_opt" "" ""
1234	# km_opt 1=old coefs, 2=tke, 3=Smagorinksy
1235	rconfig integer km_opt namelist,dynamics 1 4 irh "km_opt" "" ""
1236	rconfig integer damp_opt namelist,dynamics 1 0 irh "damp_opt" "" ""
1237	rconfig real zdamp namelist,dynamics max_domains 5000. h "zdamp" "" ""
1238	rconfig real dampcoef namelist,dynamics max_domains 0.05 h "dampcoef" "" ""
1239	rconfig real khdif namelist,dynamics max_domains 0 h "khdif" "" ""
1240	rconfig real kvdif namelist,dynamics max_domains 0 h "kvdif" "" ""
1241	rconfig real diff_6th_factor namelist,dynamics max_domains 0.2 h "diff_6th_factor" "factor that controls rate of 6th-order numerical diffusion"
1242	rconfig integer diff_6th_opt namelist,dynamics max_domains 2 irh "diff_6th_opt" "switch for 6th-order numerical diffusion"
1243	rconfig real smdiv namelist,dynamics max_domains 0.1 h "smdiv" "" ""
1244	rconfig real emdiv namelist,dynamics max_domains 0.01 h "emdiv" "" ""
1245	rconfig real epssm namelist,dynamics max_domains 0.1 h "epssm" "" ""
1246	rconfig logical non_hydrostatic namelist,dynamics max_domains .true. irh "non_hydrostatic" "" ""
1247	rconfig integer time_step_sound namelist,dynamics max_domains 6 h "time_step_sound" "" ""
1248	rconfig integer h_mom_adv_order namelist,dynamics max_domains 5 rh "h_mom_adv_order" "" ""
1249	rconfig integer v_mom_adv_order namelist,dynamics max_domains 3 rh "v_mom_adv_order" "" ""
1250	rconfig integer h_sca_adv_order namelist,dynamics max_domains 5 rh "h_sca_adv_order" "" ""
1251	rconfig integer v_sca_adv_order namelist,dynamics max_domains 3 rh "v_sca_adv_order" "" ""
1252	rconfig logical pd_moist namelist,dynamics max_domains .false. rh "pd_moist" "positive-definite RK3 transport switch" ""
1253	rconfig logical pd_chem namelist,dynamics max_domains .false. rh "pd_chem" "positive-definite RK3 transport switch" ""
1254	rconfig logical pd_scalar namelist,dynamics max_domains .true. rh "pd_scalar" "positive-definite RK3 transport switch" ""
1255	#rconfig logical pd_scalar namelist,dynamics max_domains .false. rh "pd_scalar" "positive-definite RK3 transport switch" ""
1256	rconfig logical pd_tke namelist,dynamics max_domains .false. rh "pd_tke" "positive-definite RK3 transport switch" ""
1257	rconfig logical top_radiation namelist,dynamics max_domains .false. rh "top_radiation" "" ""
1258	rconfig real mix_cr_len namelist,dynamics max_domains 200. h "mix_cr_len" "" ""
1259	rconfig real tke_upper_bound namelist,dynamics max_domains 1000. h "tke_upper_bound" "" ""
1260	rconfig real kh_tke_upper_bound namelist,dynamics max_domains 1000. h "kh_tke_upper_bound" "" ""
1261	rconfig real kv_tke_upper_bound namelist,dynamics max_domains 1000. h "kv_tke_upper_bound" "" ""
1262	rconfig real tke_drag_coefficient namelist,dynamics max_domains 0. h "tke_drag_coefficient" "" "dimensionless"
1263	rconfig real tke_heat_flux namelist,dynamics max_domains 0. h "tke_heat_flux" "" "K m s-1"
1264	rconfig logical pert_coriolis namelist,dynamics max_domains .false. irh "pert_coriolis" "" ""
1265	rconfig logical mix_full_fields namelist,dynamics max_domains .false. irh "mix_full_field" "" ""
1266	rconfig real base_pres namelist,dynamics 1 610. h "base_pres" "Mars Base state pressure, real only" "Pa"
1267	rconfig real base_temp namelist,dynamics 1 220. h "base_temp" "Mars Base state sea level temperature, real only" "K"
1268	rconfig real base_lapse namelist,dynamics 1 30. h "base_lapse" "Mars Base state temperature difference between base pres and 1/e of atm depth - do not change, real only" "K"
1269	rconfig real tiso namelist,dynamics 1 140. h "tiso" "Mars isothermal value for high alt, real only" "K"
1270
1271	# Bdy_control
1272	rconfig integer spec_bdy_width namelist,bdy_control 1 5 irh "spec_bdy_width" "" ""
1273	rconfig integer spec_zone namelist,bdy_control 1 1 irh "spec_zone" "" ""
1274	rconfig integer relax_zone namelist,bdy_control 1 4 irh "relax_zone" "" ""
1275	rconfig logical specified namelist,bdy_control max_domains .true. rh "specified" "" ""
1276	rconfig logical periodic_x namelist,bdy_control max_domains .false. rh "periodic_x" "" ""
1277	rconfig logical symmetric_xs namelist,bdy_control max_domains .false. rh "symmetric_xs" "" ""
1278	rconfig logical symmetric_xe namelist,bdy_control max_domains .false. rh "symmetric_xe" "" ""
1279	rconfig logical open_xs namelist,bdy_control max_domains .false. rh "open_xs" "" ""
1280	rconfig logical open_xe namelist,bdy_control max_domains .false. rh "open_xe" "" ""
1281	rconfig logical periodic_y namelist,bdy_control max_domains .false. rh "periodic_y" "" ""
1282	rconfig logical symmetric_ys namelist,bdy_control max_domains .false. rh "symmetric_ys" "" ""
1283	rconfig logical symmetric_ye namelist,bdy_control max_domains .false. rh "symmetric_ye" "" ""
1284	rconfig logical open_ys namelist,bdy_control max_domains .false. rh "open_ys" "" ""
1285	rconfig logical open_ye namelist,bdy_control max_domains .false. rh "open_ye" "" ""
1286	rconfig logical nested namelist,bdy_control max_domains .false. rh "nested" "" ""
1287	rconfig integer real_data_init_type namelist,bdy_control 1 1 irh "real_data_init_type" "REAL DATA INITIALIZATION OPTIONS: 1=SI, 2=MM5, 3=GENERIC" "PRE-PROCESSOR TYPES"
1288
1289	rconfig integer background_proc_id namelist,grib2 1 255 rh "background_proc_id" "Background processing id for grib2" ""
1290	rconfig integer forecast_proc_id namelist,grib2 1 255 rh "forecast_proc_id" "Analysis and forecast processing id for grib2" ""
1291	rconfig integer production_status namelist,grib2 1 255 rh "production_status" "Background processing id for grib2" ""
1292	rconfig integer compression namelist,grib2 1 40 rh "compression" "grib2 compression, 40 for JPEG2000 or 41 for PNG" ""
1293
1294	# NAMELIST DERIVED
1295	rconfig real cen_lat derived max_domains 0 - "cen_lat" "center latitude" "degrees, negative is south"
1296	rconfig real cen_lon derived max_domains 0 - "cen_lon" "central longitude" "degrees, negative is west"
1297	rconfig real truelat1 derived max_domains 0 - "true_lat1" "first standard parallel" "degrees, negative is south"
1298	rconfig real truelat2 derived max_domains 0 - "true_lat2" "second standard parallel" "degrees, negative is south"
1299	rconfig real moad_cen_lat derived max_domains 0 - "moad_cen_lat" "center latitude of the most coarse grid" "degrees, negative is south"
1300	rconfig real stand_lon derived max_domains 0 - "stand_lon" "standard longitude, parallel to j-direction, perpendicular to i-direction " "degrees, negative is west"
1301	rconfig real bdyfrq derived max_domains 0 - "bdyfrq" "lateral boundary input frequency" "seconds"
1302	rconfig integer iswater derived max_domains 0 - "iswater" "land use index of water" "index category"
1303	rconfig integer isice derived max_domains 0 - "isice" "land use index of ice" "index category"
1304	rconfig integer isurban derived max_domains 0 - "isurban" "land use index for 'urban and built-up" "index category"
1305	rconfig integer isoilwater derived max_domains 0 - "isoilwater" "land use index of water for soil" "index category"
1306	rconfig integer map_proj derived max_domains 0 - "map_proj" "domain map projection" "0=none, 1=Lambert, 2=polar, 3=Mercator"
1307	rconfig integer simulation_start_year derived 1 0 - "simulation_start_year" "start of simulation through restarts" "4-digit year"
1308	rconfig integer simulation_start_month derived 1 0 - "simulation_start_month" "start of simulation through restarts" "2-digit month"
1309	rconfig integer simulation_start_day derived 1 0 - "simulation_start_day" "start of simulation through restarts" "2-digit day"
1310	rconfig integer simulation_start_hour derived 1 0 - "simulation_start_hour" "start of simulation through restarts" "2-digit hour"
1311	rconfig integer simulation_start_minute derived 1 0 - "simulation_start_minute" "start of simulation through restarts" "2-digit minute"
1312	rconfig integer simulation_start_second derived 1 0 - "simulation_start_second" "start of simulation through restarts" "2-digit second"
1313
1314	#
1315	# Single dummy declaration to define a nodyn dyn option
1316	state integer nodyn_dummy - dyn_nodyn - - - "" "" ""
1317
1318	#
1319	#---------------------------------------------------------------------------------------------------------------------------------------
1320	# Package Declarations
1321	#
1322
1323	#key package associated package associated 4d scalars
1324	# name namelist choice state vars
1325
1326	package dyn_nodyn dyn_opt==0 - -
1327	package dyn_em dyn_opt==2 - -
1328	#package passivec1 chem_opt==0 -
1329
1330
1331	##### MARS OPTIONS
1332	##### MARS OPTIONS
1333	##### MARS OPTIONS
1334	##
1335	## CHANGING THE OPTION CHANGES TRACER REFERENCES WITHOUT ANY NEED TO RECOMPILE ENTIRELY LMD_MM
1336	## --> EVEN IF FIELDS ARE SCHEDULED FOR OUTPUTS ABOVE [WITH h], THESE ARE NOT IN wrfout* FILES
1337	## UNLESS THE CORRESPONDING mars OPTION IS SELECTED
1338	## --> HOWEVER, YOU STILL NEED TO RECOMPILE THE PHYSICS BEFORE CHANGING THE OPTION
1339	## AND PROBABLY YOU'LL ALSO HAVE TO CHANGE callphys.def
1340	## --> IF mars IS CHANGED IN THE NAMELIST, YOU HAVE TO RE-RUN real.exe
1341	##
1342	## IF YOU ADD A mars REFERENCE BELOW, YOU HAVE TO RECOMPILE ENTIRELY THE MODEL
1343	##
1344	## THE ORDER of TRACERS IN EACH OPTION IS VERY IMPORTANT
1345	##
1346	package default mars==0 - -
1347	package water mars==1 - scalar:qh2o,qh2o_ice
1348	package dust1 mars==2 - scalar:qdust
1349	package dust2eq mars==3 - scalar:qdust,qdustn
1350	package newbasis mars==10 - scalar:qco2
1351	package newwater mars==11 - scalar:qh2o,qh2o_ice,qdust,qdustn
1352	package newwatersca mars==12 - scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn
1353	package photochem mars==42 - scalar:qco2,chem_co,chem_o,chem_o1d,chem_o2,chem_o3,chem_h,chem_h2,chem_oh,chem_ho2,chem_h2o2,chem_ch4,chem_n2,chem_ar,qh2o_ice,qh2o,qdust,qdustn
1354	##### MARS OPTIONS
1355	##### MARS OPTIONS
1356	##### MARS OPTIONS
1357
1358
1359	package kesslerscheme mp_physics==1 - moist:qv,qc,qr
1360	package linscheme mp_physics==2 - moist:qv,qc,qr,qi,qs,qg
1361	package wsm3scheme mp_physics==3 - moist:qv,qc,qr
1362	package wsm5scheme mp_physics==4 - moist:qv,qc,qr,qi,qs
1363	package etampnew mp_physics==5 - moist:qv,qc,qr,qi,qs,qg;scalar:qt
1364	package wsm6scheme mp_physics==6 - moist:qv,qc,qr,qi,qs,qg
1365	package ncepcloud3 mp_physics==98 - moist:qv,qc,qr
1366	package ncepcloud5 mp_physics==99 - moist:qv,qc,qr,qi,qs
1367	package thompson mp_physics==8 - moist:qv,qc,qr,qi,qs,qg;scalar:qni
1368
1369	package rrtmscheme ra_lw_physics==1 - -
1370	package camlwscheme ra_lw_physics==3 - ozmixm:mth01,mth02,mth03,mth04,mth05,mth06,mth07,mth08,mth09,mth10,mth11,mth12;aerosolc:sul,sslt,dust1,dust2,dust3,dust4,ocpho,bcpho,ocphi,bcphi,bg,volc
1371	package gfdllwscheme ra_lw_physics==99 - -
1372
1373	package swradscheme ra_sw_physics==1 - -
1374	package gsfcswscheme ra_sw_physics==2 - -
1375	package camswscheme ra_sw_physics==3 - ozmixm:mth01,mth02,mth03,mth04,mth05,mth06,mth07,mth08,mth09,mth10,mth11,mth12;aerosolc:sul,sslt,dust1,dust2,dust3,dust4,ocpho,bcpho,ocphi,bcphi,bg,volc
1376	package gfdlswscheme ra_sw_physics==99 - -
1377
1378	package sfclayscheme sf_sfclay_physics==1 - -
1379	package myjsfcscheme sf_sfclay_physics==2 - -
1380	package gfssfcscheme sf_sfclay_physics==3 - -
1381	package slabscheme sf_surface_physics==1 - -
1382	package lsmscheme sf_surface_physics==2 - -
1383	package ruclsmscheme sf_surface_physics==3 - -
1384	package ysuscheme bl_pbl_physics==1 - -
1385	package myjpblscheme bl_pbl_physics==2 - -
1386	package gfsscheme bl_pbl_physics==3 - -
1387	package mrfscheme bl_pbl_physics==99 - -
1388
1389	package kfetascheme cu_physics==1 - -
1390	package bmjscheme cu_physics==2 - -
1391	package gdscheme cu_physics==3 - -
1392	package sasscheme cu_physics==4 - -
1393	package kfscheme cu_physics==99 - -
1394
1395	package psufddagd grid_fdda==1 - fdda3d:u_ndg_old,v_ndg_old,t_ndg_old,q_ndg_old,ph_ndg_old,u_ndg_new,v_ndg_new,t_ndg_new,q_ndg_new,ph_ndg_new;fdda2d:mu_ndg_old,mu_ndg_new
1396
1397	# only need to specify these once; not for every io_form* variable
1398	package io_intio io_form_restart==1 - -
1399	package io_netcdf io_form_restart==2 - -
1400	# Placeholders for additional packages (we can go beyond zzz
1401	# but that will entail modifying frame/module_io.F and frame/md_calls.m4)
1402	# Please note these are placeholders; HDF has not been implemented yet.
1403	package io_hdf io_form_restart==3 - -
1404	package io_phdf5 io_form_restart==4 - -
1405	package io_grib1 io_form_restart==5 - -
1406	package io_mcel io_form_restart==6 - -
1407	package io_esmf io_form_restart==7 - -
1408	package io_yyy io_form_restart==8 - -
1409	package io_zzz io_form_restart==9 - -
1410	package io_grib2 io_form_restart==10 - -
1411	package io_pnetcdf io_form_restart==11 - -
1412
1413	#---------------------------------------------------------------------------------------------------------------------------------------
1414	## communications
1415
1416	### 8. Edit the Registry file and create a halo-exchange for x_1.
1417
1418	# Halo Update Communications
1419
1420	halo HALO_EM_INIT_1 dyn_em 48:u_1,u_2,v_1,v_2,w_1,w_2,ph_1,ph_2
1421	halo HALO_EM_INIT_2 dyn_em 48:t_1,t_2,mu_1,mu_2,tke_1,tke_2,ww,phb
1422	halo HALO_EM_INIT_3 dyn_em 48:ph0,php,t_init,mub,mu0,p,al,alt,alb
1423	halo HALO_EM_INIT_4 dyn_em 48:pb,h_diabatic,msft,msfu,msfv,f,e,sina,cosa,ht,potevp,snopcx,soiltb
1424	halo HALO_EM_INIT_5 dyn_em 48:moist,chem,scalar
1425	halo HALO_EM_VINTERP_UV_1 dyn_em 8:pd_gc,pb
1426	halo HALO_EM_A dyn_em 8:ru,rv,rw,ww,php,alt,al,p,muu,muv,mut
1427	halo HALO_EM_PHYS_A dyn_em 4:u_2,v_2
1428	halo HALO_EM_PHYS_PBL dyn_em 4:rublten,rvblten
1429	halo HALO_EM_FDDA dyn_em 4:rundgdten,rvndgdten
1430	halo HALO_EM_PHYS_DIFFUSION dyn_em 4:defor11,defor22,defor12,defor13,defor23,div,xkmv,xkmh,xkmhd,xkhv,xkhh,tke_1,tke_2
1431	halo HALO_EM_TKE_ADVECT_3 dyn_em 24:tke_2
1432	halo HALO_EM_TKE_ADVECT_5 dyn_em 48:tke_2
1433	halo HALO_EM_TKE_A dyn_em 4:ph_2,phb
1434	halo HALO_EM_TKE_B dyn_em 4:z,rdz,rdzw,zx,zy
1435	halo HALO_EM_TKE_C dyn_em 8:u_2,v_2,z,zx,zy,rdz,rdzw
1436	halo HALO_EM_TKE_D dyn_em 8:defor11,defor22,defor33,defor12,defor13,defor23,div
1437	halo HALO_EM_TKE_E dyn_em 8:xkmv,xkmh,xkmhd,xkhv,xkhh,BN2,moist
1438	halo HALO_EM_TKE_3 dyn_em 24:tke_1,tke_2
1439	halo HALO_EM_TKE_5 dyn_em 48:tke_1,tke_2
1440	halo HALO_EM_TKE_7 dyn_em 80:tke_1,tke_2
1441	halo HALO_EM_TKE_F dyn_em 48:tke_1,tke_2
1442	halo HALO_EM_TKE_OLD_E_5 dyn_em 48:tke_1
1443	halo HALO_EM_TKE_OLD_E_7 dyn_em 80:tke_1
1444	halo HALO_EM_B dyn_em 4:ph_2,al,p,t_1,t_save,u_save,v_save,mu_1,mu_2,mudf,php,alt,pb
1445	halo HALO_EM_B2 dyn_em 4:ru_tend,rv_tend
1446	halo HALO_EM_C dyn_em 4:u_2,v_2
1447	halo HALO_EM_C2 dyn_em 4:ph_2,al,p,mu_2,muts,mudf
1448	halo HALO_EM_D dyn_em 24:ru_m,rv_m,ww_m,mut
1449	halo HALO_EM_D2_3 dyn_em 24:u_2,v_2,w_2,t_2,ph_2;24:moist,chem,scalar;4:mu_2,al
1450	halo HALO_EM_D2_5 dyn_em 48:u_2,v_2,w_2,t_2,ph_2;24:moist,chem,scalar;4:mu_2,al
1451	halo HALO_EM_D3_3 dyn_em 24:u_1,u_2,v_1,v_2,w_1,w_2,t_1,t_2,ph_1,ph_2,tke_1,tke_2,moist,chem,scalar;4:mu_1,mu_2
1452	halo HALO_EM_D3_5 dyn_em 48:u_1,u_2,v_1,v_2,w_1,w_2,t_1,t_2,ph_1,ph_2,tke_1,tke_2,moist,chem,scalar;4:mu_1,mu_2
1453	halo HALO_EM_E_3 dyn_em 24:u_1,u_2,v_1,v_2,w_1,w_2,t_1,t_2,ph_1,ph_2,tke_1,tke_2,;4:mu_1,mu_2
1454	halo HALO_EM_E_5 dyn_em 48:u_1,u_2,v_1,v_2,w_1,w_2,t_1,t_2,ph_1,ph_2,tke_1,tke_2,;4:mu_1,mu_2
1455	halo HALO_EM_MOIST_E_3 dyn_em 24:moist
1456	halo HALO_EM_MOIST_E_5 dyn_em 48:moist
1457	halo HALO_EM_MOIST_E_7 dyn_em 80:moist
1458	halo HALO_EM_CHEM_E_3 dyn_em 24:chem
1459	halo HALO_EM_CHEM_E_5 dyn_em 48:chem
1460	halo HALO_EM_CHEM_E_7 dyn_em 80:chem
1461	halo HALO_EM_SCALAR_E_3 dyn_em 24:scalar
1462	halo HALO_EM_SCALAR_E_5 dyn_em 48:scalar
1463	halo HALO_EM_SCALAR_E_7 dyn_em 80:scalar
1464
1465	halo HALO_EM_MOIST_OLD_E_3 dyn_em 24:moist_old
1466	halo HALO_EM_MOIST_OLD_E_5 dyn_em 48:moist_old
1467	halo HALO_EM_MOIST_OLD_E_7 dyn_em 80:moist_old
1468	halo HALO_EM_CHEM_OLD_E_3 dyn_em 24:chem_old
1469	halo HALO_EM_CHEM_OLD_E_5 dyn_em 48:chem_old
1470	halo HALO_EM_CHEM_OLD_E_7 dyn_em 80:chem_old
1471	halo HALO_EM_SCALAR_OLD_E_3 dyn_em 24:scalar_old
1472	halo HALO_EM_SCALAR_OLD_E_5 dyn_em 48:scalar_old
1473	halo HALO_EM_SCALAR_OLD_E_7 dyn_em 80:scalar_old
1474
1475	halo HALO_EM_FEEDBACK dyn_em 48:ht
1476
1477	halo HALO_EM_COUPLE_A dyn_em 24:mub,mu_1,mu_2
1478	period PERIOD_EM_COUPLE_A dyn_em 2:mub,mu_1,mu_2
1479	halo HALO_EM_COUPLE_B dyn_em 48:ph_1,ph_2,w_1,w_2,t_1,t_2,u_1,u_2,v_1,v_2,\
1480	moist,chem,scalar
1481	period PERIOD_EM_COUPLE_B dyn_em 3:ph_1,ph_2,w_1,w_2,t_1,t_2,u_1,u_2,v_1,v_2,\
1482	moist,chem,scalar
1483
1484	# For moving nests
1485	halo em_shift_halo_y dyn_em 48:imask_nostag,imask_xstag,imask_ystag,imask_xystag,u_2,v_2,t_2
1486	halo em_shift_halo_x dyn_em 48:imask_nostag,imask_xstag,imask_ystag,imask_xystag,u_2,v_2,t_2
1487
1488	# For observational nudging
1489	halo HALO_OBS_NUDGE dyn_em 24:pb,p,uratx,vratx,tratx
1490
1491	# Periodic Boundary Communications
1492
1493	period PERIOD_BDY_EM_INIT dyn_em 3:u_1,u_2,v_1,v_2,w_1,w_2,t_1,t_2,ph_1,ph_2,t_init,phb,ph0,php,pb,al,alt,alb,mu_1,mu_2,mub,mu0,ht,msft,msfu,msfv,sina,cosa,e,f
1494	period PERIOD_BDY_EM_MOIST dyn_em 3:moist
1495	period PERIOD_BDY_EM_CHEM dyn_em 3:chem
1496	period PERIOD_BDY_EM_SCALAR dyn_em 3:scalar
1497	period PERIOD_BDY_EM_MOIST2 dyn_em 3:moist
1498	period PERIOD_BDY_EM_CHEM2 dyn_em 3:chem
1499	period PERIOD_BDY_EM_SCALAR2 dyn_em 3:scalar
1500	period PERIOD_BDY_EM_A dyn_em 2:ru,rv,rw,ww,php,alt,p,muu,muv,mut,ph_2,al
1501	period PERIOD_BDY_EM_A1 dyn_em 3:rdzw,rdz,z,zx,zy
1502	period PERIOD_BDY_EM_PHY_BC dyn_em 2:rublten,rvblten,xkmh,xkmhd,xkmv,xkhh,xkhv,div,defor11,defor22,defor12,defor13,defor23,defor33,tke_2
1503	period PERIOD_BDY_EM_FDDA_BC dyn_em 2:rundgdten,rvndgdten
1504	period PERIOD_BDY_EM_B dyn_em 2:ru_tend,rv_tend,ph_2,al,p,t_1,t_save,u_save,v_save,mu_1,mu_2,mudf,php,alt,pb
1505	period PERIOD_BDY_EM_B3 dyn_em 2:ph_2,al,p,mu_2,muts,mudf
1506	period PERIOD_BDY_EM_B2 dyn_em 2:ru_tend,rv_tend
1507	period PERIOD_BDY_EM_C dyn_em 2:u_2,u_save,v_2,v_save,t_2,t_save,muv,msfv,muu,msfu
1508	period PERIOD_BDY_EM_D dyn_em 3:u_2,v_2,w_2,t_2,ph_2,mu_2,tke_2
1509	period PERIOD_BDY_EM_D3 dyn_em 3:u_1,u_2,v_1,v_2,w_1,w_2,t_1,t_2,ph_1,ph_2,tke_1,tke_2,mu_1,mu_2
1510
1511	#
1512	#swap SWAP_ETAMP_NEW dyn_em 1:dz8w,p_phy,pi_phy,rho,th_phy,moist,F_ICE_PHY,F_RAIN_PHY,F_RIMEF_PHY,RAINNC,RAINNCV,SR,LOWLYR
1513	#swap SWAP_WSM3 dyn_em 1:th_phy,moist,w_2,rho,pi_phy,p_phy,dz8w,rainnc,rainncv
1514	#cycle CYCLE_TEST dyn_em 1:xlong
1515
1516	##
1517
1518	# FDDA (Observational-nudging) Variables
1519	typedef fdob_type integer domain_tot # total number of domains to apply obs-nudging
1520	typedef fdob_type integer domain_init # domain initialization flag
1521	typedef fdob_type integer IEODI # end of obs data flag for current model step
1522	typedef fdob_type integer IWTSIG # flag for nudging on pressure surfaces
1523	typedef fdob_type integer NSTAT # number of obs stations used to nudge current model step
1524	typedef fdob_type integer KTAUR # restart model step
1525	typedef fdob_type integer SN_MAXCG # coarse domain grid dimension in south-north coordinate
1526	typedef fdob_type integer WE_MAXCG # coarse domain grid dimension in west-east coordinate
1527	typedef fdob_type integer SN_END # ending north-south grid index
1528	typedef fdob_type integer LEVIDN(max_domains) # level of nest
1529	typedef fdob_type real DS_CG # coarse domain grid size
1530	typedef fdob_type real WINDOW # time window half-period for nudging (in minutes)
1531	typedef fdob_type real RTLAST # time in hours of last obs used in current model step
1532	typedef fdob_type real DATEND # time in minutes after which data are asuumed to have ended
1533	typedef fdob_type real RINFMN # minimum radius of influence
1534	typedef fdob_type real RINFMX # maximum radius of influence
1535	typedef fdob_type real PFREE # pressure level (cb) where terrain effect becomes small
1536	typedef fdob_type real DCON # 1/DPSMX
1537	typedef fdob_type real DPSMX # max pres change (cb) allowed within infl range of surf obs
1538	typedef fdob_type real TFACI # scale factor used for ramp-down in dynamic initialization
1539	typedef fdob_type real XN # cone factor for Lambert projection
1540
1541	# table entries are of the form
1542	# <Table> <Type> <Sym> <Dims> <Use> <NumTLev> <Stagger> <IO> <DNAME> <DESCRIP> <UNITS>
1543	#Grid variables
1544	typedef fdob_type real varobs rz - 1 - - "varobs" "observational values in each variable"
1545	typedef fdob_type real errf hz - 1 - - "errf" "errors between model and obs values"
1546	typedef fdob_type real timeob z - 1 - - "timeob" "model times for each observation" "hours"
1547	typedef fdob_type real nlevs_ob z - 1 - - "nlevs_ob" "numbers of levels in sounding obs"
1548	typedef fdob_type real lev_in_ob z - 1 - - "lev_in_ob" "level in sounding-type obs"
1549	typedef fdob_type real plfo z - 1 - - "plfo" "index for type of obs-platform"
1550	typedef fdob_type real elevob z - 1 - - "elevob" "elevation of observation" "meters"
1551	typedef fdob_type real rio z - 1 - - "rio" "west-east grid coordinate"
1552	typedef fdob_type real rjo z - 1 - - "rjo" "south-north grid coordinate"
1553	typedef fdob_type real rko z - 1 - - "rko" "vertical grid coordinate"
1554
1555	state fdob_type fdob - -
1556	##

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format