WRF-CHEM，WRF-VPRM

zicheng zheng

unread,

May 28, 2024, 10:05:01 PM5/28/24

to wrf-chem

Dear Sir,

I hope this email finds you well.

I would like to thank you for reading my letter and I hope to get some help from you.

My base case is that I want to analyze the variation of CO2 using WRF4.3 version with chem_opt=16, and for this I want to use the VPRM with auxinput=15 to achieve the natural carbon flux.

At present, I have completed the operation of VPRMpreproc, obtained 7 variables related to LSWI and EVI, and merged them into the same nc. For this reason, I also read the WRF_GHG related materials. I set up namelist.input according to the data I have at hand and tried to run it. To do this, I prepared the relevant list of anthropogenic emissions, and use auxinput = 5, CO2_ANT in wrfout has change obviously, but CO2_BIO almost no change, and no other variables appeared in wrfout except CO2_BIO.

I tried to set the input files of VPRM to 0 and also got the above results, so I guess there is a problem with my VPRM input files.When I set my debug_level to 300, some prompts appear, like

NetCDF error: NetCDF: Attribute not found,

input_wrf.F reading 4d real VEGFRA_VPRM Status = -1021.

Now I don't know how to proceed to the next step, such as how to correctly merge the 7 variables of VPRM into the same one, and how to set the related parameters of namelist.input. I can't do anything about it.

Please help me. If possible, please provide me with a specific example and file of VPRM input file and namelist Settings.

Thanks!

Best wishes，

zicheng zheng

unread,

May 28, 2024, 10:06:07 PM5/28/24

to wrf-chem, zicheng zheng

my namelist.input is as follows:
&time_control
run_days = 0,
run_hours =0,
run_minutes = 0,
run_seconds = 0,
start_year = 2023, 2023
start_month = 01, 01,
start_day =01, 01,
start_hour = 00, 00,
start_minute = 00, 00,
start_second = 00, 00,
end_year = 2023, 2023 ,
end_month = 12, 12,
end_day = 31, 31,
end_hour = 00, 00,
end_minute = 00, 00,
end_second = 00, 00,
interval_seconds = 43200
input_from_file = .true.,.true.,
history_interval = 60, 60,
frames_per_outfile = 24, 24,
restart = .false.,
restart_interval = 1440,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput2 = 2,
io_form_auxinput5 = 2,
io_form_auxinput6 = 2,
io_form_auxinput7 = 2,
io_form_auxinput8 = 2,
io_form_auxinput15 = 2,
!iofields_filename = "myvprmco2_file_d01.txt","myvprmco2_file_d02.txt","myvprmco2_file_d03.txt",
auxinput5_interval_m =1440, 1440,
auxinput6_interval_m = 1440, 1440,
auxinput7_interval_m =1440, 1440,
auxinput8_interval_m =1440, 1440,
auxinput25_interval_m =11520,11520,

! frames_per_auxinput5 = 1,1,1,

debug_level = 300,

auxinput5_inname = 'wrfchemi_d<domain>.nc',
auxinput6_inname = 'wrfbiochemi_d<domain>.nc',
auxinput7_inname = 'wrffirechemi_d<domain>.nc',
auxinput8_inname = 'wrfchemi_gocart_bg_d<domain>',
auxinput15_inname = 'VPRM_input_2023_d<domain>.nc',

force_use_old_data =.true.
/

&domains
time_step = 50,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we =40, 40,
e_sn =40, 40,
e_vert = 34, 34,
p_top_requested = 5000,
num_metgrid_levels = 34,
num_metgrid_soil_levels = 4,
dx = 9000, 3000,
dy = 9000, 3000,
grid_id = 1, 2,
parent_id = 0, 1,
i_parent_start = 1, 15,
j_parent_start = 1, 5,
parent_grid_ratio = 1, 3,
parent_time_step_ratio = 1, 3,
feedback = 1,
smooth_option = 0,
sfcp_to_sfcp = .true.
/

&physics
mp_physics = 2, 2,
ra_lw_physics = 1, 1,
ra_sw_physics = 1, 1,
radt = 27, 9,
sf_sfclay_physics = 1, 1,
sf_surface_physics = 1, 1,
bl_pbl_physics = 1, 1,
bldt = 0, 0,
cu_physics = 1, 1,
cudt = 5, 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_land_cat = 21,
num_soil_layers = 2,
sf_urban_physics = 0, 0,
mfshconv = 0,
sst_update = 1,
/

&fdda

/

&dynamics
rk_ord = 3, !积分方案 3阶
w_damping = 1, !阻尼
diff_opt = 1, 1, 1, 1,
km_opt = 4, 4, 4, 4,
diff_6th_opt = 0, 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2, 0.2
khdif = 0, 0, 0, 0,
kvdif = 0, 0, 0, 0,
non_hydrostatic = .true., .true., .true., .true.,
moist_adv_opt = 1, 1, 1, 1,
scalar_adv_opt = 1, 1, 1, 1,
/

&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,
/

&chem
kemit = 6,
chem_opt = 16,16,16,
bioemdt = 60,60,60,
photdt = 0,0,0,
chemdt = 10,3.3333,1,
io_style_emissions = 2,
emiss_inpt_opt = 16,16,16,
emiss_opt = 16,16,16,
emiss_opt_vol = 0,0,0,
emiss_ash_hgt = 20000.,
chem_in_opt = 0,0,0,
phot_opt = 0,0,0,
gas_drydep_opt = 0,0,0,
aer_drydep_opt = 1,1,1,
bio_emiss_opt = 16,
ne_area = 100,
dust_opt = 0,
dmsemis_opt = 0,
seas_opt = 0,
depo_fact = 0.25,0.25,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,
vprm_opt = 'VPRM_table_US','VPRM_table_US','VPRM_table_US',
!term_opt = 'CH4_termite_OW'
!wpeat = 0.05
!wflood = 0.19
vertmix_onoff = 1,1,1,
chem_conv_tr = 0,0,0,
conv_tr_wetscav = 0,0,0,
conv_tr_aqchem = 0,0,0,
biomass_burn_opt = 0,0,0,
plumerisefire_frq = 120,120,120,
have_bcs_chem = .true.,.true.,.true.,
aer_ra_feedback = 0,0,0,
aer_op_opt = 0,
opt_pars_out = 0,
diagnostic_chem = 0,0,0,
/

zicheng zheng

unread,

May 28, 2024, 10:07:31 PM5/28/24

to wrf-chem, zicheng zheng

More information can be found in my attachment.

rls.out.jpg

namelist.input

VPRM_input_2023_d02.nc

Gabriele Pfister

unread,

May 29, 2024, 1:08:15 PM5/29/24

to zhengzic...@igsnrr.ac.cn, wrf-chem

Hi,

we are unfortunately not able to help with your question since our NCAR User Forum is focused on the by us implemented MOZART chemical schemes and WRF-chem preprocessing tools we provide. Please direct your question to the general WRF/MPAS User Forum: https://forum.mmm.ucar.edu/

Gabi

--

================================
Gabriele Pfister

Atmospheric Chemistry Observations & Modeling (ACOM) Lab

NSF National Center for Atmospheric Research (NCAR)

Email: pfi...@ucar.edu

Work Phone: +1 303 497 2915

Web: https://staff.ucar.edu/users/pfister

zicheng zheng

unread,

May 30, 2024, 10:24:49 AM5/30/24

to wrf-chem, pfi...@ucar.edu, wrf-chem, zicheng zheng

Thank you for your reply as well. I'm sorry I asked an inappropriate question
In fact, I came here from the forum, because I could not wait for a reply in the paper, probably due to his activity.
thank you.

Reply all

Reply to author

Forward