[gradsusr] ECMWF Data Issues

stuart markham stumark at hotmail.co.uk
Fri Mar 25 12:53:14 EDT 2016


Hi,

How do you download the data. I have used your scripts before and they are great. Would it be possible to view your script for downloading the ecm data?

Best regards,

Stuart

Sent from Outlook Mobile<https://aka.ms/blhgte>



On Fri, Mar 25, 2016 at 9:46 AM -0700, "AMS EmsiWx Bytča" <emsiwx at pocasie-bytca.sk<mailto:emsiwx at pocasie-bytca.sk>> wrote:

Thanks Stuart,

I have been processing ECMWF for some time, allready. I am glad to hear, all is ok now.

Take care, have a nice Easter ;)

Marian

Dňa 25. 3. 2016 o 17:01 stuart markham napísal(a):

Thanks Marian,


Adjusting the pres did the trick thanks. Out of interest i am downloading via this link ftp://data-portal.ecmwf.int/ but using an script you designed to download the gfs with, example below. Would you have an easier way of downloading data daily? Ive only come across this today. Many thanks



:024
wget -O ecmwf_24.grb ftp://wmo:essential@data-portal.ecmwf.int/20160325000000/A_HHXE50ECMF250000_C_ECMF_20160325000000_24h_gh_500hPa_global_0p5deg_grib2.bin ftp://wmo:essential@data-portal.ecmwf.int/20160325000000/A_HPXE89ECMF250000_C_ECMF_20160325000000_24h_msl_global_0p5deg_grib2.bin
IF %ERRORLEVEL% NEQ 0 @ECHO Waiting for file... & TIMEOUT /T 180 & GOTO 024



________________________________
From: gradsusr-bounces at gradsusr.org<mailto:gradsusr-bounces at gradsusr.org> <gradsusr-bounces at gradsusr.org><mailto:gradsusr-bounces at gradsusr.org> on behalf of AMS EmsiWx Bytča <emsiwx at pocasie-bytca.sk><mailto:emsiwx at pocasie-bytca.sk>
Sent: 25 March 2016 14:06
To: GrADS Users Forum
Subject: Re: [gradsusr] ECMWF Data Issues

Stuart,

This one does not seem correct to me:

'define PRESmsl = PRESmsl*0.01'

replace with e.g.

'define PRES = PRESmsl*0.01'

I am not sure if that helps you to get all time steps, but for sure on t=2 you would not get correct pressure contours.

Let me know how it goes, and please show me how you get the ecmwf grib.

M.

Dňa 25. 3. 2016 o 13:33 stuart markham napísal(a):

Hi Marian,


I've attached the full script for you, many thanks for your help.


regards,


Stuart.


________________________________
From: gradsusr-bounces at gradsusr.org<mailto:gradsusr-bounces at gradsusr.org> <gradsusr-bounces at gradsusr.org><mailto:gradsusr-bounces at gradsusr.org> on behalf of AMS EmsiWx Bytča <emsiwx at pocasie-bytca.sk><mailto:emsiwx at pocasie-bytca.sk>
Sent: 25 March 2016 11:56
To: GrADS Users Forum
Subject: Re: [gradsusr] ECMWF Data Issues

Hi Stuart,

in your script you are opening a GFS control file... Replace with (I suppose) ecmwf_grib.ctl . And I can't see where you change the t (time), Can you please show me the rest of your script?

Regards,

Marian


Dňa 25. 3. 2016 o 12:11 stuart markham napísal(a):

Hi all,


I'm currently developing charts using ecmwf free data  but I'm having an issue with the time steps. The first image is fine but the rest of the time steps use the same data as the first one. The charts are not running in sequence. Can anybody advise where my issues are. I've attached parts of my control file and script to assist.



dset ^ecmwf_grib.grib2
index ^ecmwf_grib.grib2.idx
undef 9.999E+20
title ecmwf_grib.grib2
* produced by g2ctl v0.0.8.8
* command line options: -verf ecmwf_grib.grib2
* griddef=1:0:(720 x 361):grid_template=0:winds(N/S): lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 0.500000 lon 0.000000 to 359.500000 by 0.500000 #points=259920:winds(N/S)

dtype grib2
ydef 361 linear -90.000000 0.5
xdef 720 linear 0.000000 0.500000
tdef 11 linear 00Z25mar2016 1dy
zdef 1 linear 1 1
vars 2
HGT500mb   0,100,50000   0,3,5 ** 500 mb Geopotential Height [gpm]
PRESmsl  0,101   0,3,0 ** mean sea level Pressure [Pa]
ENDVARS


----------------------------------------------------------------------------------------------------------------------------------------------------------------




'reinit'
'open gfs_grib.ctl'

* Get info from the descriptor file
'q ctlinfo'
_ctl = result
_undef = getctl(undef)
_tdef = getctl(tdef)
_zdef = getctl(zdef)

maps = 25

* Get the Time axis info
tsize = subwrd(_tdef,2)
_t1 = 1       ;* 2nd half of timeseries
_t2 = 45
tsize = _t2 - _t1 + 1
'set t '_t1' '_t2
'q dims'
times  = sublin(result,5)
_time1 = subwrd(times,6)
_time2 = subwrd(times,8)
_tdim = _time1' '_time2

tincr = subwrd(_tdef,5)
_tdef = 'tdef 'tsize' linear '_time1' 'tincr

huh = subwrd(_tdef,4)

Best regards,

Stuart.




_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr





_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr





_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160325/33e6dc33/attachment-0001.html 


More information about the gradsusr mailing list