Context Navigation

Registry.EM.newphys @ 1751

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

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

Download in other formats:

Original Format