<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Oh, I missed, this is the way I download under Win:<br>
<br>
@ECHO OFF<br>
set FileDate=%date:~9,4%%date:~6,2%%date:~3,2%<br>
set FileTime=%time:~0,2%<br>
<br>
for /f "delims=" %%a in ('wmic OS Get localdatetime ^| find "."')
do set "dt=%%a"<br>
set "YYYY=%dt:~0,4%"<br>
set "MM=%dt:~4,2%"<br>
set "DD=%dt:~6,2%"<br>
set "HH=%dt:~8,2%"<br>
set "Min=%dt:~10,2%"<br>
set "Sec=%dt:~12,2%"<br>
<br>
set FileDate=%YYYY%%MM%%DD%<br>
set timestamp=%HH%%Min%%Sec%<br>
set FileTime=%HH%<br>
set fullstamp=%YYYY%-%MM%-%DD%_%HH%-%Min%-%Sec%<br>
<br>
set FD=%FileDate%<br>
if /i %FileTime% gtr 06 (set T=00)<br>
if /i %FileTime% gtr 18 (set T=12)<br>
@echo Year = %FD%<br>
@echo Time = %FileTime%<br>
echo %T%<br>
<br>
wget --ftp-user=wmo --ftp-password=essential
ftp://data-portal.ecmwf.int/%FD%%T%0000/*%FD%%T%0000_an_msl_global_0p5deg_grib2.bin<br>
rename *%FD%%T%0000_an_msl_global_0p5deg_grib2.bin
%FD%%T%00000000_C_ECMF_%FD%%T%0000_0h_msl_global_0p5deg_grib2.bin<br>
wget --ftp-user=wmo --ftp-password=essential
ftp://data-portal.ecmwf.int/%FD%%T%0000/*%FD%%T%0000*h_msl_global_0p5deg_grib2.bin<br>
rename *%FD%%T%0000*h_msl_global_0p5deg_grib2.bin
%FD%%T%0000*h_msl_global_0p5deg_grib2.bin<br>
cat %FD%%T%0000*h_msl_global_0p5deg_grib2.bin > mslp.bin<br>
<br>
perl g2ctl.pl -verf mslp.bin > mslp.ctl<br>
gribmap -i mslp.ctl -e<br>
<br>
opengrads -lbxc 'ecmwf_tlak.gs %FD% %T%'<br>
<br>
<br>
wget --ftp-user=wmo --ftp-password=essential
ftp://data-portal.ecmwf.int/%FD%%T%0000/*%FD%%T%0000_an_t_850hPa_global_0p5deg_grib2.bin<br>
rename *%FD%%T%0000_an_t_850hPa_global_0p5deg_grib2.bin
%FD%%T%00000000_C_ECMF_%FD%%T%0000_0h_t_850hPa_global_0p5deg_grib2.bin<br>
wget --ftp-user=wmo --ftp-password=essential
ftp://data-portal.ecmwf.int/%FD%%T%0000/*%FD%%T%0000*h_t_850hPa_global_0p5deg_grib2.bin<br>
rename *%FD%%T%0000*h_t_850hPa_global_0p5deg_grib2.bin
%FD%%T%0000*h_t_850hPa_global_0p5deg_grib2.bin<br>
cat %FD%%T%0000*h_t_850hPa_global_0p5deg_grib2.bin > t850.bin<br>
<br>
perl g2ctl.pl -verf t850.bin > t850.ctl<br>
gribmap -i t850.ctl -e<br>
<br>
opengrads -lbxc 'ecmwf_t850.gs %FD% %T%'<br>
<br>
<br>
wget --ftp-user=wmo --ftp-password=essential
ftp://data-portal.ecmwf.int/%FD%%T%0000/*%FD%%T%0000_an_gh_500hPa_global_0p5deg_grib2.bin<br>
rename *%FD%%T%0000_an_gh_500hPa_global_0p5deg_grib2.bin
%FD%%T%00000000_C_ECMF_%FD%%T%0000_0h_gh_500hPa_global_0p5deg_grib2.bin<br>
wget --ftp-user=wmo --ftp-password=essential
ftp://data-portal.ecmwf.int/%FD%%T%0000/*%FD%%T%0000*h_gh_500hPa_global_0p5deg_grib2.bin<br>
rename *%FD%%T%0000*h_gh_500hPa_global_0p5deg_grib2.bin
%FD%%T%0000*h_gh_500hPa_global_0p5deg_grib2.bin<br>
cat %FD%%T%0000*h_gh_500hPa_global_0p5deg_grib2.bin >
geopot500.bin<br>
<br>
perl g2ctl.pl -verf geopot500.bin > geopot500.ctl<br>
gribmap -i geopot500.ctl -e<br>
<br>
opengrads -lbxc 'ecmwf_geopot.gs %FD% %T%'<br>
<br>
del 20*.bin<br>
<br>
<br>
M.<br>
<br>
<div class="moz-cite-prefix">Dňa 25. 3. 2016 o 17:01 stuart markham
napísal(a):<br>
</div>
<blockquote
cite="mid:AM2PR09MB0339DE381C293BA7AB0E3AA9F5830@AM2PR09MB0339.eurprd09.prod.outlook.com"
type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper"
style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Thanks <span style="font-family: Calibri, Arial, Helvetica,
sans-serif; font-size: 12pt;">Marian,</span></p>
<p><span style="font-family: Calibri, Arial, Helvetica,
sans-serif; font-size: 12pt;"><br>
</span></p>
<p><span style="font-family: Calibri, Arial, Helvetica,
sans-serif; font-size: 12pt;">Adjusting the pres did the
trick thanks. Out of interest i am downloading via this
link <a moz-do-not-send="true"
href="ftp://data-portal.ecmwf.int/">ftp://data-portal.ecmwf.int/</a> 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</span></p>
<p><span style="font-family: Calibri, Arial, Helvetica,
sans-serif; font-size: 12pt;"><br>
</span></p>
<p><span style="font-family: Calibri, Arial, Helvetica,
sans-serif; font-size: 12pt;"></span></p>
<div><br>
</div>
<div><br>
</div>
<div><span style="font-size: 12pt;">:024</span><br>
</div>
<div>wget -O ecmwf_24.grb
<a class="moz-txt-link-freetext" href="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_HHXE50ECMF250000_C_ECMF_20160325000000_24h_gh_500hPa_global_0p5deg_grib2.bin</a>
<a class="moz-txt-link-freetext" href="ftp://wmo:essential@data-portal.ecmwf.int/20160325000000/A_HPXE89ECMF250000_C_ECMF_20160325000000_24h_msl_global_0p5deg_grib2.bin">ftp://wmo:essential@data-portal.ecmwf.int/20160325000000/A_HPXE89ECMF250000_C_ECMF_20160325000000_24h_msl_global_0p5deg_grib2.bin</a></div>
<div>IF %ERRORLEVEL% NEQ 0 @ECHO Waiting for file... &
TIMEOUT /T 180 & GOTO 024</div>
<br>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt"
color="#000000" face="Calibri, sans-serif"><b>From:</b>
<a class="moz-txt-link-abbreviated" href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>
<a class="moz-txt-link-rfc2396E" href="mailto:gradsusr-bounces@gradsusr.org"><gradsusr-bounces@gradsusr.org></a> on behalf of AMS
EmsiWx Bytča <a class="moz-txt-link-rfc2396E" href="mailto:emsiwx@pocasie-bytca.sk"><emsiwx@pocasie-bytca.sk></a><br>
<b>Sent:</b> 25 March 2016 14:06<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] ECMWF Data Issues</font>
<div> </div>
</div>
<div>Stuart,<br>
<br>
This one does not seem correct to me:<br>
<br>
'define PRESmsl = PRESmsl*0.01'<br>
<br>
replace with e.g.<br>
<br>
'define PRES = PRESmsl*0.01'<br>
<br>
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.<br>
<br>
Let me know how it goes, and please show me how you get the
ecmwf grib.<br>
<br>
M.<br>
<br>
<div class="moz-cite-prefix">Dňa 25. 3. 2016 o 13:33 stuart
markham napísal(a):<br>
</div>
<blockquote type="cite">
<div id="divtagdefaultwrapper" style="font-size:12pt;
color:#000000; background-color:#FFFFFF;
font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hi Marian, </p>
<p><br>
</p>
<p>I've attached the full script for you, many thanks
for your help. </p>
<p><br>
</p>
<p>regards, </p>
<p><br>
</p>
<p>Stuart.</p>
<br>
<br>
<div style="color:rgb(0,0,0)">
<hr tabindex="-1" style="display:inline-block;
width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font
style="font-size:11pt" color="#000000"
face="Calibri, sans-serif"><b>From:</b>
<a moz-do-not-send="true"
class="moz-txt-link-abbreviated"
href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>
<a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:gradsusr-bounces@gradsusr.org"><gradsusr-bounces@gradsusr.org></a>
on behalf of AMS EmsiWx Bytča
<a moz-do-not-send="true"
class="moz-txt-link-rfc2396E"
href="mailto:emsiwx@pocasie-bytca.sk"><emsiwx@pocasie-bytca.sk></a><br>
<b>Sent:</b> 25 March 2016 11:56<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] ECMWF Data Issues</font>
<div> </div>
</div>
<div>Hi Stuart,<br>
<br>
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?<br>
<br>
Regards,<br>
<br>
Marian<br>
<br>
<br>
<div class="moz-cite-prefix">Dňa 25. 3. 2016 o 12:11
stuart markham napísal(a):<br>
</div>
<blockquote type="cite">
<div id="divtagdefaultwrapper"
style="font-size:12pt; color:#000000;
background-color:#FFFFFF;
font-family:Calibri,Arial,Helvetica,sans-serif">
<p style="">Hi all,</p>
<p style=""><br>
</p>
<p style="">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. </p>
<p style=""><br>
</p>
<p style=""><br>
</p>
<div>
<div>dset ^ecmwf_grib.grib2</div>
<div>index ^ecmwf_grib.grib2.idx</div>
<div>undef 9.999E+20</div>
<div>title ecmwf_grib.grib2</div>
<div>* produced by g2ctl v0.0.8.8</div>
<div>* command line options: -verf
ecmwf_grib.grib2</div>
<div>* 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)</div>
<div><br>
</div>
<div>dtype grib2</div>
<div>ydef 361 linear -90.000000 0.5</div>
<div>xdef 720 linear 0.000000 0.500000</div>
<div>tdef 11 linear 00Z25mar2016 1dy</div>
<div>zdef 1 linear 1 1</div>
<div>vars 2</div>
<div>HGT500mb 0,100,50000 0,3,5 ** 500 mb
Geopotential Height [gpm]</div>
<div>PRESmsl 0,101 0,3,0 ** mean sea level
Pressure [Pa]</div>
<div>ENDVARS</div>
<div><br>
</div>
</div>
<div><br>
</div>
----------------------------------------------------------------------------------------------------------------------------------------------------------------
<p style=""><br>
</p>
<p style=""><br>
</p>
<p style=""><br>
</p>
<div>'reinit'</div>
<div>'open gfs_grib.ctl'</div>
<div><br>
</div>
<div>* Get info from the descriptor file</div>
<div>'q ctlinfo'</div>
<div>_ctl = result</div>
<div>_undef = getctl(undef)</div>
<div>_tdef = getctl(tdef)</div>
<div>_zdef = getctl(zdef)</div>
<div><br>
</div>
<div>maps = 25</div>
<div><br>
</div>
<div>* Get the Time axis info</div>
<div>tsize = subwrd(_tdef,2) </div>
<div>_t1 = 1 ;* 2nd half of timeseries</div>
<div>_t2 = 45</div>
<div>tsize = _t2 - _t1 + 1</div>
<div>'set t '_t1' '_t2</div>
<div>'q dims'</div>
<div>times = sublin(result,5)</div>
<div>_time1 = subwrd(times,6) </div>
<div>_time2 = subwrd(times,8)</div>
<div>_tdim = _time1' '_time2</div>
<div><br>
</div>
<div>tincr = subwrd(_tdef,5)</div>
<div>_tdef = 'tdef 'tsize' linear '_time1'
'tincr</div>
<div><br>
</div>
<div>huh = subwrd(_tdef,4)</div>
<div><br>
</div>
<div>Best regards,</div>
<div><br>
</div>
<div>Stuart.</div>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre>_______________________________________________
gradsusr mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<br>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre>_______________________________________________
gradsusr mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<br>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
gradsusr mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<br>
</body>
</html>