Hi all,
I’m encountering an error when running the matching_radiation.csh script on WRF output files with a 1 km grid size. Interestingly, the matching_surface.csh and matching_raob.csh scripts ran successfully on the same set of WRF output files.
It seems the script is unable to download the SURFRAD data file and then encounters an array indexing error (subscript out of bounds)—possibly due to missing observation data.
Below are the matching_radiation.csh log file and the header info from the WRF output file:
matching.radiation.wrfDC1_Feb.log.txt (3.6 KB)
wrfout_d03_header.txt (39.1 KB)
Thanks in advance for your help!
Hello @hhallaji ,
what is the spatial extent of your 1 km WRF domain? Are any of the 7 SURFRAD sites listed in the station masterfile (/scratch/hhallaji/AMET_v15/obs/MET/bsrn/surfrad_sites.csv, with an online copy from the AMET code repository here) located within it?
I’m not familiar with this matching code, but the error message makes me think that none of the seven sites are within your domain, so the code just looks for a station “NA” to download and then fails when the download fails and no observational data is available for matching.
1 Like
@hogrefe.christian,
Thanks for the response.
I checked the surfrad_sites.csv file, and it seems the listed SURFRAD sites cover the entire U.S., while my finest WRF domain is much smaller. Specifically, it’s a 77×77 grid with a 1 km resolution, centered at:
CEN_LAT = 38.86612
CEN_LON = -77.07074
Given this limited spatial extent, it looks like none of the SURFRAD sites fall within the domain, which would explain why the script attempts to use a generic “NA” station and then fails when it can’t download the corresponding data file.