Stazione meteo: Davis Vantage Vue; Luogo: Molino del Piano, 12 km a NE di Firenze; Alt.: 120 m; Sito web: www.firenzemeteo.it
up
nessuna novità?
io ho riprovato con grads a fare girare i grib2 e non mi riesce più
ho visto che qui http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2/ esistono delle utility per convertire i grib 2 in grib 1, ma non so da dove partite ... è roba in fortran
cmq è il pacchetto cnvgrib
Stazione meteo: Davis Vantage Vue; Luogo: Molino del Piano, 12 km a NE di Firenze; Alt.: 120 m; Sito web: www.firenzemeteo.it
da un po' di tempo sto trascurando Grads per mille altri impegni; mi spiace non esserti d'aiuto...
Domanda: ho creato il seguente file eseguibile.
reset
open d:\Fabio\Meteo\GRADS\HR\gfs_00z.ctl
set display color white
set mpdset hires
set poli on
set gxout linefill
set csmooth linear
set lat 39 42
set lon 14 19
set t 1
set lev 850
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\850-1.png
c
set lev 900
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\900-1.png
c
set lev 950
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\950-1.png
c
set t 2
set lev 850
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\850-2.png
c
set lev 900
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\900-2.png
c
set lev 950
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\950-2.png
Ora, non capisco perchè, dopo ogni comando c si blocca, o meglio attende che io prema il tasto Invio. Come mai non prosegue?
Ultima modifica di Fabio Volcano; 04/01/2009 alle 11:38
Se invece dopo ogni comando printim faccio il reset e il settaggio completo, tutto va bene. E' così che devo procedere? Ma il file in questo modo viene immenso
reset
open d:\Fabio\Meteo\GRADS\HR\gfs_00z.ctl
set display color white
set mpdset hires
set poli on
set gxout linefill
set csmooth linear
set lat 39 42
set lon 14 19
set t 1
set lev 850
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\850-1.png
reset
open d:\Fabio\Meteo\GRADS\HR\gfs_00z.ctl
set display color white
set mpdset hires
set poli on
set gxout linefill
set csmooth linear
set lat 39 42
set lon 14 19
set t 1
set lev 900
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\900-1.png
reset
open d:\Fabio\Meteo\GRADS\HR\gfs_00z.ctl
set display color white
set mpdset hires
set poli on
set gxout linefill
set csmooth linear
set lat 39 42
set lon 14 19
set t 1
set lev 950
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\950-1.png
reset
open d:\Fabio\Meteo\GRADS\HR\gfs_00z.ctl
set display color white
set mpdset hires
set poli on
set gxout linefill
set csmooth linear
set lat 39 42
set lon 14 19
set t 2
set lev 850
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\850-2.png
reset
open d:\Fabio\Meteo\GRADS\HR\gfs_00z.ctl
set display color white
set mpdset hires
set poli on
set gxout linefill
set csmooth linear
set lat 39 42
set lon 14 19
set t 2
set lev 900
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\900-2.png
reset
open d:\Fabio\Meteo\GRADS\HR\gfs_00z.ctl
set display color white
set mpdset hires
set poli on
set gxout linefill
set csmooth linear
set lat 39 42
set lon 14 19
set t 2
set lev 950
d tmp-273.15
printim d:\Fabio\Meteo\GRADS\HR\950-2.png
1) Perché hai creato un eseguibile?? Devi creare uno script con estensione.gs che lanci da dos in questo modo:
grads -lxc "run script_grads.gs" (dove script_grads.gs è lo script che hai creato)
I comandi dello script devono essere all'interno di 2 apostrofi (tipo: 'set lev 600' ). Ecco un esempio di script da pochi kB:
'open gfs.t06z.master.grbf180.ctl'
'query time'
curtime=subwrd(result,3)
ora_v=substr(curtime,1,2)
giorno_v=substr(curtime,4,2)
mese_v=substr(curtime,6,3)
anno_v=substr(curtime,9,4)
gg_set=subwrd(result,6)
'set display color white'
'set mpdraw off'
'set grads off'
'set lat 42.3 44.5'
'set lon 9.6 12.4'
'run script_tmp2m.gs'
'set line 15'
'draw title Init: Sun 04Jan2009 at 06Z Valid: 'gg_set' 'giorno_v''mese_v''anno_v' at 'ora_v'Z \2m Temperature (C)'
'draw string 2.7 0.3 Dati GFS High Resolution - Elaborazione: www.firenzemeteo.it'
'printim tuscany_2mT_' 180'.png'
'c'
'set display color white'
'set mpdraw off'
'set grads off'
'set lat 42.3 44.5'
'set lon 9.6 12.4'
'run script_0level.gs'
'set line 15'
'draw title Init: Sun 04Jan2009 at 06Z Valid: 'gg_set' 'giorno_v''mese_v''anno_v' at 'ora_v'Z \0C Isotherm Level (m) and 3-Hours Precipitation (mm)'
'draw string 2.7 0.3 Dati GFS High Resolution - Elaborazione: www.firenzemeteo.it'
'printim tuscany_0C_level_' 180'.png'
Come vedi il file grib viene aperto una sola volta ed è possibile lanciare altri script .gs
2) Ma alla fine ce l'hai fatta a convertire i grib2 in grib1???
Stazione meteo: Davis Vantage Vue; Luogo: Molino del Piano, 12 km a NE di Firenze; Alt.: 120 m; Sito web: www.firenzemeteo.it
Segnalibri