anthro_emis: failed to read namelist; error = 19
366 views
Skip to first unread message
Gonzalo Ferrada
Jun 27, 2019, 4:41:23 PM6/27/19
to wrf-chem-anthro_emiss
Hi everyone,
I am getting the following output when I run anthro_emis:
./anthro_emis < anthro_emis.inp
anthro_emis: failed to read namelist; error = 19
And no other hint is given. Has anyone experienced this issue?
I verified that wrfinput_d01 exists and the paths are correct.
I am using EPA2014v2 data downloaded from: https://www.acom.ucar.edu/Models/EPA/ and being sure that I have the sectorlist_2014fd_nata on it.
I am attaching my namelist:
&CONTROL
anthro_dir = '/Shared/CGRER-Data/gonzalo/EPA_EMISSIONS/cmaq_cb6'
wrf_dir = '/Shared/CGRER-Scratch/gonzalo/FIREX/FORECAST/8km_WRFCHEM_3.9_NCAR_tracers'
src_lon_dim_name = 'COL'
src_lat_dim_name = 'ROW'
domains = 1
sec_file_prefix = 'emis_mole_'
sec_file_suffix = '_12US2_nobeis_2014fd_nata_cb6_14j_nohap.nc4'
stk_file_prefix = 'inln_mole_'
stk_file_suffix = '_12US2_cmaq_cb6_2014fd_nata_cb6_14j_nohap.nc4'
stk_grp_file_suffix = '_12US2_2014fd_nata_cb6_14j_nohap.nc4'
currDateStr = '20140619'
src_names(1:6) = 'all(1):epa-sector','cmv_c3(1):epa-stack','othpt(1):epa-stack','ptegu(1):epa-stack','ptnonipm(1):epa-stack','pt_oilgas(1):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:26) = 'XYLMN','ISOP','TERP','SULF','ETHY','CH4'
sub_categories(27:33) = 'PMOTHR','PEC','POC','PMC','PSO4','PNO3','PNH4'
cat_var_prefix = ' '
serial_output = .true.
start_output_time = '2014-06-19_12:00:00'
stop_output_time = '2014-06-23_12:00:00'
data_yrs_offset = 0
emissions_zdim_stag = 10
emis_map(1) = 'CO->all(CO)+cmv_c3(CO)+othpt(CO)+ptegu(CO)+ptnonipm(CO)+pt_oilgas(CO)'
emis_map(2) = 'CH4->all(CH4)+cmv_c3(CH4)+othpt(CH4)+ptegu(CH4)+ptnonipm(CH4)+pt_oilgas(CH4)'
emis_map(3) = 'NO->all(NO)+cmv_c3(NO)+othpt(NO)+ptegu(NO)+ptnonipm(NO)+pt_oilgas(NO)'
emis_map(4) = 'NO2->all(NO2)+cmv_c3(NO2)+othpt(NO2)+ptegu(NO2)+ptnonipm(NO2)+pt_oilgas(NO2)'
emis_map(5) = 'SO2->all(SO2)+cmv_c3(SO2)+othpt(SO2)+ptegu(SO2)+ptnonipm(SO2)+pt_oilgas(SO2)'
emis_map(6) = 'NH3->all(NH3)+cmv_c3(NH3)+othpt(NH3)+ptegu(NH3)+ptnonipm(NH3)+pt_oilgas(NH3)'
emis_map(7) = 'C2H5OH->all(ETOH)+cmv_c3(ETOH)+othpt(ETOH)+ptegu(ETOH)+ptnonipm(ETOH)+pt_oilgas(ETOH)'
emis_map(8) = 'BIGALK->.2*all(PAR)+.2*cmv_c3(PAR)+.2*othpt(PAR)+.2*ptegu(PAR)+.2*ptnonipm(PAR)+.2*pt_oilgas(PAR)'
emis_map(9) = 'BIGENE->all(IOLE)+cmv_c3(IOLE)+othpt(IOLE)+ptegu(IOLE)+ptnonipm(IOLE)+pt_oilgas(IOLE)'
emis_map(10) = 'C2H4->all(ETH)+cmv_c3(ETH)+othpt(ETH)+ptegu(ETH)+ptnonipm(ETH)+pt_oilgas(ETH)'
emis_map(11) = 'C2H6->all(ETHA)+cmv_c3(ETHA)+othpt(ETHA)+ptegu(ETHA)+ptnonipm(ETHA)+pt_oilgas(ETHA)'
emis_map(12) = 'C3H6->all(OLE)+cmv_c3(OLE)+othpt(OLE)+ptegu(OLE)+ptnonipm(OLE)+pt_oilgas(OLE)'
emis_map(13) = 'C3H8->all(PRPA)+cmv_c3(PRPA)+othpt(PRPA)+ptegu(PRPA)+ptnonipm(PRPA)+pt_oilgas(PRPA)'
emis_map(14) = 'CH2O->all(FORM)+cmv_c3(FORM)+othpt(FORM)+ptegu(FORM)+ptnonipm(FORM)+pt_oilgas(FORM)'
emis_map(15) = 'CH3CHO>all(ALD2+ALDX)+cmv_c3(ALD2+ALDX)+othpt(ALD2+ALDX)+ptegu(ALD2+ALDX)+ptnonipm(ALD2+ALDX)+pt_oilgas(ALD2+ALDX)'
emis_map(16) = 'CH3COCH3->all(ACET)+cmv_c3(ACET)+othpt(ACET)+ptegu(ACET)+ptnonipm(ACET)+pt_oilgas(ACET)'
emis_map(17) = 'CH3OH->all(MEOH)+cmv_c3(MEOH)+othpt(MEOH)+ptegu(MEOH)+ptnonipm(MEOH)+pt_oilgas(MEOH)'
emis_map(18) = 'MEK->all(KET)+cmv_c3(KET)+othpt(KET)+ptegu(KET)+ptnonipm(KET)+pt_oilgas(KET)'
emis_map(19) = 'TOLUENE->all(TOL)+cmv_c3(TOL)+othpt(TOL)+ptegu(TOL)+ptnonipm(TOL)+pt_oilgas(TOL)'
emis_map(20) = 'BENZENE->all(BENZ)+cmv_c3(BENZ)+othpt(BENZ)+ptegu(BENZ)+ptnonipm(BENZ)+pt_oilgas(BENZ)'
emis_map(21) = 'XYLENE->all(XYLMN)+cmv_c3(XYLMN)+othpt(XYLMN)+ptegu(XYLMN)+ptnonipm(XYLMN)+pt_oilgas(XYLMN)'
emis_map(22) = 'ISOP->all(ISOP)+cmv_c3(ISOP)+othpt(ISOP)+ptegu(ISOP)+ptnonipm(ISOP)+pt_oilgas(ISOP)'
emis_map(23) = 'C10H16->all(TERP)+cmv_c3(TERP)+othpt(TERP)+ptegu(TERP)+ptnonipm(TERP)+pt_oilgas(TERP)'
emis_map(24) = 'sulf->all(SULF)+cmv_c3(SULF)+othpt(SULF)+ptegu(SULF)+ptnonipm(SULF)+pt_oilgas(SULF)'
emis_map(25) = 'C2H2->all(ETHY)+cmv_c3(ETHY)+othpt(ETHY)+ptegu(ETHY)+ptnonipm(ETHY)+pt_oilgas(ETHY)'
emis_map(26) = 'PM_25(A)>all(PMOTHR)+cmv_c3(PMOTHR)+othpt(PMOTHR)+ptegu(PMOTHR)+ptnonipm(PMOTHR)+pt_oilgas(PMOTHR)'
emis_map(27) = 'BC(A)->all(PEC)+cmv_c3(PEC)+othpt(PEC)+ptegu(PEC)+ptnonipm(PEC)+pt_oilgas(PEC)'
emis_map(28) = 'OC(A)->all(POC)+cmv_c3(POC)+othpt(POC)+ptegu(POC)+ptnonipm(POC)+pt_oilgas(POC)'
emis_map(29) = 'PM_10(A)->all(PMC)+cmv_c3(PMC)+othpt(PMC)+ptegu(PMC)+ptnonipm(PMC)+pt_oilgas(PMC)'
emis_map(30) = 'SO4I(A)->.15*all(PSO4)+.15*cmv_c3(PSO4)+.15*othpt(PSO4)+.15*ptegu(PSO4)+.15*ptnonipm(PSO4)+.15*pt_oilgas(PSO4)'
emis_map(31) = 'SO4J(A)->.85*all(PSO4)+.85*cmv_c3(PSO4)+.85*othpt(PSO4)+.85*ptegu(PSO4)+.85*ptnonipm(PSO4)+.85*pt_oilgas(PSO4)'
emis_map(32) = 'ECI(A)->.15*all(PEC)+.15*cmv_c3(PEC)+.15*othpt(PEC)+.15*ptegu(PEC)+.15*ptnonipm(PEC)+.15*pt_oilgas(PEC)'
emis_map(33) = 'ECJ(A)->.85*all(PEC)+.85*cmv_c3(PEC)+.85*othpt(PEC)+.85*ptegu(PEC)+.85*ptnonipm(PEC)+.85*pt_oilgas(PEC)'
emis_map(34) = 'ORGI(A)->.15*all(POC)+.15*cmv_c3(POC)+.15*othpt(POC)+.15*ptegu(POC)+.15*ptnonipm(POC)+.15*pt_oilgas(POC)'
emis_map(35) = 'ORGJ(A)->.85*all(POC)+.85*cmv_c3(POC)+.85*othpt(POC)+.85*ptegu(POC)+.85*ptnonipm(POC)+.85*pt_oilgas(POC)'
emis_map(36) = 'NO3I(A)->.15*all(PNO3)+.15*cmv_c3(PNO3)+.15*othpt(PNO3)+.15*ptegu(PNO3)+.15*ptnonipm(PNO3)+.15*pt_oilgas(PNO3)'
emis_map(37) = 'NO3J(A)->.85*all(PNO3)+.85*cmv_c3(PNO3)+.85*othpt(PNO3)+.85*ptegu(PNO3)+.85*ptnonipm(PNO3)+.85*pt_oilgas(PNO3)'
emis_map(38) = 'NH4I(A)->.15*all(PNH4)+.15*cmv_c3(PNH4)+.15*othpt(PNH4)+.15*ptegu(PNH4)+.15*ptnonipm(PNH4)+.15*pt_oilgas(PNH4)'
emis_map(39) = 'NH4J(A)->.85*all(PNH4)+.85*cmv_c3(PNH4)+.85*othpt(PNH4)+.85*ptegu(PNH4)+.85*ptnonipm(PNH4)+.85*pt_oilgas(PNH4)'
/
Gabriele Pfister
Jun 27, 2019, 5:21:18 PM6/27/19
to Gonzalo Ferrada, wrf-chem-anthro_emiss
Gonzalo
this indicates that you have some incorrect character or typos or similar in your namelist. One thing I can see after a quick look is that you have:
emis_map(15) = 'CH3CHO>all(ALD2+ALDX)+...
but it should be:
emis_map(15) = 'CH3CHO->all(ALD2+ALDX)+...
this indicates that you have some incorrect character or typos or similar in your namelist. One thing I can see after a quick look is that you have:
emis_map(15) = 'CH3CHO>all(ALD2+ALDX)+...
but it should be:
emis_map(15) = 'CH3CHO->all(ALD2+ALDX)+...
Gabi
================================
Gabriele Pfister, Deputy Director
Atmospheric Chemistry Observations & ModelingGabriele Pfister, Deputy Director
Suqian
Feb 9, 2022, 2:40:43 PM2/9/22
to wrf-chem-anthro_emiss, Gabriele Pfister, wrf-chem-anthro_emiss, gonzal...@gmail.com
Hi Gabi,
src_lon_dim_name = 'COL'
src_lat_dim_name = 'ROW'
domains = 1
sec_file_prefix = 'emis_mole_'
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'
I got the same issue and I found it would work only when I comment all the three lines:
currDateStr = '20170101'
serial_output = .true.
data_yrs_offset = 0
I am still confused about the solution to this problem. Do you have any idea why these three lines seem not work well? Is there any other solution?
Here is my input and it is based on the example input from the websie https://www.acom.ucar.edu/wrf-chem/anthro_emis_MOZCART-T1_EPA2017_web.inp
Many thanks in advance!
Suqian
&CONTROL
! Note: Chemical speciation is for MOZCRT-T1 but the mapping also included MOSAIC aerosol species
anthro_dir = '../DATA/EPA_2017'
wrf_dir = '../wrfinput/20170701_d01-EUS27km'
! Note: Chemical speciation is for MOZCRT-T1 but the mapping also included MOSAIC aerosol species
anthro_dir = '../DATA/EPA_2017'
wrf_dir = '../wrfinput/20170701_d01-EUS27km'
src_lon_dim_name = 'COL'
src_lat_dim_name = 'ROW'
domains = 1
sec_file_prefix = 'emis_mole_'
sec_file_suffix = '_12US1_nobeis_norwc_2017gb_17j.ncf'
stk_file_prefix = 'inln_mole_'
stk_file_suffix = '_12US1_cmaq_cb6ae7_2017gb_17j.ncf'
stk_grp_file_suffix = '_12US1_cmaq_cb6ae7_2017gb_17j.ncf'
!currDateStr = '20170101'
src_names(1:8) = 'all(1):epa-sector','rwc(1):epa-sector','cmv_c3_12(1):epa-stack','cmv_c1c2_12(1):epa-stack','othpt(1):epa-stack','ptegu(1):epa-stack','ptnonipm(1):epa-stack','pt_oilgas(1):epa-stack'
!src_names(1:4) = 'all(1):epa-sector','rwc(1):epa-sector','cmv_c3_12(1):epa-stack','cmv_c1c2_12(1):epa-stack'
!src_names(5:8) = 'othpt(1):epa-stack','ptegu(1):epa-stack','ptnonipm(1):epa-stack','pt_oilgas(1):epa-stack'
stk_file_prefix = 'inln_mole_'
stk_file_suffix = '_12US1_cmaq_cb6ae7_2017gb_17j.ncf'
stk_grp_file_suffix = '_12US1_cmaq_cb6ae7_2017gb_17j.ncf'
!currDateStr = '20170101'
src_names(1:8) = 'all(1):epa-sector','rwc(1):epa-sector','cmv_c3_12(1):epa-stack','cmv_c1c2_12(1):epa-stack','othpt(1):epa-stack','ptegu(1):epa-stack','ptnonipm(1):epa-stack','pt_oilgas(1):epa-stack'
!src_names(1:4) = 'all(1):epa-sector','rwc(1):epa-sector','cmv_c3_12(1):epa-stack','cmv_c1c2_12(1):epa-stack'
!src_names(5:8) = 'othpt(1):epa-stack','ptegu(1):epa-stack','ptnonipm(1):epa-stack','pt_oilgas(1):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:27) = 'XYLMN','ISOP','TERP','SULF','ETHY','HONO','ACROLEIN'
sub_categories(28:34) = 'PMOTHR','PEC','POC','PMC','PSO4','PNO3','PNH4'
cat_var_prefix = ' '
!serial_output = .true.
start_output_time = '2017-01-01_00:00:00'
stop_output_time = '2017-01-02_00:00:00'
!data_yrs_offset = 0
emissions_zdim_stag = 10
emis_map(1) = 'CO->all(CO)+cmv_c3_12(CO)+othpt(CO)+ptegu(CO)+cmv_c1c2_12(CO)+pt_oilgas(CO)+ptnonipm(CO)+rwc(CO)'
emis_map(2) = 'NO->all(NO)+cmv_c3_12(NO)+othpt(NO)+ptegu(NO)+cmv_c1c2_12(NO)+pt_oilgas(NO)+ptnonipm(NO)+rwc(NO)'
emis_map(3) = 'NO2->all(NO2)+cmv_c3_12(NO2)+othpt(NO2)+ptegu(NO2)+cmv_c1c2_12(NO2)+pt_oilgas(NO2)+ptnonipm(NO2)+rwc(NO2)'
emis_map(4) = 'SO2->all(SO2)+cmv_c3_12(SO2)+othpt(SO2)+ptegu(SO2)+cmv_c1c2_12(SO2)+pt_oilgas(SO2)+ptnonipm(SO2)+rwc(SO2)'
emis_map(5) = 'NH3->all(NH3)+cmv_c3_12(NH3)+othpt(NH3)+ptegu(NH3)+cmv_c1c2_12(NH3)+pt_oilgas(NH3)+ptnonipm(NH3)+rwc(NH3)'
emis_map(6) = 'C2H5OH->all(ETOH)+cmv_c3_12(ETOH)+othpt(ETOH)+ptegu(ETOH)+cmv_c1c2_12(ETOH)+pt_oilgas(ETOH)+ptnonipm(ETOH)+rwc(ETOH)'
emis_map(7) = 'BIGALK->.2*all(PAR)+.2*cmv_c3_12(PAR)+.2*othpt(PAR)+.2*ptegu(PAR)+.2*cmv_c1c2_12(PAR)+.2*pt_oilgas(PAR)+.2*ptnonipm(PAR)+.2*rwc(PAR)'
emis_map(8) = 'BIGENE->all(IOLE)+cmv_c3_12(IOLE)+othpt(IOLE)+ptegu(IOLE)+cmv_c1c2_12(IOLE)+pt_oilgas(IOLE)+ptnonipm(IOLE)+rwc(IOLE)'
emis_map(9) = 'C2H4->all(ETH)+cmv_c3_12(ETH)+othpt(ETH)+ptegu(ETH)+cmv_c1c2_12(ETH)+pt_oilgas(ETH)+ptnonipm(ETH)+rwc(ETH)'
emis_map(10) = 'C2H6->all(ETHA)+cmv_c3_12(ETHA)+othpt(ETHA)+ptegu(ETHA)+cmv_c1c2_12(ETHA)+pt_oilgas(ETHA)+ptnonipm(ETHA)+rwc(ETHA)'
emis_map(11) = 'C3H6->all(OLE)+cmv_c3_12(OLE)+othpt(OLE)+ptegu(OLE)+cmv_c1c2_12(OLE)+pt_oilgas(OLE)+ptnonipm(OLE)+rwc(OLE)'
emis_map(12) = 'C3H8->all(PRPA)+cmv_c3_12(PRPA)+othpt(PRPA)+ptegu(PRPA)+cmv_c1c2_12(PRPA)+pt_oilgas(PRPA)+ptnonipm(PRPA)+rwc(PRPA)'
emis_map(13) = 'CH2O->all(FORM)+cmv_c3_12(FORM)+othpt(FORM)+ptegu(FORM)+cmv_c1c2_12(FORM)+pt_oilgas(FORM)+ptnonipm(FORM)+rwc(FORM)'
emis_map(14) = 'CH3CHO->all(ALD2+ALDX)+cmv_c3_12(ALD2+ALDX)+othpt(ALD2+ALDX)+ptegu(ALD2+ALDX)+cmv_c1c2_12(ALD2+ALDX)+pt_oilgas(ALD2+ALDX)+ptnonipm(ALD2+ALDX)+rwc(ALD2+ALDX)'
emis_map(15) = 'CH3COCH3->all(ACET)+cmv_c3_12(ACET)+othpt(ACET)+ptegu(ACET)+cmv_c1c2_12(ACET)+pt_oilgas(ACET)+ptnonipm(ACET)+rwc(ACET)'
emis_map(16) = 'CH3OH->all(MEOH)+cmv_c3_12(MEOH)+othpt(MEOH)+ptegu(MEOH)+cmv_c1c2_12(MEOH)+pt_oilgas(MEOH)+ptnonipm(MEOH)+rwc(MEOH)'
emis_map(17) = 'MEK->all(KET)+cmv_c3_12(KET)+othpt(KET)+ptegu(KET)+cmv_c1c2_12(KET)+pt_oilgas(KET)+ptnonipm(KET)+rwc(KET)'
emis_map(18) = 'TOLUENE->all(TOL)+cmv_c3_12(TOL)+othpt(TOL)+ptegu(TOL)+cmv_c1c2_12(TOL)+pt_oilgas(TOL)+ptnonipm(TOL)+rwc(TOL)'
emis_map(19) = 'BENZENE->all(BENZ)+cmv_c3_12(BENZ)+othpt(BENZ)+ptegu(BENZ)+cmv_c1c2_12(BENZ)+pt_oilgas(BENZ)+ptnonipm(BENZ)+rwc(BENZ)'
emis_map(20) = 'XYLENE->all(XYLMN)+cmv_c3_12(XYLMN)+othpt(XYLMN)+ptegu(XYLMN)+cmv_c1c2_12(XYLMN)+pt_oilgas(XYLMN)+ptnonipm(XYLMN)+rwc(XYLMN)'
emis_map(21) = 'ISOP->all(ISOP)+cmv_c3_12(ISOP)+othpt(ISOP)+ptegu(ISOP)+cmv_c1c2_12(ISOP)+pt_oilgas(ISOP)+ptnonipm(ISOP)+rwc(ISOP)'
emis_map(22) = 'APIN->all(TERP)+cmv_c3_12(TERP)+othpt(TERP)+ptegu(TERP)+cmv_c1c2_12(TERP)+pt_oilgas(TERP)+ptnonipm(TERP)+rwc(TERP)'
emis_map(23) = 'sulf->all(SULF)+cmv_c3_12(SULF)+othpt(SULF)+ptegu(SULF)+cmv_c1c2_12(SULF)+pt_oilgas(SULF)+ptnonipm(SULF)+rwc(SULF)'
emis_map(24) = 'C2H2->all(ETHY)+cmv_c3_12(ETHY)+othpt(ETHY)+ptegu(ETHY)+cmv_c1c2_12(ETHY)+pt_oilgas(ETHY)+ptnonipm(ETHY)+rwc(ETHY)'
emis_map(25) = 'PM_25(A)->all(PMOTHR)+cmv_c3_12(PMOTHR)+othpt(PMOTHR)+ptegu(PMOTHR)+cmv_c1c2_12(PMOTHR)+pt_oilgas(PMOTHR)+ptnonipm(PMOTHR)+rwc(PMOTHR)'
emis_map(26) = 'BC(A)->all(PEC)+cmv_c3_12(PEC)+othpt(PEC)+ptegu(PEC)+cmv_c1c2_12(PEC)+pt_oilgas(PEC)+ptnonipm(PEC)+rwc(PEC)'
emis_map(27) = 'OC(A)->all(POC)+cmv_c3_12(POC)+othpt(POC)+ptegu(POC)+cmv_c1c2_12(POC)+pt_oilgas(POC)+ptnonipm(POC)+rwc(POC)'
emis_map(28) = 'PM_10(A)->all(PMC)+cmv_c3_12(PMC)+othpt(PMC)+ptegu(PMC)+cmv_c1c2_12(PMC)+pt_oilgas(PMC)+ptnonipm(PMC)+rwc(PMC)'
emis_map(29) = 'SO4I(A)->.15*all(PSO4)+.15*cmv_c3_12(PSO4)+.15*othpt(PSO4)+.15*ptegu(PSO4)+.15*cmv_c1c2_12(PSO4)+.15*pt_oilgas(PSO4)+.15*ptnonipm(PSO4)+.15*rwc(PSO4)'
emis_map(30) = 'SO4J(A)->.85*all(PSO4)+.85*cmv_c3_12(PSO4)+.85*othpt(PSO4)+.85*ptegu(PSO4)+.85*cmv_c1c2_12(PSO4)+.85*pt_oilgas(PSO4)+.85*ptnonipm(PSO4)+.85*rwc(PSO4)'
emis_map(31) = 'ECI(A)->.15*all(PEC)+.15*cmv_c3_12(PEC)+.15*othpt(PEC)+.15*ptegu(PEC)+.15*cmv_c1c2_12(PEC)+.15*pt_oilgas(PEC)+.15*ptnonipm(PEC)+.15*rwc(PEC)'
emis_map(32) = 'ECJ(A)->.85*all(PEC)+.85*cmv_c3_12(PEC)+.85*othpt(PEC)+.85*ptegu(PEC)+.85*cmv_c1c2_12(PEC)+.85*pt_oilgas(PEC)+.85*ptnonipm(PEC)+.85*rwc(PEC)'
emis_map(33) = 'ORGI(A)->.15*all(POC)+.15*cmv_c3_12(POC)+.15*othpt(POC)+.15*ptegu(POC)+.15*cmv_c1c2_12(POC)+.15*pt_oilgas(POC)+.15*ptnonipm(POC)+.15*rwc(POC)'
emis_map(34) = 'ORGJ(A)->.85*all(POC)+.85*cmv_c3_12(POC)+.85*othpt(POC)+.85*ptegu(POC)+.85*cmv_c1c2_12(POC)+.85*pt_oilgas(POC)+.85*ptnonipm(POC)+.85*rwc(POC)'
emis_map(35) = 'NO3I(A)->.15*all(PNO3)+.15*cmv_c3_12(PNO3)+.15*othpt(PNO3)+.15*ptegu(PNO3)+.15*cmv_c1c2_12(PNO3)+.15*pt_oilgas(PNO3)+.15*ptnonipm(PNO3)+.15*rwc(PNO3)'
emis_map(36) = 'NO3J(A)->.85*all(PNO3)+.85*cmv_c3_12(PNO3)+.85*othpt(PNO3)+.85*ptegu(PNO3)+.85*cmv_c1c2_12(PNO3)+.85*pt_oilgas(PNO3)+.85*ptnonipm(PNO3)+.85*rwc(PNO3)'
emis_map(37) = 'NH4I(A)->.15*all(PNH4)+.15*cmv_c3_12(PNH4)+.15*othpt(PNH4)+.15*ptegu(PNH4)+.15*cmv_c1c2_12(PNH4)+.15*pt_oilgas(PNH4)+.15*ptnonipm(PNH4)+.15*rwc(PNH4)'
emis_map(38) = 'NH4J(A)->.85*all(PNH4)+.85*cmv_c3_12(PNH4)+.85*othpt(PNH4)+.85*ptegu(PNH4)+.85*cmv_c1c2_12(PNH4)+.85*pt_oilgas(PNH4)+.85*ptnonipm(PNH4)+.85*rwc(PNH4)'
emis_map(39) = 'HONO->all(HONO)+cmv_c3_12(HONO)+othpt(HONO)+ptegu(HONO)+cmv_c1c2_12(HONO)+pt_oilgas(HONO)+ptnonipm(HONO)+rwc(HONO)'
emis_map(40) = 'MACR->all(ACROLEIN)+ptegu(ACROLEIN)+pt_oilgas(ACROLEIN)+ptnonipm(ACROLEIN)+rwc(ACROLEIN)'
/
sub_categories(28:34) = 'PMOTHR','PEC','POC','PMC','PSO4','PNO3','PNH4'
cat_var_prefix = ' '
!serial_output = .true.
start_output_time = '2017-01-01_00:00:00'
stop_output_time = '2017-01-02_00:00:00'
!data_yrs_offset = 0
emissions_zdim_stag = 10
emis_map(1) = 'CO->all(CO)+cmv_c3_12(CO)+othpt(CO)+ptegu(CO)+cmv_c1c2_12(CO)+pt_oilgas(CO)+ptnonipm(CO)+rwc(CO)'
emis_map(2) = 'NO->all(NO)+cmv_c3_12(NO)+othpt(NO)+ptegu(NO)+cmv_c1c2_12(NO)+pt_oilgas(NO)+ptnonipm(NO)+rwc(NO)'
emis_map(3) = 'NO2->all(NO2)+cmv_c3_12(NO2)+othpt(NO2)+ptegu(NO2)+cmv_c1c2_12(NO2)+pt_oilgas(NO2)+ptnonipm(NO2)+rwc(NO2)'
emis_map(4) = 'SO2->all(SO2)+cmv_c3_12(SO2)+othpt(SO2)+ptegu(SO2)+cmv_c1c2_12(SO2)+pt_oilgas(SO2)+ptnonipm(SO2)+rwc(SO2)'
emis_map(5) = 'NH3->all(NH3)+cmv_c3_12(NH3)+othpt(NH3)+ptegu(NH3)+cmv_c1c2_12(NH3)+pt_oilgas(NH3)+ptnonipm(NH3)+rwc(NH3)'
emis_map(6) = 'C2H5OH->all(ETOH)+cmv_c3_12(ETOH)+othpt(ETOH)+ptegu(ETOH)+cmv_c1c2_12(ETOH)+pt_oilgas(ETOH)+ptnonipm(ETOH)+rwc(ETOH)'
emis_map(7) = 'BIGALK->.2*all(PAR)+.2*cmv_c3_12(PAR)+.2*othpt(PAR)+.2*ptegu(PAR)+.2*cmv_c1c2_12(PAR)+.2*pt_oilgas(PAR)+.2*ptnonipm(PAR)+.2*rwc(PAR)'
emis_map(8) = 'BIGENE->all(IOLE)+cmv_c3_12(IOLE)+othpt(IOLE)+ptegu(IOLE)+cmv_c1c2_12(IOLE)+pt_oilgas(IOLE)+ptnonipm(IOLE)+rwc(IOLE)'
emis_map(9) = 'C2H4->all(ETH)+cmv_c3_12(ETH)+othpt(ETH)+ptegu(ETH)+cmv_c1c2_12(ETH)+pt_oilgas(ETH)+ptnonipm(ETH)+rwc(ETH)'
emis_map(10) = 'C2H6->all(ETHA)+cmv_c3_12(ETHA)+othpt(ETHA)+ptegu(ETHA)+cmv_c1c2_12(ETHA)+pt_oilgas(ETHA)+ptnonipm(ETHA)+rwc(ETHA)'
emis_map(11) = 'C3H6->all(OLE)+cmv_c3_12(OLE)+othpt(OLE)+ptegu(OLE)+cmv_c1c2_12(OLE)+pt_oilgas(OLE)+ptnonipm(OLE)+rwc(OLE)'
emis_map(12) = 'C3H8->all(PRPA)+cmv_c3_12(PRPA)+othpt(PRPA)+ptegu(PRPA)+cmv_c1c2_12(PRPA)+pt_oilgas(PRPA)+ptnonipm(PRPA)+rwc(PRPA)'
emis_map(13) = 'CH2O->all(FORM)+cmv_c3_12(FORM)+othpt(FORM)+ptegu(FORM)+cmv_c1c2_12(FORM)+pt_oilgas(FORM)+ptnonipm(FORM)+rwc(FORM)'
emis_map(14) = 'CH3CHO->all(ALD2+ALDX)+cmv_c3_12(ALD2+ALDX)+othpt(ALD2+ALDX)+ptegu(ALD2+ALDX)+cmv_c1c2_12(ALD2+ALDX)+pt_oilgas(ALD2+ALDX)+ptnonipm(ALD2+ALDX)+rwc(ALD2+ALDX)'
emis_map(15) = 'CH3COCH3->all(ACET)+cmv_c3_12(ACET)+othpt(ACET)+ptegu(ACET)+cmv_c1c2_12(ACET)+pt_oilgas(ACET)+ptnonipm(ACET)+rwc(ACET)'
emis_map(16) = 'CH3OH->all(MEOH)+cmv_c3_12(MEOH)+othpt(MEOH)+ptegu(MEOH)+cmv_c1c2_12(MEOH)+pt_oilgas(MEOH)+ptnonipm(MEOH)+rwc(MEOH)'
emis_map(17) = 'MEK->all(KET)+cmv_c3_12(KET)+othpt(KET)+ptegu(KET)+cmv_c1c2_12(KET)+pt_oilgas(KET)+ptnonipm(KET)+rwc(KET)'
emis_map(18) = 'TOLUENE->all(TOL)+cmv_c3_12(TOL)+othpt(TOL)+ptegu(TOL)+cmv_c1c2_12(TOL)+pt_oilgas(TOL)+ptnonipm(TOL)+rwc(TOL)'
emis_map(19) = 'BENZENE->all(BENZ)+cmv_c3_12(BENZ)+othpt(BENZ)+ptegu(BENZ)+cmv_c1c2_12(BENZ)+pt_oilgas(BENZ)+ptnonipm(BENZ)+rwc(BENZ)'
emis_map(20) = 'XYLENE->all(XYLMN)+cmv_c3_12(XYLMN)+othpt(XYLMN)+ptegu(XYLMN)+cmv_c1c2_12(XYLMN)+pt_oilgas(XYLMN)+ptnonipm(XYLMN)+rwc(XYLMN)'
emis_map(21) = 'ISOP->all(ISOP)+cmv_c3_12(ISOP)+othpt(ISOP)+ptegu(ISOP)+cmv_c1c2_12(ISOP)+pt_oilgas(ISOP)+ptnonipm(ISOP)+rwc(ISOP)'
emis_map(22) = 'APIN->all(TERP)+cmv_c3_12(TERP)+othpt(TERP)+ptegu(TERP)+cmv_c1c2_12(TERP)+pt_oilgas(TERP)+ptnonipm(TERP)+rwc(TERP)'
emis_map(23) = 'sulf->all(SULF)+cmv_c3_12(SULF)+othpt(SULF)+ptegu(SULF)+cmv_c1c2_12(SULF)+pt_oilgas(SULF)+ptnonipm(SULF)+rwc(SULF)'
emis_map(24) = 'C2H2->all(ETHY)+cmv_c3_12(ETHY)+othpt(ETHY)+ptegu(ETHY)+cmv_c1c2_12(ETHY)+pt_oilgas(ETHY)+ptnonipm(ETHY)+rwc(ETHY)'
emis_map(25) = 'PM_25(A)->all(PMOTHR)+cmv_c3_12(PMOTHR)+othpt(PMOTHR)+ptegu(PMOTHR)+cmv_c1c2_12(PMOTHR)+pt_oilgas(PMOTHR)+ptnonipm(PMOTHR)+rwc(PMOTHR)'
emis_map(26) = 'BC(A)->all(PEC)+cmv_c3_12(PEC)+othpt(PEC)+ptegu(PEC)+cmv_c1c2_12(PEC)+pt_oilgas(PEC)+ptnonipm(PEC)+rwc(PEC)'
emis_map(27) = 'OC(A)->all(POC)+cmv_c3_12(POC)+othpt(POC)+ptegu(POC)+cmv_c1c2_12(POC)+pt_oilgas(POC)+ptnonipm(POC)+rwc(POC)'
emis_map(28) = 'PM_10(A)->all(PMC)+cmv_c3_12(PMC)+othpt(PMC)+ptegu(PMC)+cmv_c1c2_12(PMC)+pt_oilgas(PMC)+ptnonipm(PMC)+rwc(PMC)'
emis_map(29) = 'SO4I(A)->.15*all(PSO4)+.15*cmv_c3_12(PSO4)+.15*othpt(PSO4)+.15*ptegu(PSO4)+.15*cmv_c1c2_12(PSO4)+.15*pt_oilgas(PSO4)+.15*ptnonipm(PSO4)+.15*rwc(PSO4)'
emis_map(30) = 'SO4J(A)->.85*all(PSO4)+.85*cmv_c3_12(PSO4)+.85*othpt(PSO4)+.85*ptegu(PSO4)+.85*cmv_c1c2_12(PSO4)+.85*pt_oilgas(PSO4)+.85*ptnonipm(PSO4)+.85*rwc(PSO4)'
emis_map(31) = 'ECI(A)->.15*all(PEC)+.15*cmv_c3_12(PEC)+.15*othpt(PEC)+.15*ptegu(PEC)+.15*cmv_c1c2_12(PEC)+.15*pt_oilgas(PEC)+.15*ptnonipm(PEC)+.15*rwc(PEC)'
emis_map(32) = 'ECJ(A)->.85*all(PEC)+.85*cmv_c3_12(PEC)+.85*othpt(PEC)+.85*ptegu(PEC)+.85*cmv_c1c2_12(PEC)+.85*pt_oilgas(PEC)+.85*ptnonipm(PEC)+.85*rwc(PEC)'
emis_map(33) = 'ORGI(A)->.15*all(POC)+.15*cmv_c3_12(POC)+.15*othpt(POC)+.15*ptegu(POC)+.15*cmv_c1c2_12(POC)+.15*pt_oilgas(POC)+.15*ptnonipm(POC)+.15*rwc(POC)'
emis_map(34) = 'ORGJ(A)->.85*all(POC)+.85*cmv_c3_12(POC)+.85*othpt(POC)+.85*ptegu(POC)+.85*cmv_c1c2_12(POC)+.85*pt_oilgas(POC)+.85*ptnonipm(POC)+.85*rwc(POC)'
emis_map(35) = 'NO3I(A)->.15*all(PNO3)+.15*cmv_c3_12(PNO3)+.15*othpt(PNO3)+.15*ptegu(PNO3)+.15*cmv_c1c2_12(PNO3)+.15*pt_oilgas(PNO3)+.15*ptnonipm(PNO3)+.15*rwc(PNO3)'
emis_map(36) = 'NO3J(A)->.85*all(PNO3)+.85*cmv_c3_12(PNO3)+.85*othpt(PNO3)+.85*ptegu(PNO3)+.85*cmv_c1c2_12(PNO3)+.85*pt_oilgas(PNO3)+.85*ptnonipm(PNO3)+.85*rwc(PNO3)'
emis_map(37) = 'NH4I(A)->.15*all(PNH4)+.15*cmv_c3_12(PNH4)+.15*othpt(PNH4)+.15*ptegu(PNH4)+.15*cmv_c1c2_12(PNH4)+.15*pt_oilgas(PNH4)+.15*ptnonipm(PNH4)+.15*rwc(PNH4)'
emis_map(38) = 'NH4J(A)->.85*all(PNH4)+.85*cmv_c3_12(PNH4)+.85*othpt(PNH4)+.85*ptegu(PNH4)+.85*cmv_c1c2_12(PNH4)+.85*pt_oilgas(PNH4)+.85*ptnonipm(PNH4)+.85*rwc(PNH4)'
emis_map(39) = 'HONO->all(HONO)+cmv_c3_12(HONO)+othpt(HONO)+ptegu(HONO)+cmv_c1c2_12(HONO)+pt_oilgas(HONO)+ptnonipm(HONO)+rwc(HONO)'
emis_map(40) = 'MACR->all(ACROLEIN)+ptegu(ACROLEIN)+pt_oilgas(ACROLEIN)+ptnonipm(ACROLEIN)+rwc(ACROLEIN)'
/
Stacy Walters
Feb 9, 2022, 4:24:06 PM2/9/22
to Suqian, wrf-chem-anthro_emiss, Gabriele Pfister, gonzal...@gmail.com
Suqian,
Can you please send me the anthro_emis.inp you are using as an attachment?
Stacy
0 new messages