Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

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

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

saurabh4115

New member
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.
 
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.
 
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.
 
Top