Anthro_emis error = incrementing the emission inventory

[Luke Conibear]

Hi,

I have a recurring error of anthro_emis trying to increment the emission inventory.

When having this error in the past, I've been able to fix it via ensuring data_yrs_offset is set correctly, the date/datesec variables are set correctly, the inventory files have the correct units, the inventory files are readable/executable, and beginning anthro_emis for the 1st of the month. However, the error is persisting despite these. I've attached a simplified version of my anthro_emis.inp, along with the anthro_emis.out file, and a ncdump -h of one emission inventory file.

I'd be very grateful for any assistance you can provide.

Best,

Luke

[Stacy Walters]

Luke,

Would you please issue the follwing command:

ncdump -v date,datesec CAMS-GLOB-ANT_Glb_0.1x0.1_anthro_co_v4.1-R1.1_2016 > ncdump.log

and return the ncdump.log to me?

Stacy

[Luke Conibear]

Hi Stacy,

Sure, attached. I set the date as an integer for the start of the emission year e.g. date = 20160101.

Best,

Luke

[Stacy Walters]

Luke,

The date and datesec variables should be integer arrays, length time.

They should be date(time) and datesec(time).  The source file CAMS_GLOB...

looks like it represents a collection of "monthly" values.

You need to set the date and datesec arrays.

Stacy

[Luke Conibear]

Hi Stacy,

That worked, great. Thanks very much for your help.

Best,

Luke

[Caterina Mogno]

Dear Stacy, 

I am facing the same issue as Luke, with anthro_emiss trying to increment the emission inventory. I am using the EDGAR-HTAP inventory provided in the page WRF-Chem tools for the community.

I checked that the variable date and datesec are proper arrays, I also checked that the data_yrs_offset was correct, but still it doesn't work. 

I attach here my inputs, the log file of anthro_emis and the ncdump for the inventory.

Thanks in advance for help!

Best,

Caterina

[Stacy Walters]

Caterina,

You should change the variable data_yrs_offset from:

data_yrs_offset = -7

to

data_yrs_offset = 7

The variables data_yrs_offset is added to the time in the data file to produce the corresponding

time in the model simulation.  In this case, the data is for year 2010 whereas the simulation time

is for year 2017.

Stacy

[Ailish]

Hi Stacy,

I am also struggling to fix this same error (where anthro_emis is trying to increment the emission inventory) using the standard EDGAR-HTAP inventory in December 2019. I have run for  earlier months of the same year successfully though (August to December). 

I am running 2019-12-07 to 2019-12-14 so there should be no issue with a change in the data_yrs_offset. I've cheked and data_yrs_offset is set correctly (to 9) and the datesec (0,...,0) and data (20100101-20101201) arrays are arrays. 

Do you have any further suggestions of things I could try to fix this issue? I have attached my anthro_emis namelist and an ncdump of the emission file. 

Many Thanks,

Ailish

On Monday, January 27, 2020 at 4:27:16 PM UTC, Stacy Walters wrote:

Caterina,

You should change the variable data_yrs_offset from:

data_yrs_offset = -7

to

data_yrs_offset = 7

The variables data_yrs_offset is added to the time in the data file to produce the corresponding

time in the model simulation.  In this case, the data is for year 2010 whereas the simulation time

is for year 2017.

Stacy

[Rajesh Kumar]

Hi Ailish,

This is because the dates in the monthly files are set to the first day of the month. Because of this anthro_emis thinks that the last available date is 1 Dec. Thus, when you try to run anthro_emis for a later date in December as is your case, it cannot find the data in EDGAR files for that date. Since EDGAR emissions do not change from day to day in December, the easiest way would be to create two 12 hour wrfchemi files corresponding to io_style_emissions = 1 for 1 Dec and then use them in your runs. However, if you need hourly wrfchemi* files, you can either generate 24 files for 1 Dec and copy them to other days (don't forget to change the timestamp in wrfchemi* files for other days if you follow this method). Another way would be to change the last date in EDGAR-HTAP files to 20101231. 

Hope this helps,

Rajesh

--

Rajesh Kumar, Project Scientist II

Research Applications Laboratory 

National Center for Atmospheric Research 

E-mail: rku...@ucar.edu

Phone: +1-303-497-1413

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

[Ailish]

Hi Rajesh,

Thanks for the quick response. I thought this may be the case, I just wasn't sure of the solution to the problem. I've changed the last date as you suggested and anthro_emis has now successfully run. 

Thanks very much for your help, 

Ailish

[Catalyst26]

Hello Dr. Stacy,

I have generated hourly emissions files using anthro_emiss tool. I want to merge the hourly files into daily files. Can you please guide me how to achieve it?

Thanks and regards

Sulaymon

[Catalyst26]

Solved