Hi all,I am trying to run wrf-chem with chem_opt= 301, I have created anthro, biogenic, and fire emissions using anthro_emiss, bio_emiss and fire_emiss utilities.
While running real.exe, i am getting the error:
(base) [wcrg23@login03 em_real]$ tail rsl.error.0000
d01 2022-11-01_00:00:00 DEBUG wrf_timetoa(): returning with str = [2022-11-01_00:00:00]
mediation_integrate: med_read_wrf_chem_fireemissions: Open file wrffirechemi_d01
d01 2022-11-01_00:00:00 module_io.F: in wrf_open_for_read
d01 2022-11-01_00:00:00 NetCDF error: No such file or directory
d01 2022-11-01_00:00:00 NetCDF error in wrf_io.F90, line 1077
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: mediation_integrate.b LINE: 2885
med_read_wrf_chem_fireemissions: error opening wrffirechemi_d01
-------------------------------------------
Abort(1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
Seems like real program is looking for wrffirechemi_d01 instead of wrffirechemi_d<domain>_<date> even though
auxinput7_inname = 'wrffirechemi_d<domain>_<date>',is set in the namelist.
I have tried running the real by linking the first hour wrffirechemi_d01_2022-11-01_00:00:00 to wrffirechemi_d01, which strangely works even withauxinput7_interval_m = 60,frames_per_auxinput7 = 1,
I could not figure out whether the problem is related to the chemistry settings in the namelist.input file or the emissions input files.
Please guide if you have expertise on this.Thank you
Hi Gabi,
Thank you so much for your feedback and observations!
I'm currently following the instructions from WRF-Chem Exercise 6 along with a YouTube video on WRF-Chem/mozbc:
Regarding real.exe, you're correct that I initially ran it just for meteorology to generate the wrfinput_d01 file. The subsequent runs were aimed at incorporating chemistry data and updating it for certain gaseous species and particulate matter using the mozbc utility.
Earlier, I used chem_opt=112 for the same emission dataset, where the second run of real.exe was successful, and the WRF output seemed fine.
During troubleshooting, I turned off some variables one at a time to pinpoint the issue and accidentally left phot_opt=0. I’ll be using phot_opt=1 for the final simulations, as photolysis is definitely necessary.
Interestingly, the WRF output seems to run fine without updating wrfinput and boundary files for chemical inputs, so I’m starting to think the issue might lie in the namelist settings.
I would really appreciate your thoughts on this.
Thanks again for your help!
Best regards,
Vikas