Context Navigation

Registry.EM @ 55

Last change on this file since 55 was 55, checked in by aslmd, 14 years ago

LMD_MM_MARS: element manquant pour runs traceurs avec nouvelle physique
--> avec mars=0 et non actif, le modele fonctionne toujours bien [avec starts JBM et conditions de ses runs]
--> mettre les traceurs avec mars=11 [actifs ou non] semble faire crasher le modele: pourquoi ?

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq1.inc
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq2.inc
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq3.inc
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq4.inc
M 54 mars/libf/phymars/meso_physiq.F
le tableau tnom est defini dans module_lmd_driver en fonction de config_flag%mars [MARS_MODE]

et passe dans meso_physiq [pour ensuite servir dans initracer]

--> le COMMON advtrac.h est alors necessaire

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
config_flag%mars defini une serie de traceurs ordonnee dans le tableau scalar

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/share/module_model_constants.F
on passe au cp et R du GCM martien pour une complete correspondance

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq.inc
ce fichier est en fait inutile et pourrait etre supprime

M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
M 54 mesoscale/LMD_MM_MARS/SRC/WRFV2/main/real_em.F
modifications pour prise en compte des traceurs avec la nouvelle physique
--> avantage des traceurs generiques dans la physique
--> pour l'instant dans SCALAR l'indice 1 est dummy, 2/3 water vapor/ice et dernier CO2

M 54 mesoscale/LMD_MM_MARS/SRC/PREP_MARS/readmeteo_newphys.F90
M 54 mesoscale/LMD_MM_MARS/SRC/WPS/wps_mars/metgrid/METGRID.TBL.ARW_MarsBase_newphys
modifications pour possibilites de passer les tableaux necessaires pour

initialiser et guider les traceurs dans la nouvelle physique

M mars/libf/phymars/dimradmars.h
M mars/libf/phymars/callradite.F
version traceurs non actifs pour tests basiques

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

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

Download in other formats:

Original Format