Error while running ANTHRO-EMISS

750 views
Skip to first unread message

Yagni Rami

unread,
Apr 28, 2020, 12:34:42 PM4/28/20
to wrf-chem-anthro_emiss
Hello,

I am Yagni, a PhD student working on Air Quality modelling using WRF- Chem ..

I have successfully run real.exe and wrf.exe (without Chemistry)..

Now I am trying to run Anthro- Emissions.. Generated log file is attached..
Error:  anthro_emis: anthro source file
                path_to/EDGAR-HTAP/MOZART/EDGAR_HTAP_emi_CO_2010.0.1x0.1.nc
 anthro_emis: does not exist

This is what I am doing..
1)  cbmz-mozaic.inp is input file
2) ./run.anthro.csh file (./anthro_emis < cbmz-mosaic.inp > abc.log)

which gives above error..  Error displayed on terminal is STOP NETCDF err or STOP FILE err... I have tried multiple times to run and it display either of this..Preview attachment Header_EDGAR_HTAP_emi_CO_2010.0.1x0.1.nc.txt
Header_EDGAR_HTAP_emi_CO_2010.0.1x0.1.nc.txt
abc.log
cbmz-mosaic.inp (1).rami

Stacy Walters

unread,
Apr 28, 2020, 2:12:55 PM4/28/20
to Yagni Rami, wrf-chem-anthro_emiss
Yagni,

Please change the namelist variable anthro_dir from:

'/vols/fs1/abakus/rami/anthro_emis/EDGAR-HTAP/MOZART'

to:

'/vols/fs1/abakus/rami/anthro_emis/EDGAR-HTAP/MOZART/'

(Add the forward slash character, / , to the end of the anthro_dir namelist variable.

Then try again.  If that fails we'll dive a little deeper.

Stacy

Ankan Sarkar

unread,
Mar 14, 2023, 3:48:58 AM3/14/23
to wrf-chem-anthro_emiss, Stacy Walters, wrf-chem-anthro_emiss, Yagni Rami
Dear all,
I want to perform WRF-Chem simulations using the CBMZ-MOSAIC chemistry option (e.g., chem_opt = 32) to study aerosol-cloud interaction. I also want to use anthropogenic emissions from EDGAR-HTAP (a global monthly emission inventory for 2010) or  EDGARv5_MOZART_dataset (from the EDGAR v5.0 emissions inventory speciated for the MOZART chemical mechanism). I noticed that only the sample input files for MOZCART and MOZART-MOSAIC chemistry options come with utilities from NCAR. 

But I want to use the CBMZ-MOSAIC chemistry option, as Yagni did. However, the 'src_names' variables set in the input files provided by Yagni are slightly different from those in my EDGAR-HTAP dataset. For example, in my case, NOx(30), while in Yagni's case, it is NOx (31.6). And, also in my case, there are no 'CH3COCH3MEK(65)' ,' PM2.5other(1)', and 'PMcoarse(1) like in Yagni's case. That is why the 'emis_map' variable should be somewhat different in my case. I searched online a lot, but I couldn't find any guides on setting the 'emis_map' variable in the input file for running anthro_emis.
Can anyone please guide me on setting the 'emis_map' variable according to the 'src_names' of my EDGAR-HTAP dataset? 
 For your convenience, I am giving the 'src_names' in the input file for the EDGAR-HTAP dataset as follows:
src_names ='CO(28)','NOx(30)','SO2(64)','NH3(17)','BC(12)','OC(12)','PM2.5(1)','PM10(1)',
             'BIGALK(72)','BIGENE(56)','C2H4(28)','C2H5OH(46)','C2H6(30)','CH2O(30)',
             'CH3CHO(44)','CH3COCH3(58)','CH3OH(32)','MEK(72)','TOLUENE(92)','C3H6(42)','C3H8(44)',
             'BENZENE(78)','XYLENE(106)'
And for the EDGARv5_MOZART_dataset 'src_names' in the input file is as follows:
 src_names = 'CO(28)','NOx(30)','SO2(64)','NH3(17)','BC(12)','OC(12)','PM2.5(1)','PM10(1)','BENZENE(78)','BIGALK(72)','BIGENE(56)','C2H2(26)','C2H4(28)',
             'C2H5OH(46)','C2H6(30)','C3H6(42)','C3H8(44)','CH2O(30)',
             'CH3CHO(44)','CH3COCH3(58)','CH3OH(32)','CH3COOH(60)','HCOOH(46)',
             'MEK(72)','TOLUENE(92)','XYLENES(106)'

Any help on this will be greatly appreciated. Thank you for your time and consideration.
With regards,
Ankan


Ankan Sarkar

unread,
Mar 17, 2023, 10:34:14 AM3/17/23
to wrf-chem-anthro_emiss, Ankan Sarkar, Stacy Walters, wrf-chem-anthro_emiss, Yagni Rami, Gabriele Pfister
Dear all,
I checked the Registry/registry.chem, and I found the following lines related to anthropogenic emissions and CBMZ-MOSAIC:
# Anthropogenic emissions
#
state   real     -            i+jf     emis_ant     -         -      -       -                     "Anthropogenic Emissions"         ""
state    real  e_iso          i+jf     emis_ant     1         Z      i5r    "E_ISO"               "Isoprene EMISSIONS (Anth. for RADM/RACM, Anth+Bio for CBMZ)"  "mol km^-2 hr^-1"
state    real  e_so2          i+jf     emis_ant     1         Z      i5r    "E_SO2"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_no           i+jf     emis_ant     1         Z      i5r    "E_NO"                "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_co           i+jf     emis_ant     1         Z      i5r    "E_CO"                "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_co2          i+jf     emis_ant     1         Z      i5r    "E_CO2"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_eth          i+jf     emis_ant     1         Z      i5r    "E_ETH"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_hc3          i+jf     emis_ant     1         Z      i5r    "E_HC3"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_hc5          i+jf     emis_ant     1         Z      i5r    "E_HC5"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_hc8          i+jf     emis_ant     1         Z      i5r    "E_HC8"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_xyl          i+jf     emis_ant     1         Z      i5r    "E_XYL"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_ol2          i+jf     emis_ant     1         Z      i5r    "E_OL2"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_olt          i+jf     emis_ant     1         Z      i5r    "E_OLT"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_oli          i+jf     emis_ant     1         Z      i5r    "E_OLI"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_tol          i+jf     emis_ant     1         Z      i5r    "E_TOL"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_csl          i+jf     emis_ant     1         Z      i5r    "E_CSL"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_hcho         i+jf     emis_ant     1         Z      i5r    "E_HCHO"              "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_ald          i+jf     emis_ant     1         Z      i5r    "E_ALD"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_ket          i+jf     emis_ant     1         Z      i5r    "E_KET"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_ora2         i+jf     emis_ant     1         Z      i5r    "E_ORA2"              "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_nh3          i+jf     emis_ant     1         Z      i5r    "E_NH3"               "EMISSIONS"          "mol km^-2 hr^-1"
state    real  e_pm_25        i+jf     emis_ant     1         Z      i5r    "E_PM_25"             "EMISSIONS"          "ug/m3 m/s"
state    real  e_pm_10        i+jf     emis_ant     1         Z      i5r    "E_PM_10"             "EMISSIONS"          "ug/m3 m/s"
state    real  e_pm25i        i+jf     emis_ant     1         Z      i5r    "E_PM25I"             "EMISSION RATE OF UNIDEN. PM2.5 MASS"  "ug/m3 m/s"
state    real  e_pm25j        i+jf     emis_ant     1         Z      i5r    "E_PM25J"             "EMISSION RATE OF J-MODE UNIDEN. PM2.5 MASS"   "ug/m3 m/s"
state    real  e_eci          i+jf     emis_ant     1         Z      i5r    "E_ECI"               "EMISSION RATE OF I-MODE EC"          "ug/m3 m/s"
state    real  e_ecj          i+jf     emis_ant     1         Z      i5r    "E_ECJ"               "EMISSION RATE OF J-MODE EC"          "ug/m3 m/s"
state    real  e_naai          i+jf     emis_ant     1         Z      i5r   "E_NAAI"               "EMISSION RATE OF I-MODE EC"          "ug/m3 m/s"
state    real  e_naaj          i+jf     emis_ant     1         Z      i5r   "E_NAAJ"               "EMISSION RATE OF J-MODE EC"          "ug/m3 m/s"
state    real  e_orgi         i+jf     emis_ant     1         Z      i5r     "E_ORGI"              "EMISSION RATE OF I-MODE ORG. AER."   "ug/m3 m/s"
state    real  e_orgj         i+jf     emis_ant     1         Z      i5r     "E_ORGJ"              "EMISSION RATE OF J-MODE ORG. AER."   "ug/m3 m/s"
state    real  e_so4i         i+jf     emis_ant     1         Z      i5r     "E_SO4I"              "EMISSION RATE OF I-MODE ORG. AER."   "ug/m3 m/s"
state    real  e_so4j         i+jf     emis_ant     1         Z      i5r     "E_SO4J"              "EMISSION RATE OF J-MODE ORG. AER."   "ug/m3 m/s"
state    real  e_no3i         i+jf     emis_ant     1         Z      i5r     "E_NO3I"              "EMISSION RATE OF I-MODE ORG. AER."   "ug/m3 m/s"
state    real  e_no3j         i+jf     emis_ant     1         Z      i5r     "E_NO3J"              "EMISSION RATE OF J-MODE ORG. AER."   "ug/m3 m/s"
# Additional CBMZ and MOSAIC emission variables...
state    real  e_no2          i+jf     emis_ant     1         Z      i5r     "E_NO2"                "EMISSIONS NO2"                 "mol km^-2 hr^-1"
state    real  e_ch3oh        i+jf     emis_ant     1         Z      i5     "E_CH3OH"              "EMISSIONS CH3OH"               "mol km^-2 hr^-1"
state    real  e_c2h5oh       i+jf     emis_ant     1         Z      i5     "E_C2H5OH"             "EMISSIONS C2H5OH"              "mol km^-2 hr^-1"
state    real  e_so4c         i+jf     emis_ant     1         Z      i5     "E_SO4C"               "EMISSIONS COARSE SO4 AER"      "ug/m3 m/s"
state    real  e_no3c         i+jf     emis_ant     1         Z      i5     "E_NO3C"               "EMISSIONS COARSE NO3 AER"      "ug/m3 m/s"
state    real  e_nh4c         i+jf     emis_ant     1         Z      i5     "E_NH4C"               "EMISSIONS COARSE NH4 AER"      "ug/m3 m/s"
state    real  e_nac          i+jf     emis_ant     1         Z      i5     "E_NAC"                "EMISSIONS COARSE NA AER"      "ug/m3 m/s"
state    real  e_clc          i+jf     emis_ant     1         Z      i5     "E_CLC"                "EMISSIONS COARSE CL AER"      "ug/m3 m/s"
state    real  e_orgc         i+jf     emis_ant     1         Z      i5     "E_ORGC"               "EMISSIONS COARSE ORG AER"      "ug/m3 m/s"
state    real  e_ecc          i+jf     emis_ant     1         Z      i5     "E_ECC"                "EMISSIONS COARSE EC AER"       "ug/m3 m/s"
state    real  e_orgi_a       i+jf     emis_ant     1         Z      i5r     "E_ORGI_A"                "EMISSIONS ORGI_A AER"       "ug/m3 m/s"
state    real  e_orgj_a       i+jf     emis_ant     1         Z      i5r     "E_ORGJ_A"                "EMISSIONS ORGJ_A AER"       "ug/m3 m/s"
state    real  e_orgi_bb      i+jf     emis_ant     1         Z      i5r     "E_ORGI_BB"                "EMISSIONS ORGI_BB AER"       "ug/m3 m/s"
state    real  e_orgj_bb      i+jf     emis_ant     1         Z      i5r     "E_ORGJ_BB"                "EMISSIONS ORGJ_BB AER"       "ug/m3 m/s"
state    real  e_co_a         i+jf     emis_ant     1         Z      i5     "E_CO_A"                "EMISSIONS CO_A AER"       "ug/m3 m/s"
state    real  e_co_bb        i+jf     emis_ant     1         Z      i5     "E_CO_BB"                "EMISSIONS CO_BB AER"       "ug/m3 m/s"
But I couldn't figure out how to prepare 'emis_map'  for running the 'anthro_emis' utility by looking at this while considering the 'src_names' of the EDGARV5-MOZART or EDGAR-HTAP datasets (mentioned in the previous message) for CBMZ-MOSAIC (chem_opt = 31 and 32). Can anyone please provide me the 'emis_map' or input file for CBMZ-MOSAIC? That will be very helpful. Thank you for your time and consideration.
With regards,
Ankan

arsen...@gmail.com

unread,
Sep 9, 2024, 11:04:31 AM9/9/24
to wrf-chem-anthro_emiss, Ankan Sarkar, wrf-chem-anthro_emiss, Yagni Rami
Dear all,
I also want to run a WRF-Chem simulation including aerosol-cloud interaction. However, I am confused about how to generate anthropogenic emissions using EDGAR-HTAP for CBMZ-MOSAIC. The anthro_emiss seems to be developing for MOZART mechanism rather than CBMZ (The .inp file is also not clear about name mapping). I wonder which tool you choose to generate the anthropogenic emissions. If you use anthro_emis, could you share your anthro_emis.inp file? Thanks for your patient help!
With regards,
Wenshuo
Reply all
Reply to author
Forward
0 new messages