Fwd: Anthro Error

[Saurabh Sinha]

---------- Forwarded message ---------
From: Saurabh Sinha <saurabh...@gmail.com>
Date: Sun, 7 Feb, 2021, 23:41
Subject: Anthro Error
To: <wrfchem-...@acom.ucar.edu>

Hi,

I am executing ./anthro_emis < MOZART-MOSAIC.inp to prepare the anthropogenic emission but ends with the following error:

 data_file_init: Initializing type for src emission file EDGAR_HTAP_emi_CO_2010.0.1x0.1.nc

 data_file_init:  nlon_src, nlat_src =         3600        1800
 data_file_init: data_dx,dx,has_area_map =    11117.5771       45000.0000     T
 data_file_init: file EDGAR_HTAP_emi_CO_2010.0.1x0.1.nc is a new grid
data_file_init: xcen_src(1,2)  =  -179.949996948242     -179.850006103516
data_file_init: xedge(1,2) =  -179.999992370605     -179.900001525879
data_file_init: dx =                     10
data_file_init: nlon_src,nlat_src =   3600   1800
data_file_init: ycen_src  =   89.8499984741211      89.9499969482422
data_file_init: yedge_src =   89.8999938964844      89.9999961853027
 ====================================================
 get_units: con_fac(1,2) =    3.59999996E+12   1.00000000E+09
 ====================================================

 next_flnm; trying to increment file EDGAR_HTAP_emi_CO_2010.0.1x0.1.nc
 next_flnm; il, iu =           30          30
 next_flnm; file_number =            1
 next_flnm; new file = EDGAR_HTAP_emi_CO_2010.0.1x0.0.nc
 data_file_timing_init : Failed to open /beegfs/cair/sk20ahk/EDGAR-HTAP/MOZART_MOSAIC/EDGAR_HTAP_emi_CO_2010.0.1x0.0.nc
 No such file or directory

STOP Netcdf error

Kindly help me,

I don't know why it's searching for EDGAR_HTAP_emi_CO_2010.0.1x0.0.nc??.

I have checked the MOZART-MOSAIC.inp file.

Saurabh Kumar

[Rajesh Kumar]

Hi Saurabh,

This is related to how you set start_output_time in the anthro_emis. namelist.  start_output_time has to be greater than or equal to the time in the EDGAR_HTAP January file and less than or equal to the time in  the EDGAR_HTAP December file. 

Best, 

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

[Saurabh Sinha]

Respected sir,

At present, I am using following in the MOZART-MOSAIC.inp

 start_output_time = '2017-12-07_00:00:00'
 stop_output_time = '2017-12-09_00:00:00'

I am getting the following error:

next_flnm; trying to increment file EDGAR_HTAP_emi_CO_2010.0.1x0.1.nc
 next_flnm; il, iu =           30          30
 next_flnm; file_number =            1
 next_flnm; new file = EDGAR_HTAP_emi_CO_2010.0.1x0.0.nc
 data_file_timing_init : Failed to open /beegfs/cair/sk20ahk/EDGAR-HTAP/MOZART_MOSAIC/EDGAR_HTAP_emi_CO_2010.0.1x0.0.nc
 No such file or directory

STOP Netcdf error

Kindly suggest me,

Regards,

Saurabh Kumar

[Rajesh Kumar]

Hi Saurabh,

Could you please share your anthro_emis namelist so that I can better understand your configuration? 

Thanks,

Rajesh

[Saurabh Sinha]

I am hereby sending the MOZART-MOSAIC.inp file.

kindly help me.

&CONTROL

 anthro_dir = '/EDGAR-HTAP/MOZART_MOSAIC'
 wrf_dir    = '/WRFV421/test/em_real'
 domains    = 4,
 src_file_prefix = 'EDGAR_HTAP_emi_'
 src_file_suffix = '_2010.0.1x0.1.nc'
 src_names = 'CO(28)','NOx(30)','SO2(64)','NH3(17)','BC(12)','OC(12)','PM2.5(1)','PM10(1)',
             'BIGALK(72)','BIGENE(56)','C2H4(28)','C2H5OH(46)','C2H6(30)','CH2O(30)',
             'CH3CHO(44)','CH3COCH3(58)','CH3OH(32)','MEK(72)','TOLUENE(92)','C3H6(42)','C3H8(44)',
             'BENZENE(78)','XYLENE(106)'
 sub_categories  = 'emis_tot'
 cat_var_prefix  = ''
 cat_var_suffix  = ''
 serial_output   = .false.

 start_output_time = '2017-12-07_00:00:00'
 stop_output_time = '2017-12-09_00:00:00'

 data_yrs_offset   = 7,
 emissions_zdim_stag = 1,
 output_interval = 86400,
 emis_map = 'CO->CO','NO->NOx','SO2->SO2','NH3->NH3','BIGALK->BIGALK','BIGENE->BIGENE',
            'C2H4->C2H4','C2H5OH->C2H5OH','C2H6->C2H6','CH2O->CH2O','CH3CHO->CH3CHO',
            'CH3COCH3->CH3COCH3','CH3OH->CH3OH','MEK->MEK','TOLUENE->TOLUENE',
            'C3H6->C3H6','C3H8->C3H8','NO2->0.0*NOx','ISOP->0.0*CO','C10H16->0.0*CO',
            'SULF->0.0*SO2','C2H2->0.00561790*CO','BENZENE->BENZENE','XYLENE->XYLENE',
            'GLY->0.0*CO','MACR->0.0*CO','MGLY->0.0*CO','MVK->0.0*CO',
            'HCOOH->0.0*CO','HONO->0.0*CO',
            'ECI(a)->0.2*BC','ECJ(a)->0.8*BC','ORGI(a)->0.2*OC','ORGJ(a)->0.8*OC','PM25I(a)->0.2*PM2.5',
            'PM25J(a)->0.8*PM2.5','PM_10(a)->PM10','SO4I(a)->0.0*PM10','SO4J(a)->0.0*PM10','NO3I(a)->0.0*PM10',
            'NO3J(a)->0.0*PM10','NH4I(a)->0.0*PM10','NH4J(a)->0.0*PM10','NAI(a)->0.0*PM10','NAJ(a)->0.0*PM10',
            'CLI(a)->0.0*PM10','CLJ(a)->0.0*PM10','CO_A->CO','CO_BB->0.0*CO','ORGI_A(a)->0.0*PM10',
            'ORGI_BB(a)->0.0*PM10','ORGJ_A(a)->0.0*PM10','ORGJ_BB(a)->0.0*PM10','VOCA->0.0*PM10','VOCBB->0.0*PM10'

/

Thanks,

Saurabh  Kumar

[Rajesh Kumar]

Hi Saurabh,

You are getting this problem because your start_output_time is after the last time in the input data, which is 01 Dec 2010. Since EDGAR-HTAP does not include daily variability, the emissions are going to be the same for all days of December. So, I suggest the following: 

1. set start_output_time = '2017-12-01_00:00:00'

2. Remove stop_output_time from the namelist

That should generate two 12 hour wrfchemi* files which you can use with any day of December in your simulation. 

cheers,

Rajesh

[Stephanie Martínez Morales]

Hello everyone,

I'm facing the same proble with the error:

 next_flnm; trying to increment file edgar_v5_2015_CO_0.1x0.1.nc
 next_flnm; il, iu =           24          24
 next_flnm; file_number =            1
 next_flnm; new file = edgar_v5_2015_CO_0.1x0.0.nc
 data_file_timing_init : Failed to open /home/WRF_CHEM_Tools/EDGAR-HTAP/MOZART_V5/edgar_v5_2015_CO_0.1x0.0.nc

 No such file or directory                                                      
STOP Netcdf error

I've tried the Rajesh suggestion but didn't work for me, and the problem persist.

Could you share a solution for this issue?

Thank you in advance.

- Stephanie M

[Ankan Sarkar]

Hi Stephanie,

This error probably came because of the date and datesec variables in your emission datasets. If you read those variables, I think end date of the emission file is 20101201. You have to change '20101201' to '20101231' if you want to simulate beyond 1 December. So, please try and let me know if it works for you.

With regards,

Ankan

--
To unsubscribe from this group and stop receiving emails from it, send an email to wrf-chem+u...@ucar.edu.