Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Why is the "qndrop" not a constant in the case of progn=0?

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

cxyiyh

Member
Hi, friends,
When I set mp_physics = 10(Morrison-2-moment) and progn = 0,the cloud droplet should be constant,however, I change the Registry.COMMON to output in the case of progn=0, and the variable "qndrop" is not a constant.
Here is a picture shows regional average qndrop changes over time.
cloud-water.png
In my opnion, the variable "qndrop" represents cloud droplet number concentration.

Best Regards!
yours,
cxyiyh
 
Ming Chen said:
Please send me your namelist.input to take a look. thanks.
Thank you for your reply!Here is my namelist.input
&time_control
run_days = 3,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2014, 2014, 2014,
start_month = 04, 04, 04,
start_day = 22, 22, 22,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2014, 2014, 2014,
end_month = 04, 04, 04,
end_day = 25, 25, 25,
end_hour = 00, 00, 00,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,
history_interval = 360, 180, 90
frames_per_outfile = 1, 1, 1,
restart = .false.,
restart_interval = 100000,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,

io_form_auxinput5 = 2,
io_form_auxinput7 = 2,
io_form_auxinput8 = 2,

auxinput5_inname = 'wrfchemi_d<domain>_<date>',
auxinput7_inname = 'wrffirechemi_d<domain>_<date>',
auxinput8_inname = 'wrfchemi_gocart_bg_d<domain>_<date>',
auxinput5_interval_m = 14400,14400,14400,
auxinput7_interval_m = 14400,14400,14400,
auxinput8_interval_m = 14400,14400,14400,

frames_per_auxinput5 = 1, 1, 1,
frames_per_auxinput7 = 1, 1, 1,
frames_per_auxinput8 = 1, 1, 1,
/

&dfi_control
/

&domains
time_step = 90,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
s_we = 1, 1, 1,
e_we = 230,181,181,
s_sn = 1, 1, 1,
e_sn = 180,151,151,
e_vert = 34, 34, 34,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 27000, 9000, 3000,
dy = 27000, 9000, 3000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 83, 55,
j_parent_start = 1, 65, 54,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
p_top_requested = 5000,
feedback = 1,
smooth_option = 0
zap_close_levels = 50
interp_type = 1
t_extrap_type = 2
force_sfc_in_vinterp = 0
use_levels_below_ground = .true.
use_surface = .true.
lagrange_order = 1
/

&physics
mp_physics = 10, 10, 10,
progn = 0, 0, 0,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
radt = 30, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 3, 3, 3,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 3, 0, 0,
cu_diag = 1, 0, 0,
cudt = 5, 5, 5,
ishallow = 0,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-12
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,
num_land_cat = 21,
/

&fdda
/

&dynamics
rk_ord = 3,
w_damping = 0,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 2, 2, 2,
scalar_adv_opt = 2, 2, 2,
chem_adv_opt = 2, 2, 2,
tke_adv_opt = 2, 2, 2,
time_step_sound = 4, 4, 4,
h_mom_adv_order = 5, 5, 5,
v_mom_adv_order = 3, 3, 3,
h_sca_adv_order = 5, 5, 5,
v_sca_adv_order = 3, 3, 3,
/

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

&grib2
/

&chem
kemit = 1,
chem_conv_tr = 1, 0, 0,
chem_opt = 8, 8, 8,
bioemdt = 30, 30, 30,
photdt = 30, 30, 30,
chemdt = 1.5, 1.5, 1.5,
io_style_emissions = 2,
emiss_opt = 5, 5, 5,
emiss_inpt_opt = 1, 1, 1,
emiss_opt_vol = 0, 0, 0,
emiss_ash_hgt = 20000.,
chem_in_opt = 0, 0, 0,
phot_opt = 1, 1, 1,
gas_drydep_opt = 1, 1, 1,
aer_drydep_opt = 301, 301, 301,
bio_emiss_opt = 0, 0, 0,
ne_area = 0,
dust_opt = 3,
dmsemis_opt = 1,
seas_opt = 1,
depo_fact = 0.25,
gas_bc_opt = 1, 1, 1,
gas_ic_opt = 1, 1, 1,
aer_bc_opt = 1, 1, 1,
aer_ic_opt = 1, 1, 1,
gaschem_onoff = 1, 1, 1,
aerchem_onoff = 1, 1, 1,
wetscav_onoff = 0, 0, 0,
cldchem_onoff = 0, 0, 0,
vertmix_onoff = 1, 1, 1,
conv_tr_wetscav = 1, 0, 0,
conv_tr_aqchem = 0, 0, 0,
biomass_burn_opt = 1, 1, 1,
plumerisefire_frq = 30, 30, 30,
have_bcs_chem = .false., .false., .false.,
aer_ra_feedback = 1, 1, 1,
aer_op_opt = 1, 1, 1,
opt_pars_out = 1,
diagnostic_chem = 0, 0, 0,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
 
I looked at your namelist.input, which indicates that you are running wrf-chem. Please let me know if I am wrong.
Looking at the code for wrf-chem, qndrop is calculated and updated in module_mp_morr_two_moment.F. Below is a piece of code that updates qndrop:

! wrf-chem
IF (flag_qndrop .AND. PRESENT( qndrop )) THEN
qndrop(i,k,j) = nc1d(k)
!jdf CSED3D(I,K,J) = CSED(K)
END IF

Note that some physics schemes may work differently between WRF-ARW and WRF-Chem.
 
Top