[Dart-dev] [4046] DART/trunk/models/wrf/work/input.nml: Made the field names match the new defaults, updated the radar
nancy at ucar.edu
nancy at ucar.edu
Wed Sep 9 15:29:04 MDT 2009
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20090909/478d6285/attachment.html
-------------- next part --------------
Modified: DART/trunk/models/wrf/work/input.nml
===================================================================
--- DART/trunk/models/wrf/work/input.nml 2009-09-09 21:03:27 UTC (rev 4045)
+++ DART/trunk/models/wrf/work/input.nml 2009-09-09 21:29:04 UTC (rev 4046)
@@ -65,9 +65,10 @@
restart_in_file_name = 'smoother_ics',
restart_out_file_name = 'smoother_restart' /
+# cutoff is in radians; 0.05 is about 100km
&assim_tools_nml
filter_kind = 1,
- cutoff = 0.0004,
+ cutoff = 0.05,
sort_obs_inc = .false.,
spread_restoration = .false.,
sampling_error_correction = .false.,
@@ -96,33 +97,35 @@
&obs_kind_nml
assimilate_these_obs_types = 'RADIOSONDE_TEMPERATURE',
- 'METAR_U_10_METER_WIND',
- 'METAR_V_10_METER_WIND',
- 'METAR_TEMPERATURE_2_METER',
'RADIOSONDE_U_WIND_COMPONENT',
'RADIOSONDE_V_WIND_COMPONENT',
'SAT_U_WIND_COMPONENT',
'SAT_V_WIND_COMPONENT',
- 'GPSRO_REFRACTIVITY',
evaluate_these_obs_types = 'RADIOSONDE_SPECIFIC_HUMIDITY' /
# Notes for obs_def_radar_mod_nml:
-# (1) Reflectivity limit can be applied to observations and/or state (forward operator).
-# (2) Default lowest_reflectivity values DART will use (if apply_reflectivity_limit = .true.)
-# is missing_r8. If you want to use the default, delete the line of respective
-# lowest_reflectivity.
-# (3) As it is not clear how to assimilate Z (for now), "convert_to_dbz" is reset to .true.
-# even if you set it to .false. here in the namelist.
+# (1) Reflectivity limit can be applied to observations and/or forward operator.
+# (2) The default constants below match the WRF defaults. They will need to
+# be changed for other cases, depending on which microphysics scheme is used.
+#
&obs_def_radar_mod_nml
- convert_to_dbz = .true. ,
- dbz_threshold = 0.001 ,
- apply_ref_limit_to_obs = .false. ,
- reflectivity_limit_obs = 0.0 ,
- lowest_reflectivity_obs = 0.0 ,
- apply_ref_limit_to_state = .false. ,
- reflectivity_limit_state = 0.0 ,
- lowest_reflectivity_state = 0.0 /
+ apply_ref_limit_to_obs = .false.,
+ reflectivity_limit_obs = -10.0_r8,
+ lowest_reflectivity_obs = -10.0_r8,
+ apply_ref_limit_to_fwd_op = .false.,
+ reflectivity_limit_fwd_op = -10.0_r8,
+ lowest_reflectivity_fwd_op = -10.0_r8,
+ max_radial_vel_obs = 1000000,
+ allow_wet_graupel = .false.,
+ dielectric_factor = 0.224_r8,
+ n0_rain = 8.0e6_r8,
+ n0_graupel = 4.0e6_r8,
+ n0_snow = 3.0e6_r8,
+ rho_rain = 1000.0_r8,
+ rho_graupel = 400.0_r8,
+ rho_snow = 100.0_r8
+ /
&assim_model_nml
write_binary_restart_files = .true. /
@@ -134,7 +137,7 @@
# 3 = height
# set default_state_variables to .false. to use the explicit list.
-# otherwise it uses a hardcoded default list: U, V, W, PH, T, MU only.
+# otherwise it uses a hardcoded default list: U, V, W, PH, T, MU, QV only.
# other possible common fields are:
# wrf_state_variables = 'U','KIND_U_WIND_COMPONENT','TYPE_U','UPDATE','999',
# 'V','KIND_V_WIND_COMPONENT','TYPE_V','UPDATE','999',
@@ -154,26 +157,31 @@
&model_nml
default_state_variables = .true.,
- wrf_state_variables = 'U', 'KIND_U_WIND_COMPONENT', 'TYPE_U', 'UPDATE','999',
- 'V', 'KIND_V_WIND_COMPONENT', 'TYPE_V', 'UPDATE','999',
- 'W', 'KIND_VERTICAL_VELOCITY', 'TYPE_W', 'UPDATE','999',
- 'PH','KIND_GEOPOTENTIAL_HEIGHT', 'TYPE_GZ','UPDATE','999',
- 'T', 'KIND_POTENTIAL_TEMPERATURE','TYPE_T', 'UPDATE','999',
- 'MU','KIND_PRESSURE', 'TYPE_MU','UPDATE','999',
+ wrf_state_variables = 'U', 'KIND_U_WIND_COMPONENT', 'TYPE_U', 'UPDATE','999',
+ 'V', 'KIND_V_WIND_COMPONENT', 'TYPE_V', 'UPDATE','999',
+ 'W', 'KIND_VERTICAL_VELOCITY', 'TYPE_W', 'UPDATE','999',
+ 'PH', 'KIND_GEOPOTENTIAL_HEIGHT', 'TYPE_GZ','UPDATE','999',
+ 'T', 'KIND_POTENTIAL_TEMPERATURE','TYPE_T', 'UPDATE','999',
+ 'MU', 'KIND_PRESSURE', 'TYPE_MU','UPDATE','999',
+ 'QVAPOR','KIND_VAPOR_MIXING_RATIO', 'TYPE_QV','UPDATE','999',
wrf_state_bounds = 'QVAPOR','0.0','NULL','CLAMP',
- 'QRAIN','0.0','NULL','CLAMP',
+ 'QRAIN', '0.0','NULL','CLAMP',
'QCLOUD','0.0','NULL','CLAMP',
- output_state_vector = .false.,
- num_moist_vars = 3,
- num_domains = 1,
- calendar_type = 3,
- surf_obs = .true.,
+ output_state_vector = .false.,
+ num_moist_vars = 3,
+ num_domains = 1,
+ calendar_type = 3,
+ surf_obs = .true.,
assimilation_period_seconds = 21600,
adv_mod_command = './wrf.exe',
vert_localization_coord = 3,
center_search_half_length = 500000.,
center_spline_grid_scale = 10 /
+# vert_normal is amount of X equiv to 1 radian in horiz distance.
+# adjust to give a different localization in the vertical, if
+# horiz_dist_only is set to .false.. the default nlon/nlat should
+# be good for most grids; nlon must be an odd number.
&location_nml
horiz_dist_only = .true.,
vert_normalization_pressure = 6666666.7,
@@ -212,16 +220,16 @@
# Acceptable longitudes range from [ 0, Inf]
&obs_diag_nml
- obs_sequence_name = 'obs_seq.final',
- first_bin_center = 2003, 1, 1, 0, 0, 0 ,
- last_bin_center = 2003, 1, 2, 0, 0, 0 ,
- bin_separation = 0, 0, 0,12, 0, 0 ,
- bin_width = 0, 0, 0, 6, 0, 0 ,
- time_to_skip = 0, 0, 0, 0, 0, 0 ,
- max_num_bins = 1000,
+ obs_sequence_name = 'obs_seq.final',
+ first_bin_center = 2003, 1, 1, 0, 0, 0 ,
+ last_bin_center = 2003, 1, 2, 0, 0, 0 ,
+ bin_separation = 0, 0, 0,12, 0, 0 ,
+ bin_width = 0, 0, 0, 6, 0, 0 ,
+ time_to_skip = 0, 0, 0, 0, 0, 0 ,
+ max_num_bins = 1000,
rat_cri = 5000.0,
input_qc_threshold = 4.0,
- Nregions = 1,
+ Nregions = 4,
lonlim1 = 0.0, 0.0, 0.0, 235.0,
lonlim2 = 360.0, 360.0, 360.0, 295.0,
latlim1 = 20.0, -80.0, -20.0, 25.0,
@@ -250,19 +258,19 @@
&obs_sequence_tool_nml
- num_input_files = 1,
- filename_seq = 'obs_seq.out',
- filename_out = 'obs_seq.processed',
- first_obs_days = -1,
+ num_input_files = 1,
+ filename_seq = 'obs_seq.out',
+ filename_out = 'obs_seq.processed',
+ first_obs_days = -1,
first_obs_seconds = -1,
- last_obs_days = -1,
- last_obs_seconds = -1,
- obs_types = '',
- keep_types = .false.,
- print_only = .false.,
- min_lat = -90.0,
- max_lat = 90.0,
- min_lon = 0.0,
- max_lon = 360.0
+ last_obs_days = -1,
+ last_obs_seconds = -1,
+ obs_types = '',
+ keep_types = .false.,
+ print_only = .false.,
+ min_lat = -90.0,
+ max_lat = 90.0,
+ min_lon = 0.0,
+ max_lon = 360.0
/
More information about the Dart-dev
mailing list