Re: anthro_emis issue

[Gabriele Pfister]

Hello Gonzalo, 

Below you find an example namelist for EPA_anthro_emis. You will have to adapt it for your configuration and chemical and aerosol schemes. 

Please also send all future requests to our Discussion Forum (https://www2.acom.ucar.edu/wrf-chem/discussion-forum) to be shared with the community. 

Gabi

&CONTROL
anthro_dir = '/data14a/pfister/EPA_2014/cmaq_cb6'
wrf_dir    = '/home/pfister/ANTHRO_EPA/WRF'
src_lon_dim_name = 'COL'
src_lat_dim_name = 'ROW'
domains = 2
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     = '20140810'
src_names(1:6) = 'all(1):epa-sector','cmv_c3(1):epa-stack','othpt(1):epa-stack','ptegu(1):epastack','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-08-10_06:00:00'
stop_output_time  = '2014-08-13_06: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)'

On Wed, Jun 26, 2019 at 4:24 PM Ferrada, Gonzalo A <gonzalo...@uiowa.edu> wrote:

Hi Gabi, Stacy,

I just downloaded the anthro_emis tool to process EPA 2014 v2 data for WRF inputs from this link:

https://www.acom.ucar.edu/wrf-chem/download.shtml

In the compressed file only comes a directory called src with the source code. However, in the README file you talk about a namelist file that I assume it’s not the same that I am using for my WRF simulations since it says that I have to specify some variables that I have not seen before:

******************************************************************************************************

*** The following variables must be set in the input namelist file or anthro_emis will error halt. ***

******************************************************************************************************

        start_output_time

        src_names

        sub_categories

        emis_map

        If you want to use EPA 2D sector emissions you will need to set the following:

        sec_file_prefix

        sec_file_suffix

        If you want to use EPA point emissions you will need to set the following:

        stk_file_prefix

        stk_file_suffix

        stk_grp_file_suffix

        Although none of the remaining namelist variables need be set it is most likely you will 

        want to set some of the remaining variables.

My guess is that you forgot to include a sample of that namelist in the compressed file “EPA_ANTHRO_EMIS.tgz”. See in the image attached the contents of it when uncompressed.

Can you please send me a sample namelist?

Thanks,

--
Gonzalo A. Ferrada, Ph.D. (c)
Center for Global and Regional Environmental Research
Department of Chemical and Biochemical Engineering
The University of Iowa

--

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

Atmospheric Chemistry Observations & Modeling

National Center for Atmospheric Research

Email: pfi...@ucar.edu

Phone: +1 303 497 2915

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