Pagina 3 di 8 PrimaPrima 12345 ... UltimaUltima
Risultati da 21 a 30 di 77

Discussione: Domanda su grads

  1. #21
    Burrasca
    Data Registrazione
    22/02/07
    Località
    Molino del Piano(FI)
    Messaggi
    5,587
    Menzionato
    0 Post(s)

    Predefinito Re: Domanda su grads

    Citazione Originariamente Scritto da Fabio Volcano Visualizza Messaggio
    Non va
    Unexpected else. Incorrect if block nesting
    Error occurred scanning line 244
    io mi ricordo di averli fatti plottare i grib2 con grads 1.8 ...

    hai il file wgrib2 vero? e i grib sono del tipo gfs.thhz.pgrb2fHH, vero?
    Stazione meteo: Davis Vantage Vue; Luogo: Molino del Piano, 12 km a NE di Firenze; Alt.: 120 m; Sito web: www.firenzemeteo.it

  2. #22
    Burrasca L'avatar di Fabio Volcano
    Data Registrazione
    22/12/03
    Località
    Triggiano (BA)
    Età
    41
    Messaggi
    6,277
    Menzionato
    2 Post(s)

    Predefinito Re: Domanda su grads

    Citazione Originariamente Scritto da Controcorrente Visualizza Messaggio
    io mi ricordo di averli fatti plottare i grib2 con grads 1.8 ...

    hai il file wgrib2 vero? e i grib sono del tipo gfs.thhz.pgrb2fHH, vero?
    esattamente
    Ma quel messaggio mi esce quando faccio andare lo script

  3. #23
    Burrasca L'avatar di Fabio Volcano
    Data Registrazione
    22/12/03
    Località
    Triggiano (BA)
    Età
    41
    Messaggi
    6,277
    Menzionato
    2 Post(s)

    Predefinito Re: Domanda su grads

    up

  4. #24
    Burrasca
    Data Registrazione
    22/02/07
    Località
    Molino del Piano(FI)
    Messaggi
    5,587
    Menzionato
    0 Post(s)

    Predefinito Re: Domanda su grads

    Citazione Originariamente Scritto da Fabio Volcano Visualizza Messaggio
    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

  5. #25
    train
    Ospite

    Predefinito Re: Domanda su grads

    da un po' di tempo sto trascurando Grads per mille altri impegni; mi spiace non esserti d'aiuto...

  6. #26
    Burrasca L'avatar di Fabio Volcano
    Data Registrazione
    22/12/03
    Località
    Triggiano (BA)
    Età
    41
    Messaggi
    6,277
    Menzionato
    2 Post(s)

    Predefinito Re: Domanda su grads

    Citazione Originariamente Scritto da Controcorrente Visualizza Messaggio
    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
    Grazie per la dritta, me lo vedo questo pacchetto anche perchè so programmare in Fortran (anche se 77)

  7. #27
    Burrasca L'avatar di Fabio Volcano
    Data Registrazione
    22/12/03
    Località
    Triggiano (BA)
    Età
    41
    Messaggi
    6,277
    Menzionato
    2 Post(s)

    Predefinito Re: Domanda su grads

    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

  8. #28
    Burrasca L'avatar di Fabio Volcano
    Data Registrazione
    22/12/03
    Località
    Triggiano (BA)
    Età
    41
    Messaggi
    6,277
    Menzionato
    2 Post(s)

    Predefinito Re: Domanda su grads

    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

  9. #29
    Burrasca
    Data Registrazione
    22/02/07
    Località
    Molino del Piano(FI)
    Messaggi
    5,587
    Menzionato
    0 Post(s)

    Predefinito Re: Domanda su grads

    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

  10. #30
    Burrasca L'avatar di Fabio Volcano
    Data Registrazione
    22/12/03
    Località
    Triggiano (BA)
    Età
    41
    Messaggi
    6,277
    Menzionato
    2 Post(s)

    Predefinito Re: Domanda su grads

    Citazione Originariamente Scritto da Controcorrente Visualizza Messaggio
    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???
    Ok grazie.
    No, a dire il vero non ci ho nemmeno provato.

Segnalibri

Permessi di Scrittura

  • Tu non puoi inviare nuove discussioni
  • Tu non puoi inviare risposte
  • Tu non puoi inviare allegati
  • Tu non puoi modificare i tuoi messaggi
  •