[Wrf-users] Grell 3-D Subgrid Cloud Output

[Robert Hahn]

I'm trying to run the new Grell 3-D cumulus scheme with sub-grid cloud
information output (WRF version 3.4). Setting "cu_diag" to 1 in the
namelist should output the sub-grid cloud information (as seen by this
statement in the registry):

package g3tave cu_diag==1 -
state:GD_CLOUD,GD_CLOUD2,GD_CLDFR,GD_CLOUD_A,GD_CLOUD_B,GD_CLOUD2_A,GD_CLOUD2_B,kbcon_deep,ktop_deep,k22_deep

However, that alone did not produce output, so I have forced the model to
output the variables using the "iofields_filename" and associated text file
to output these variables. The variables appear in the wrfout files,
however they contain null (zero) values even after model spin-up.

Has anyone else encountered or solved this problem? My namelist.input is
below.

Robert

&time_control
run_days = 5,
run_hours = 12,
run_minutes = 0,
run_seconds = 0,
start_year = 2004, 2004, 2004, 2004,
start_month = 12, 12, 12, 12,
start_day = 15, 15, 15, 15,
start_hour = 12, 12, 12, 12,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2004, 2004, 2004, 2004,
end_month = 12, 12, 12, 12,
end_day = 21, 21, 21, 21,
end_hour = 00, 00, 00, 00,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00. 00,
interval_seconds = 10800
input_from_file = .true.,.true.,.true.,.true.
history_interval = 60, 60, 60, 60,
frames_per_outfile = 12, 12, 12, 12,
restart = .false.,
restart_interval = 5000,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
nocolons = .false.,
io_form_auxinput4 = 2
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval = 180
auxinput4_end_h = 9999
iofields_filename = "iofield_list.txt"
/

&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 472, 166, 201
e_sn = 312, 184, 201
e_vert = 35, 35, 35
num_metgrid_levels = 40,
p_top_requested = 5000,
eta_levels = 1.000, 0.995, 0.990, 0.985, 0.980,
0.970, 0.960, 0.950, 0.940, 0.930,
0.920, 0.910, 0.900, 0.880, 0.860,
0.840, 0.820, 0.800, 0.770, 0.740,
0.700, 0.650, 0.600, 0.550, 0.500,
0.450, 0.400, 0.350, 0.300, 0.250,
0.200, 0.150, 0.100, 0.050, 0.000
dx = 12000, 4000, 1000
dy = 12000, 4000, 1000
grid_id = 1, 2, 3
parent_id = 0, 1, 2
i_parent_start = 0,202, 75
j_parent_start = 0, 35, 40
parent_grid_ratio = 1, 3, 4
parent_time_step_ratio = 1, 3, 4
feedback = 1
smooth_option = 0
numtiles = 1,
nproc_x = 6,
nproc_y = 4,
use_adaptive_time_step = .true.,
step_to_output_time = .true.,
target_cfl = 0.7, 0.7, 0.5, 1.0,
max_step_increase_pct = 5, 11, 11, 22,
starting_time_step = 108, 36, 12, -1,
max_time_step = 216, 72, 15, -1,
min_time_step = 1, 1, 1, -1,
adaptation_domain = 1,

/

&physics
mp_physics = 10, 10, 10
ra_lw_physics = 4, 4, 4
ra_sw_physics = 4, 4, 4
radt = 20, 20, 20
sf_sfclay_physics = 7, 7, 7
sf_surface_physics = 7, 7, 7
bl_pbl_physics = 7, 7, 7
bldt = 0, 0, 0
cu_diag = 1,
cu_physics = 5,
cudt = 0, 0, 0
convtrans_avglen_m = 60
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 2,
sst_update = 1,
sst_skin = 0,
pxlsm_smois_init = 0,0,0
mp_zero_out = 0,
slope_rad = 1,1,1
topo_shading = 1,1,1
shadlen = 25000.,
num_land_cat = 40
prec_acc_dt = 60
/

&fdda
grid_fdda = 1
grid_sfdda = 1
pxlsm_soil_nudge = 1
sgfdda_inname = "wrfsfdda_d<domain>"
sgfdda_end_h = 144
sgfdda_interval_m = 180
sgfdda_interval = 10800
gfdda_end_h = 144
gfdda_interval_m = 180
gfdda_inname = "wrffdda_d<domain>"
fgdt = 0
if_no_pbl_nudging_uv = 1
if_no_pbl_nudging_t = 1
if_no_pbl_nudging_q = 1
if_zfac_uv = 0
k_zfac_uv = 13
if_zfac_t = 0
k_zfac_t = 13
if_zfac_q = 0
k_zfac_q = 13
guv = 0.0003
gt = 0.0003
gq = 0.0000
guv_sfc = 0.0003
gt_sfc = 0.0000
gq_sfc = 0.0000
if_ramping = 1,
dtramp_min = 60.0
rinblw = 250.0
obs_nudge_opt = 0
/

&dynamics
w_damping = 1,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 2,2,2
diff_6th_factor = 0.12,0.12,0.24
base_temp = 290.
damp_opt = 3,
zdamp = 5000.
dampcoef = 0.01
khdif = 0
kvdif = 0
non_hydrostatic = .true.,.true.,.true.
moist_adv_opt = 2, 2, 2
tke_adv_opt = 2, 2, 2
scalar_adv_opt = 2, 2, 2
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false., .false., .false.
nested = .false., .true., .true., .true.
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20120618/c19d0e78/attachment-0001.html