[gradsusr] ECMWF open data in GrADS

[Bill_Reilly]

Davide:

I found an easier way using CDO by getting CDO to rename the variables, saving me from having to create and then edit a descriptor file:

C:\cygwin64\bin\cdo -f nc -select,name=10u,10v,msl C:\Progra~1\PCGrADS\ecmwf\grib\ecmwf_000.grib2 C:\Progra~1\PCGrADS\win32e\ecmwf_000_temp.nc
C:\cygwin64\bin\cdo -chname,10u,u10 -chname,10v,v10 C:\Progra~1\PCGrADS\win32e\ecmwf_000_temp.nc C:\Progra~1\PCGrADS\win32e\ecmwf_000.nc

Problem solved!... Thanks again for your reply.

Bill

--

Bill Reilly

(+34) 686-110-748

bill at passageweather.com
www.passageweather.com



 

    On Monday, October 16, 2023 at 04:02:36 PM GMT+2, Bill_Reilly <bill_reilly at compuserve.com> wrote:  
 
 Davide:

Thank you for your reply, and for the suggestion... That is exactly what I ended up doing.
I used CDO to only extract the surface variables that I need (10u,10v,msl):
cdo -f nc select,name=10u,10v,msl C:/Progra~1/PCGrADS/win32e/20231015120000-9h-oper-fc.grib2 C:/Progra~1/PCGrADS/win32e/ecmwf_sfc_009.nc

Then I had CDO create a GrADS control file:
 
cdo gradsdes C:/Progra~1/PCGrADS/win32e/ecmwf_sfc_009.nc

I then have to manually rename the variables (and add the z dimension for the wind variables) in the control file:

VARS  3
msl               0  t,y,x  Mean sea level pressure  [Pa]
10v=>v10          0  t,z,y,x  10 metre V wind component  [m s**-1]
10u=>u10          0  t,z,y,x  10 metre U wind component  [m s**-1]
ENDVARS

It now works perfectly!

Thanks again for your help.

Bill
--Bill Reilly

(+34) 686-110-748

bill at passageweather.com
www.passageweather.com




  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20231016/65aef144/attachment.html>