Error running convert_emiss with NEI11 emissions

Ideas, problems or suggestions regarding the initialization of WRF Chemistry

Error running convert_emiss with NEI11 emissions

Postby Tornerdo » Tue Apr 04, 2017 1:16 pm

I am trying to use NEI 2011 emissions in WRF-Chem 3.8.1. I have ran emiss_v04 and have generated the 2 binary files. However, when I run convert_emiss, I get the following error:

USING FILE: wrfem_00to12z_d01 1
Reading data from file: 1028663218
forrtl: severe (67): input statement requires too much data, unit 91, file /glade/scratch/srahimi/myWRFCHEM/WRFV3/test/em_real/wrfem_00to12z_d01
Image PC Routine Line Source
libirc.so 00002AC66DD70A7A Unknown Unknown Unknown
libirc.so 00002AC66DD6F576 Unknown Unknown Unknown
convert_emiss.exe 0000000003524E10 Unknown Unknown Unknown
convert_emiss.exe 00000000034C4D6E Unknown Unknown Unknown
convert_emiss.exe 00000000034C42AF Unknown Unknown Unknown
convert_emiss.exe 00000000034E9919 Unknown Unknown Unknown
convert_emiss.exe 00000000034E7704 Unknown Unknown Unknown
convert_emiss.exe 0000000000465A0B Unknown Unknown Unknown
convert_emiss.exe 000000000040BC4C Unknown Unknown Unknown
libc.so.6 00002AC66E57CD5D Unknown Unknown Unknown
convert_emiss.exe 000000000040BB49 Unknown Unknown Unknown

Here is the emiss_v04 file has the following user specifications:

PARAMETER(IX2=150,JX2=150,KX=44,KP=KX+1)
PARAMETER(IPRIM=7, IVOC=50, IPM25=19)
PARAMETER(NAL2DO=56, NRADM=31, NAMFILE2=64)
PARAMETER(KWIN=44,KWINP=KWIN+1)
INTEGER :: hemi = 1
INTEGER :: iproj = 1
REAL :: rekm = 6371.
REAL :: dx = 6.E3
REAL :: xlatc = 35.83001 !ref_lat in namelist.wps
REAL :: xlonc = -98.03 !ref_lon in namelist.wps
REAL :: clat1 = 60.0 !truelat2 in namelist.wps
REAL :: clat2 = 30.0 !truelat1 in namelist.wps
INTEGER :: inest1 = 0
REAL :: xnesstr = 1.00
REAL :: ynesstr = 1.00
REAL :: dxbigdo = 6.E3
INTEGER :: il = 151
INTEGER :: jl = 151
INTEGER :: starthr = 01
INTEGER :: maxhr = 24

I generated REFWZ and ZFA (which are nearly identical) from the geopotential and HGT_M information. The only difference between REFWZ and ZFA is that the REFWZ array starts at zero. My WSPD array was read in from an external file [using average windspeeds from a WRF (no Chem) case] in the AL2RADM2 subroutine.

My namelist is:

&time_control
run_days = 0,
run_hours = 23,
run_minutes = 0,
run_seconds = 0,
start_year = 2015, 2000, 2000,
start_month = 07, 01, 01,
start_day = 14, 24, 24,
start_hour = 00, 12, 12,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2015, 2000, 2000,
end_month = 07, 01, 01,
end_day = 14, 25, 25,
end_hour = 23, 12, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 3600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1, 1000, 1000,
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval_m = 360, 360
restart = .false.,
restart_interval = 5000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2

auxinput6_inname = 'wrfbiochemi_d01',
auxinput7_inname = 'wrffirechemi',
auxinput8_inname = 'wrfchemi_gocart_bg',
auxinput12_inname = 'wrf_chem_input',
auxinput5_interval_m = 60, 60, 60,
auxinput7_interval = 86400, 60, 60,
auxinput8_interval = 86400, 60, 60,
io_form_auxinput2 = 2,
io_form_auxinput5 = 2,
io_form_auxinput6 = 0,
io_form_auxinput7 = 0,
io_form_auxinput8 = 0,
io_form_auxinput12 = 0,

auxinput5_inname = "wrfchemi_d<domain>_<date>",
frames_per_auxinput5 = 12

debug_level = 200
/

&domains
time_step = 30,
time_step_fract_den = 1,
max_dom = 1,
e_we = 151, 112, 94,
e_sn = 151, 97, 91,
e_vert = 45, 30, 30,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 6000, 10000, 3333.33,
dy = 6000, 10000, 3333.33,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 31, 30,
j_parent_start = 1, 17, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/


&physics
mp_physics = 10, 3, 3,
ra_lw_physics = 4, 1, 1,
ra_sw_physics = 4, 1, 1,
radt = 0, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 5, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 3,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
/

&fdda
/

&dynamics
w_damping = 0,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&chem

kemit = 19,
chem_opt = 7,
bioemdt = 30,
photdt = 30,
chemdt = 5.00,
io_style_emissions = 1,
emiss_opt = 3, 3,
chem_in_opt = 0, 0,
phot_opt = 2, 0,
gas_drydep_opt = 1, 1,
aer_drydep_opt = 1, 1,
bio_emiss_opt = 1, 0,
dust_opt = 0,
dmsemis_opt = 0,
seas_opt = 0,
gas_bc_opt = 1, 0,
gas_ic_opt = 1, 0,
aer_bc_opt = 1, 0,
aer_ic_opt = 1, 0,
gaschem_onoff = 1, 0,
aerchem_onoff = 1, 0,
wetscav_onoff = 0, 0,
cldchem_onoff = 0, 0,
vertmix_onoff = 1, 0,
chem_conv_tr = 0, 0,
biomass_burn_opt = 0, 0,
plumerisefire_frq = 30, 0,
have_bcs_chem = .false., .false., .false.,
aer_ra_feedback = 0,
aer_op_opt = 0,
opt_pars_out = 0,
diagnostic_chem = 0,
/


What is going on here???

-Stefan
Tornerdo
 
Posts: 9
Joined: Wed Jun 17, 2015 9:25 pm
Location: Laramie, WY

Re: Error running convert_emiss with NEI11 emissions

Postby Tornerdo » Thu Apr 06, 2017 3:17 pm

Update to problem....

The wrfchemi_00z file is produced. However, I would like this for all hours in a day....The error occurs after this first file is produced.

-Stefan
Tornerdo
 
Posts: 9
Joined: Wed Jun 17, 2015 9:25 pm
Location: Laramie, WY

Re: Error running convert_emiss with NEI11 emissions

Postby Tornerdo » Thu Apr 13, 2017 6:08 pm

PROBLEM FIXED!!!!

I had the wrong ix2 and jx2 values. Check these....If you have any questions, please comment on this thread, and I can help you!

-Stefan
Tornerdo
 
Posts: 9
Joined: Wed Jun 17, 2015 9:25 pm
Location: Laramie, WY

Re: Error running convert_emiss with NEI11 emissions

Postby cloudphysics » Thu Oct 26, 2017 3:14 am

Dear, Stefan.

Hello! I'm having a same problem with you.
but I'm bit confused with the correct value for IX2 and JX2.

in your emiss_v04 file, you use il=151 and IX2=150.
isn't this correct value for IX2??

how did you change IX2=150 to what value?

Thank you in advance!
cloudphysics
 
Posts: 1
Joined: Thu Oct 26, 2017 3:06 am


Return to Initialization

Who is online

Users browsing this forum: No registered users and 1 guest

cron