[gradsusr] Converting and displaying a data subset
SteynAS at ufs.ac.za
Sun Jul 28 16:48:39 EDT 2019
I am using grib2ctl.pl and GrADS in order to convert ERA5 reanalysis data (e.g. 2m temperature at 12Z) from a .grib format to a .dat format.
The grib file contains global data (x=1440 and y=721 for a 0.25x0.25 degree grid).
When I use convert.gs to create the .dat file for the whole dataset everything works fine and the data is displayed correctly thereafter.
However, I don't have the computing power to ingest such a huge dataset and would rather work with a subset for a window area.
Now, when I run a modified convert.gs to create the .dat file for the smaller domain (e.g. x=96 to 128 and y=234 to 264) I don't get any error messages and everything seems to still work fine.
But upon displaying the data I get funny zonal isotherms for temp and constant field displays for wind.
Any ideas please?
Modified script for converting for smaller window area (original domain in comments):
'set x 96 128'
** was originally 'set x 1 1440'
** in order to set lon 0.00 359.75
'set y 234 264'
** was originally 'set y 1 721'
** in order to set lat -90.00 90.00
'set time 'tlabel''
'set gxout fwrite'
'set fwrite /run/media/stephan/SS_Transcend/ERA5/'fname2''
while (tyd <= 31)
'set t 'tyd
*** surface 10 metre U wind component [m/s]
*** surface 10 metre V wind component [m/s]
*** surface 2 metre dewpoint temperature [K]
*** surface 2 metre temperature [K]
University of the Free State: This message and its contents are subject to a disclaimer.
Please refer to http://www.ufs.ac.za/disclaimer for full details.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gradsusr