WRF Chem Anthropogenic Emission Error

[Parker Yang]

Hi Everyone, I have been trying to run anthro emission using wrf_chem. I was able to create an emission file however, been stuck when trying to run wrf.exe with it. I get the following error: 

[kwyang0331@ln002 em_real]$ vi rsl.error.0000

      ra_lw_physics:      4

      ra_sw_physics:      4

     bl_pbl_physics:      2

  sf_sfclay_physics:      2

 sf_surface_physics:      2

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

  Domain # 1: dx =   900.000 m

WRF V4.6.0 MODEL

git commit 0a11865f97680fdd6865b278ea29d910e5db3ed7 1 file changed, 24 insertions(+), 24 deletions(-)

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

 Parent domain

 ids,ide,jds,jde            1         299           1         299

 ims,ime,jms,jme           -4          37          -4          26

 ips,ipe,jps,jpe            1          30           1          19

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

DYNAMICS OPTION: Eulerian Mass Coordinate

   alloc_space_field: domain            1,               49837048 bytes allocated

  med_initialdata_input: calling input_input

   Input data is acceptable to use: wrfinput_d01

 CURRENT DATE          = 2017-06-17_00:00:00

 SIMULATION START DATE = 2017-06-17_00:00:00

Timing for processing wrfinput file (stream 0) for domain        1:   16.26828 elapsed seconds

Max map factor in domain 1 =  1.00. Scale the dt in the model accordingly.

 D01: Time step                              =    5.000000     (s)

 D01: Grid Distance                          =   0.9000000     (km)

 D01: Grid Distance Ratio dt/dx              =    5.555556     (s/km)

 D01: Ratio Including Maximum Map Factor     =    5.556793     (s/km)

 D01: NML defined reasonable_time_step_ratio =    6.000000

 Climate GHG input from file from year         1765 to         2499

 CO2   range =    277.913000000000        579.264000000000      ppm

 N2O   range =    274.372000000000        359.798000000000      ppb

 CH4   range =    738.986000000000        997.311000000000      ppb

 CFC11 range =   0.000000000000000E+000  1.400000000000000E-002 ppt

 CFC12 range =   0.000000000000000E+000   2.88100000000000      ppt

 Normal ending of CAMtr_volume_mixing_ratio file

 GHG annual values from CAM trace gas file

 Year =         2017, Julian day =          168

 CO2   =   4.057855861263989E-004 volume mixing ratio

 N2O   =   3.297912028356296E-007 volume mixing ratio

 CH4   =   1.873907776629597E-006 volume mixing ratio

 CFC11 =   2.262880692754965E-010 volume mixing ratio

 CFC12 =   5.087825408042348E-010 volume mixing ratio

INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"

 LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          61 CATEGORIES           2 SEASONS WATER CATEGORY =           17 SNOW CATEGORY =           15

INITIALIZE THREE Noah LSM RELATED TABLES

Skipping over LUTYPE = USGS

 LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20 CATEGORIES

 INPUT SOIL TEXTURE CLASSIFICATION = STAS

 SOIL TEXTURE CLASSIFICATION = STAS FOUND          19 CATEGORIES

ThompMP: read qr_acr_qg_V4.dat     instead of computing

ThompMP: read qr_acr_qsV2.dat instead of computing

ThompMP: read freezeH2O.dat instead of computing

 mediation_integrate.G        1242 DATASET=HISTORY

 mediation_integrate.G        1243  grid%id            1  grid%oid            1

Timing for Writing wrfout_d01_2017-06-17_00:00:00 for domain        1:   47.73076 elapsed seconds

d01 2017-06-17_00:00:00  Input data is acceptable to use: wrfchemi_d01_2017-06-17_00:00:00

forrtl: severe (64): input conversion error, unit -5, file Internal Formatted Read

It seems that my wrfchemi file is acceptable but constant input conversion error. I was wondering if anyone is familiar with such error. Thank you in advance! 

[Gabriele Pfister]

Hi 

I am not familiar with the issue but could you share the emission file and your namelist.input for information?

Thanks

Gabi

--

================================
Gabriele Pfister

Atmospheric Chemistry Observations & Modeling (ACOM) Lab

NSF National Center for Atmospheric Research (NCAR)

Email: pfi...@ucar.edu

Work Phone: +1 303 497 2915

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

[Parker Yang]

Hi Gabi,

Thank you for the quick reply. I have attached the wrfchemi file and namelist.input to this email. I wanted to let you know that I have used the following code to generate my wrfchemi file: https://github.com/quishqa/PyChEmiss/blob/master/README.md. Please let me know if you need any other information. Again, thank you for your help.

Best,

Parker Yang

 namelist.input

 wrfchemi_d01_2017-06-17_00:00:00

[Gabriele Pfister]

Hi, 

I am not at all familiar with this tool but looking at the wrfchemi file, I see that the dimensions record is called emissions_zdim and not emissions_zdim_stag. It might be fine but what can cause an issue is that your emission_zdim is 1 in the wrfchemi file but you use kemit=19 in your namelist. Try setting kemit=1. 

Otherwise the wrfchemi file seems to have the needed fields included. The aerosol units are incorrect but I do not think that WRF-Chem cares about that. 

In your namelist I noticed another typos: 

 force_use_old_data                  = .trule., instead of .true.

It's seems a rather strange setup you are using and the emissions only cover a part of your domain, but I assume you have a very specific application in mind for this given that you also run at quite high resolution.

Hope this helps, 

Gabi