Hourly wrffirechemi files

Any issues with running WRF Chemistry

Hourly wrffirechemi files

Postby Eowyn » Fri Aug 10, 2012 1:53 pm

I have hourly fire emissions of the form:
wrffirechemi_d01_2008-08-16_11:00:00
wrffirechemi_d01_2008-08-16_12:00:00
wrffirechemi_d01_2008-08-16_13:00:00
etc

I'm running with biomass burning and plumerise ON. In the &time_control namelist, I have the following:
auxinput7_inname = 'wrffirechemi_d01',
io_form_auxinput7 = 2,
auxinput7_interval = 3600,

WRF/Chem fails to find the hourly emissions data (it is present via symbolic link in the run dir). It runs just fine if one of the hourly files is renamed wrffirechemi_d01. Presumably it considers the data static. I believe the namelist options above are being ignored, for example it runs without complaint if I say:

auxinput7_inname = 'bozo_the_cloud',

so long as there is a wrffirechemi_d01 file available.



The error when trying to use hourly emissions is:

-------------- FATAL CALLED ---------------
med_read_wrf_chem_fireemissions: error opening wrffirechemi_d01
-------------------------------------------

And the debug file (debug level set to 500) ends with:

READING EMISSIONS DATA OPT 3
d01 2008-08-16_00:00:00 DEBUG wrf_timetoa(): returning with str = [2008-08-16_00:00:00]
mediation_integrate: med_read_wrf_chem_fireemissions: Open file wrffirechemi_d01
d01 2008-08-16_00:00:00 module_io.F: in wrf_open_for_read
d01 2008-08-16_00:00:00 NetCDF error: No such file or directory
d01 2008-08-16_00:00:00 NetCDF error in wrf_io.F90, line 1025
FATAL CALLED FROM FILE: mediation_integrate.b LINE: 2647


Thank you in advance for your help.

Regards,
Eowyn Baughman, Univ of Washington
Eowyn
 
Posts: 7
Joined: Fri Aug 10, 2012 1:39 pm

Re: Hourly wrffirechemi files

Postby Eowyn » Mon Aug 13, 2012 7:11 pm

Resolved! Yay!

&time_control
frames_per_auxinput7 = 1,
auxinput7_inname = 'wrffirechemi_d<domain>_<date>',
auxinput7_interval_m = 60,
io_form_auxinput7 = 2,

&chem
biomass_burn_opt = 2 or 3 or 4

biomass_burn_opt CANNOT be 1.

This is all for wrffirechemi files produced by the fire_emis tool from NCAR ACD.
Eowyn
 
Posts: 7
Joined: Fri Aug 10, 2012 1:39 pm

Re: Hourly wrffirechemi files

Postby Benjamin » Fri Aug 17, 2012 5:32 am

Hi Eowyn,

I have the same problem as you had, but unfortunately what you propose didn't solve the issue for me..

here is basically what i have:
&time_control
run_days = 2,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2010, 2010, 2000,
start_month = 01, 07, 01,
start_day = 01, 14, 24,
start_hour = 00, 00, 12,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2010, 2010, 2000,
end_month = 01, 07, 01,
end_day = 03, 16, 25,
end_hour = 00, 00, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true.,.true.,.false.,
history_interval = 60, 60, 60,
frames_per_outfile = 12, 1000, 1000,
restart = .false.,
frames_per_auxinput7 = 1,
restart_interval = 0,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
auxinput6_inname = 'wrfbiochemi_d01',
auxinput7_inname = 'wrffirechemi_d<domain>_<date>',
auxinput8_inname = 'wrfchemi_gocart_bg_d<domain>',
auxinput12_inname = 'wrf_chem_input',
auxinput13_inname = 'wrfchemv_d<domain>',
auxinput5_interval_m = 86400, 86400, 60,
auxinput7_interval_m = 1440,
auxinput8_interval_m = 86400, 86400, 60,
auxinput13_interval_m = 86400, 86400, 60,
io_form_auxinput2 = 2,
io_form_auxinput5 = 2,
io_form_auxinput6 = 0,
io_form_auxinput7 = 2,
io_form_auxinput8 = 2,
io_form_auxinput12 = 0,
io_form_auxinput13 = 0,
frames_per_auxinput7 = 1,
debug_level = 0,
auxinput1_inname = "met_em.d<domain>.<date>",
/


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

Do you see anything different that in your name;list (except i want an update every day)
Thanks in advance,
Benjamin
Benjamin
 
Posts: 6
Joined: Tue Jul 17, 2012 10:00 am

Re: Hourly wrffirechemi files

Postby Eowyn » Fri Aug 17, 2012 6:04 pm

Hi Benjamin,

Is the error occurring at real.exe or wrf.exe?

-Eowyn
Eowyn
 
Posts: 7
Joined: Fri Aug 10, 2012 1:39 pm

Re: Hourly wrffirechemi files

Postby aghedo » Thu Feb 07, 2013 2:03 pm

Benjamin,
I know this is quite a while ago...have you been able to resolve the problem? I have found that inconsistencies in the namelist controlling WRF-Chem emissions could cause errors. Sometimes, the problem may have occurred before the error messages written in between the "FATAL CALLED". This is such that even though WRF-Chem point to a particular emissions file, the problem may have occurred prior to reading this file, and may be with another file for another emissions category.

Now to the specific of your namelist, I see some possible inconsistencies already:

    I don't know how you set up your emissions, but your "auxinput*_interval_m" look rather suspicious. For example, combining "auxinput7_interval_m = 1440", "frames_per_auxinput7 = 1", and "auxinput7_inname = 'wrffirechemi_d<domain>_<date>'" will mean you have at least 2 files, each containing a daily mean fire emissions. Unfortunately, the problem with this setup is that wrf-chem will try to look for fire emissions at every <date>, which in this case will be "every day"..is this what you intended?

Also,check the following for consistency:

    io_style_emissions = 2, (this means you have anthropogenic emissions for every hour of your simulations. This will mean for your 2 days run, you must have 48 files. If this is not the case, then WRF-Chem will complain).

    set "auxinput5_inname" for your anthropogenic emissions, right now, the model will look for files with names: wrfchemi_d<domain>_<date>. Again because your io_style_emissions=2
    I suggest you define "frames_per_auxinput5" (if they are indeed available for all hours of your simulation), and "frames_per_auxinput8", "frames_per_auxinput13"
    I will remove:
    io_form_auxinput2 = 2,
    io_form_auxinput5 = 2,
    Are you sure about these?
    io_form_auxinput6 = 0
    io_form_auxinput12 = 0,
    io_form_auxinput12 = 0,

Hope this helps

~Tutu Aghedo
aghedo
 
Posts: 2
Joined: Wed Feb 06, 2013 11:17 am

Re: Hourly wrffirechemi files

Postby FC_stl » Wed Dec 09, 2020 1:55 am

Also I would recommend checking your settings on your "auxinput ?_interval_m". The "m" in the end means the numbers you type in are with the unit of "minute". 86400 minutes is 60 days, which seems not right.
FC_stl
 
Posts: 2
Joined: Sat Sep 21, 2019 1:26 am


Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 1 guest