problem with convert_emiss.exe using wrf chem

Introduction to the forum

problem with convert_emiss.exe using wrf chem

Postby Christian813 » Thu Nov 17, 2016 7:12 pm

HI everyone, my name is Christian and I'm facing a problem when I use ./convert_emiss.exe.

I'm trying to simulate a dust storm so I need to use the wrf-chem. After using the DomainWizard to define my domains, then I modify the namelist.input on /run adding the chem section, and also modify the namelist.input on /test/em_real. When I tried convert_emiss.exe, appears this:

"WARNING: dust option 1 currently works only with the GOCART aerosol option."
and also: " READ GOCART BACKGROUND DATA
At line 383 of file convert_emiss.f90 (unit = 19, file = 'wrf_gocart_backg')
Fortran runtime error: End of file"

In the first exercise of the wrf-chem tutorial web page, in a step says "verify that geogrid.log says :
Processing field 25 of 27 (EROD)
Processing field 26 of 27 (CLAYFRAC)
Processing field 27 of 27 (SANDFRAC)

I read the geogrid.log and these lines are missed, so I don't know which is the source of the problem.
This is the final part of the geogrid.log

2016-11-16 09:46:33.942 --- Processing field 1 of 25 (LANDUSEF)
2016-11-16 09:46:34.235 --- Processing field 2 of 25 (LU_INDEX)
2016-11-16 09:46:34.235 --- Processing field 3 of 25 (HGT_M)
2016-11-16 09:46:34.245 --- Processing field 4 of 25 (SOILTEMP)
2016-11-16 09:46:34.248 --- Processing field 5 of 25 (SOILCTOP)
2016-11-16 09:46:34.540 --- Processing field 6 of 25 (SCT_DOM)
2016-11-16 09:46:34.540 --- Processing field 7 of 25 (SOILCBOT)
2016-11-16 09:46:34.832 --- Processing field 8 of 25 (SCB_DOM)
2016-11-16 09:46:34.833 --- Processing field 9 of 25 (ALBEDO12M)
2016-11-16 09:46:34.905 --- Processing field 10 of 25 (GREENFRAC)
2016-11-16 09:46:34.979 --- Processing field 11 of 25 (LAI12M)
2016-11-16 09:46:34.989 --- Processing field 12 of 25 (SNOALB)
2016-11-16 09:46:34.991 --- Processing field 13 of 25 (SLOPECAT)
2016-11-16 09:46:35.000 --- Processing field 14 of 25 (CON)
2016-11-16 09:46:35.012 --- Processing field 15 of 25 (VAR)
2016-11-16 09:46:35.024 --- Processing field 16 of 25 (OA1)
2016-11-16 09:46:35.038 --- Processing field 17 of 25 (OA2)
2016-11-16 09:46:35.052 --- Processing field 18 of 25 (OA3)
2016-11-16 09:46:35.067 --- Processing field 19 of 25 (OA4)
2016-11-16 09:46:35.082 --- Processing field 20 of 25 (OL1)
2016-11-16 09:46:35.094 --- Processing field 21 of 25 (OL2)
2016-11-16 09:46:35.105 --- Processing field 22 of 25 (OL3)
2016-11-16 09:46:35.117 --- Processing field 23 of 25 (OL4)
2016-11-16 09:46:35.129 --- Processing field 24 of 25 (VAR_SSO)
2016-11-16 09:46:35.133 --- Processing field 25 of 25 (LAKE_DEPTH)
2016-11-16 09:46:35.475 --- *** Successful completion of program geogrid.exe ***

Then if I go to geogrid.TBL I find that the EROD, CLAYFRAC and SANDFRAC variables are not contempled, so maybe they need to be added or something, but I don't have a clear idea of how.

My namelist is the next:

&time_control
run_days = 6,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2007, 2007, 2007,
start_month = 10, 10, 10,
start_day = 21, 21, 21,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2007, 2007, 2007,
end_month = 10, 10, 10,
end_day = 23, 23, 23,
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 = 120, 30, 10,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
io_form_auxinput2 = 0,
io_form_auxinput4 = 0,
io_form_auxinput5 = 0,
io_form_auxinput6 = 0,
io_form_auxinput7 = 0,
io_form_auxinput8 = 0,
io_form_auxinput12 = 0,
io_form_auxinput13 = 0,
/

&domains
time_step = 48,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 100, 88, 76,
e_sn = 100, 88, 76,
e_vert = 29, 29, 29,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 9000, 3000, 1000,
dy = 9000, 3000, 1000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 36, 32,
j_parent_start = 1, 36, 32,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0,
/

&physics
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 30, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,
/

&fdda
/

&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 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 = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
chem_adv_opt = 2, 0, 0,
tke_adv_opt = 2, 0, 0,
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 = 0,
chem_opt = 401,
bioemdt = 0,
photdt = 0,
chemdt = 5,
io_style_emissions = 0,
emiss_opt = 0, 0,
emiss_opt_vol = 0, 0,
emiss_ash_hgt = 0.,
chem_in_opt = 0, 0,
phot_opt = 0, 0,
gas_drydep_opt = 0, 0,
aer_drydep_opt = 0, 0,
bio_emiss_opt = 0, 0,
ne_area = 0,
dust_opt = 1,
dmsemis_opt = 0,
seas_opt = 0,
depo_fact = 0,
gas_bc_opt = 0, 0,
gas_ic_opt = 0, 0,
aer_bc_opt = 0, 0,
aer_ic_opt = 0, 0,
gaschem_onoff = 0, 0,
aerchem_onoff = 0, 0,
wetscav_onoff = 0, 0,
cldchem_onoff = 0, 0,
vertmix_onoff = 0, 0,
chem_conv_tr = 0, 0,
conv_tr_wetscav = 0, 0,
conv_tr_aqchem = 0, 0,
biomass_burn_opt = 0, 0,
plumerisefire_frq = 0, 0,
have_bcs_chem = .false., .false., .false.,
aer_ra_feedback = 0,
aer_op_opt = 0,
opt_pars_out = 0,
diagnostic_chem = 0,
/


&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

I hope someone could help me with this, I have a long time stuck on this. Thank you very much!!!
Christian813
 
Posts: 11
Joined: Thu Nov 17, 2016 6:54 pm

Re: problem with convert_emiss.exe using wrf chem

Postby Christian813 » Fri Nov 18, 2016 4:22 pm

Hi everyone, I solved, I just needed to add the lines:
===============================
name = EROD
priority = 1
dest_type = continuous
masked=water
z_dim_name=dust_erosion_dimension
fill_missing=0.
interp_option = default:average_4pt
rel_path = default:erod/
===============================
name = CLAYFRAC
priority = 1
dest_type = continuous
masked=water
fill_missing=0.
interp_option = default:four_pt+average_4pt+average_16pt+search
rel_path = default:clayfrac_5m/
===============================
name = SANDFRAC
priority = 1
dest_type = continuous
masked=water
fill_missing=0.
interp_option = default:four_pt+average_4pt+average_16pt+search
rel_path = default:sandfrac_5m/

in the GEOGRID.TBL file to make it work with dust_opt= 1
Christian813
 
Posts: 11
Joined: Thu Nov 17, 2016 6:54 pm


Return to Welcome to the WRF User's Forum

Who is online

Users browsing this forum: No registered users and 1 guest