EPA anthro_emis error messages

[Lena Low]

Hello everyone, 

After downloading EPA 2014v2 NEI input data, I used both the README file and an example namelist script that I found in this User Forum to write two slightly different scripts based on the two different sources. I attached them in the file below. When I run EPA anthro_emis, I get these error messages, each from a different script: 

anthro_emis: failed to read namelist; error =           18

anthro_emis: failed to read namelist; error =           19

Does anyone have any experience with these error messages and know what they mean? I am wondering if the numbers correspond with the line in the makelist file or something else entirely? 

Thank you,

Lena

[Stacy Walters]

Lena,

Please email one of the input namelist files you are using.

From "googling" for fortran namelist error codes I get the following for codes 18, 19:

18

severe (18): Too many values for NAMELIST variable FOR$IOS_TOOMANVAL. An attempt was made to assign too many values to a variable during a namelist READ statement.

19

severe (19): Invalid reference to variable in NAMELIST input FOR$IOS_INVREFVAR. One of the following conditions occurred: The variable was not a member of the namelist group. An attempt was made to subscript a scalar variable. A subscript of the array variable was out-of-bounds. An array variable was specified with too many or too few subscripts for the variable. An attempt was made to specify a substring of a non-character variable or array name. A substring specifier of the character variable was out-of-bounds. A subscript or substring specifier of the variable was not an integer constant. An attempt was made to specify a substring by using an unsubscripted array variable

[Zhiqi Yang]

Hi Lena,

Do you solve your problem? I have the same problem. Could you give me some suggestions? Here is my input:

&CONTROL

 anthro_dir = '/Users/WRFC/NEI2017'

 wrf_dir    = '/Users/WRFC/WRFchemv3/UTILINP'

 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     = '20170322'

 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'

 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-03-22_00:00:00'

 stop_output_time  = '2017-03-23_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)'

/

Thank you,

Zhiqi

[Stacy Walters]

Zhiqi,

Have you tried to break up the line starting with 'src_names(1:8)=...' into two lines:

src_names(1:4)= '...'

src_names(5:8)='...'

Stacy

[Zhiqi Yang]

Hi Stacy,

It still doesn't work. It shows:

 anthro_emis: failed to read namelist; error =           19

Zhiqi

[Suqian]

Hello Zhiqi,

I got the same issue and my input is almost the same as yours. Have you figured it out? 

Any help would be much appreciated!

Cheers,

Suqian