Format of AOD data is not supported with RRTMG aer_opt = 2

[saurabh4115]

Hi there,

I am trying to run a regional climate simulation using wrf with fnl data set and I would like to incorporate aerosol data in the simulation from MERRA2. I understand that using by aer_opt = 2, with RRTMG I should be able to provide the aod5502d data to wrf using a netcdf-wrf supported file via auxiliary input stream 15.

I have created the necessary files however, wrf returns the following error:

open_hist_w : error opening wrfout_d01_2021-01-05_00:00:00 for writing. 100
Timing for Writing wrfout_d01_2021-01-05_00:00:00 for domain 1: 0.00646 elapsed seconds
d01 2021-01-05_00:00:00 MMINLU error on input
1 input_wrf: wrf_get_next_time current_date: 2021-01-05_00:00:00 Status = -4
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 906
... Could not find matching time in input file wrfinput_d01
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0


I don't know Why I am getting this error. Should the auxiliary input stream 15 be used by real.exe? Is there a compilation step that I am missing to incorporate this data? Any additional information or advice you can provide on how I get wrf to read this data would be appreciated.

The relevant namelist options and metadata from my aux input file are below:

auxinput15_inname = 'wrfaerinp_d<domain>_<date>',
frames_per_auxinput15 = 23
auxinput15_interval_m = 360,
io_form_auxinput15 = 2,



ncdump -h wrfaerinp_d01_2021-01-05.nc
netcdf wrfaerinp_d01_2021-01-05 {
dimensions:
Time = 23 ;
DateStrLen = 19 ;
west_east = 349 ;
south_north = 318 ;
variables:
float Time(Time) ;
Time:long_name = "synthesized time coordinate from Times(time)" ;
Time:units = "minutes since 2021-01-05 00:00:00" ;
Time:_CoordinateAxisType = "Time" ;
Time:calendar = "noleap" ;
char Times(Time, DateStrLen) ;
float XLONG(south_north, west_east) ;
XLONG:FieldType = 104. ;
XLONG:MemoryOrder = "XY " ;
XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG:units = "degree_east" ;
XLONG:stagger = "" ;
XLONG:coordinates = "XLONG XLAT" ;
float XLAT(south_north, west_east) ;
XLAT:FieldType = 104. ;
XLAT:MemoryOrder = "XY " ;
XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT:units = "degree_north" ;
XLAT:stagger = "" ;
XLAT:coordinates = "XLONG XLAT" ;
double AOD5502D(Time, south_north, west_east) ;
AOD5502D:FieldType = 104. ;
AOD5502D:MemoryOrder = "XY " ;
AOD5502D:description = "Total aerosol optical depth at 550 nm" ;
AOD5502D:units = "" ;
AOD5502D:stagger = "" ;
AOD5502D:coordinates = "XLONG XLAT" ;

// global attributes:
:Title = "Aerosol Properties" ;
}

Thank you for your time.

 

[Ming Chen]

Hi,
Everything looks fine to me except that you may need to change

auxinput15_inname = 'wrfaerinp_d<domain>_<date>'

to

auxinput15_inname = 'wrfaerinp_d<domain>.<date>'

Please try and let me know whether it works.

 

[saurabh4115]

Dear Ming Chen,


Thank you for your time and reply.

The problem is solved now. I have used directly my AOD data file name in auxfile15_inname option because the file is in the same folder.

auxinput15_inname = "wrfaerinp_d01_2021-01-05.nc"


Thanks again.