(EPA_)anthro_emis: failed to read namelist; error = 5010

[elham baranizadeh]

Hi,

I'm trying to generate wrf--chem-ready emissions from SMOKE (output) emissions (both gridded/sector and inline/stack) using the NEW EPA_anthro_emiss but I'm getting the following error when executing the anthro_emiss:

anthro_emis: failed to read namelist; error =         5010

any idea why I'm getting this error?

Here is the namelist settings (please note I just used the same emis_map as in the README of anthro_emis):

&CONTROL
 anthro_dir = '/XYZ1/EPA_Anthro_Emiss/src/emissData'

 wrf_dir = '/XYZ2/WRF4p3/WRF/run'
 src_lon_dim_name = 'COL'
 src_lat_dim_name = 'ROW'
 domains = 3
 
 cat_var_prefix = ' '
 output_interval = 3600

 sec_file_prefix = 'Merged_Base2018_OGVSMK3D_'
 sec_file_suffix = '.ncf'

 stk_file_prefix = 'inlnts_l_'
 stk_file_suffix = '_4km_Base2018_Point.ncf'
 stk_grp_file_suffix = '_4km_Base2018_Point.ncf''

 src_names(1:6) = 'GID619:epa-sector','aircraft:epa-stack','wRxBurn:epa-stack'
 
 sub_categories(1:10)   = 'CO','NO','NO2','SO2','NH3','ETOH','PAR','IOLE','ETH','ETHA'
 sub_categories(11:20)  = 'OLE','PRPA','FORM','ALD2','ALDX','ACET','MEOH','KET','TOL','BENZ'
 sub_categories(21:25)  = 'XYLMN','ISOP','TERP','SULF','ETHY'
 sub_categories(26:32)  = 'PMOTHR','PEC','POC','PMC','PSO4','PNO3','PNH4'

 start_output_time = '2018-01-06_00:00:00'
 stop_output_time  = '2018-01-06_12:00:00'
 emissions_zdim_stag = 10

 emis_map(1) = 'CO->GID619(CO)+aircraft(CO)+wRxBurn(CO)'
 emis_map(2) = 'NO->GID619(NO)+aircraft(NO)+wRxBurn(NO)'
 emis_map(3) = 'NO2->GID619(NO2)+aircraft(NO2)+wRxBurn(NO2)'
 emis_map(4) = 'SO2->GID619(SO2)+aircraft(SO2)+wRxBurn(SO2)'
 emis_map(5) = 'NH3->GID619(NH3)+aircraft(NH3)+wRxBurn(NH3)'
 emis_map(6) = 'C2H5OH->GID619(ETOH)+aircraft(ETOH)+wRxBurn(ETOH)'
 emis_map(7) = 'BIGALK->GID619(PAR)+aircraft(PAR)+wRxBurn(PAR)'
 emis_map(8) = 'BIGENE->GID619(IOLE)+aircraft(IOLE)+wRxBurn(IOLE)'
 emis_map(9) = 'C2H4->GID619(ETH)+aircraft(ETH)+wRxBurn(ETH)'
 emis_map(10) = 'C2H6->GID619(ETHA)+aircraft(ETHA)+wRxBurn(ETHA)'
 emis_map(11) = 'C3H6->GID619(OLE)+aircraft(OLE)+wRxBurn(OLE)'
 emis_map(12) = 'C3H8->GID619(PRPA)+aircraft(PRPA)+wRxBurn(PRPA)'
 emis_map(13) = 'CH2O->GID619(FORM)+aircraft(FORM)+wRxBurn(FORM)'
 emis_map(14) = 'CH3CHO->GID619(ALD2+ALDX)+aircraft(ALD2+ALDX)+wRxBurn(ALD2+ALDX)'
 emis_map(15) = 'CH3COCH3->GID619(ACET)+aircraft(ACET)+wRxBurn(ACET)'
 emis_map(16) = 'CH3OH->GID619(MEOH)+aircraft(MEOH)+wRxBurn(MEOH)'
 emis_map(17) = 'MEK->GID619(KET)+aircraft(KET)+wRxBurn(KET)'
 emis_map(18) = 'TOLUENE->GID619(TOL)+aircraft(TOL)+wRxBurn(TOL)'
 emis_map(19) = 'BENZENE->GID619(BENZ)+aircraft(BENZ)+wRxBurn(BENZ)'
 emis_map(20) = 'XYLENE->GID619(XYLMN)+aircraft(XYLMN)+wRxBurn(XYLMN)'
 emis_map(21) = 'ISOP->GID619(ISOP)+aircraft(ISOP)+wRxBurn(ISOP)'
 emis_map(22) = 'C10H16->GID619(TERP)+aircraft(TERP)+wRxBurn(TERP)'
 emis_map(23) = 'sulf->GID619(SULF)+aircraft(SULF)+wRxBurn(SULF)'
 emis_map(24) = 'C2H2->GID619(ETHY)+aircraft(ETHY)+wRxBurn(ETHY)'
 emis_map(25) = 'PM_25(A)->GID619(PMOTHR)+aircraft(PMOTHR)+wRxBurn(PMOTHR)'
 emis_map(26) = 'BC(A)->GID619(PEC)+aircraft(PEC)+wRxBurn(PEC)'
 emis_map(27) = 'OC(A)->GID619(POC)+aircraft(POC)+wRxBurn(POC)'
 emis_map(28) = 'PM_10(A)->GID619(PMC)+aircraft(PMC)+wRxBurn(PMC)'
 emis_map(29) = 'SO4I(A)->GID619(PSO4)+aircraft(PSO4)+wRxBurn(PSO4)'
 emis_map(30) = 'SO4J(A)->GID619(PSO4)+aircraft(PSO4)+wRxBurn(PSO4)'
 emis_map(31) = 'ECI(A)->GID619(PEC)+aircraft(PEC)+wRxBurn(PEC)'
 emis_map(32) = 'ECJ(A)->GID619(PEC)+aircraft(PEC)+wRxBurn(PEC)'
 emis_map(33) = 'ORGI(A)->GID619(POC)+aircraft(POC)+wRxBurn(POC)'
 emis_map(34) = 'ORGJ(A)->GID619(POC)+aircraft(POC)+wRxBurn(PO)'
 emis_map(35) = 'NO3I(A)->GID619(PNO3)+aircraft(PNO3)+wRxBurn(PNO3)'
 emis_map(36) = 'NO3J(A)->GID619(PNO3)+aircraft(PNO3)+wRxBurn(PNO3)'
 emis_map(37) = 'NH4I(A)->GID619(PNH4)+aircraft(PNH4)+wRxBurn(PNH4)'
 emis_map(38) = 'NH4J(A)->GID619(PNH4)+aircraft(PNH4)+wRxBurn(PNH4)'
/

and here is the my  sectorlist_2018fd_nata:

sector,sectorcase,sectbaseyr,mrgapproach,prevyrspinup,endzip,speciation,mergesector,projectroot
"GID619",2018fd_nata_cb6_14j","2018","all","SectBaseYr","Y","","Y",""
"wRxBurn","2018fd_nata_cb6_14j","2018","all","SectBaseYr","Y","","N",""
"aircraft","2018fd_nata_cb6_14j","2018","all","SectBaseYr","Y","","N",""

Thank you,

Elham

[Gabriele Pfister]

Dear Elham

it seems you are using your own emission inventory rather than one that is provided by us.  Since we do not have access to this inventory we unfortunately cannot help you. In theory the tool should work with any SMOKE CMAQ ready output as long as you make the necessary adjustments to the namelist.

Gabriele

--

================================
Gabriele Pfister, Deputy Director

Atmospheric Chemistry Observations & Modeling

National Center for Atmospheric Research (NCAR)

Email: pfi...@ucar.edu

Temporary Home Office Phone: +1 303 241 1329

Work Phone: +1 303 497 2915

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