Titolo enigmatico

sicuramente tra voi realizzatori di SWPI ci saranno molti appassionati di Meteo ( voglio sperare ) e quindi voglio condividere quello che sto studiando in questi giorni.
Come certamente saprete i vari modelli meteorologici a larga scala ( tipo il GFS americano ) forniscono dati fruibili da tutti ma con risoluzione , in Europa , di 48 km.
Tale risoluzione è assolutamente insufficiente per previsioni locali dove fattori morfologici e di uso del suolo hanno un impatto rilevante.

Per chi, come gli amanti del volo libero e del volo a vela ( ma non solo ), ha invece bisogno di previsioni più dettagliate l'unica soluzione è quella di far girare un proprio modello metereologico tipo il WFR su una risoluzione maggiore ed inizializzato con i dati GFS.

Un forte contributo all'uso di questi modelli ( da far girare in casa ) è stato quello di un americano ( Jack Glendening ) che nel 2008 ha sviluppato e pubblicato un pacchetto software facilmente installabile su macchine linux chiamato RASP ( nulla a che vedere con in RASPBerry ) e che non è altro che una serie di script che automatizzano dei run del WRF ( con qualche aggiunta legata a dei parametri di volo )

Gli operatori che fanno girare i propri modelli RASP locali sono molti oramai ( vedete la cartina nel sito ) e nelle ultime 2 settimane mi sono divertito ad installare una macchina RASP e a produrre delle previsione per il Cucco con risoluzione di 2 km ( sono quelle che vedete nel menu meteo alla voce windgram o blipmap )

La cosa è relativamente semplice ma ovviamente richiede una macchina dedicata abbastanza potente ( il WRF è un pò pesantino ) accesa 24/24 con tutte le implicazioni connesse. Impossibile farlo girare sul Raspberry

Ieri mi sono poi imbattuto in questo fantastico progetto di un francese : openmeteodata.org/
Invece di far girare tanti modelli farne girare 1 solo per l'Europa e distribuire l'output del modello con licenza Open Data !!!!

In questo momento stò studiando come accedere a quei dati mediante OPeNDAP e produrre degli output grafici con il Raspberry. Potrebbe essere un plugin interessante per SWPI che a quel punto fornirebbe anche delle previsioni iper-locali.

TT

Se la segui attentamente funzia tutto . Io ho usato un libuntu 13.4 su macchina virtuale Virtual Box

Allegato un script con i passaggi che ho seguito .. ti può aiutare. E' anche numerato secondo la guida di Jack.

OMD è semplicissimo. Quando fai girare il WRF ottieni dei file netCDF ( si chiamano wrfout.....nc ) con tutti i parametri di output del modello. Alcuni sono direttamente i parametri che trovi nei bipmap ( come lo spessore del BL ) altri invece ti permettono di calcolare gli altri parametri RASP ( come la quota massima in termica )
La parte pesante del lavoro è il modello mentre il post-processing può girare tranquillamente anche sul RaspBerry.
OMD fornisce i file output del modello ed ( per ora ) anche i file post-procesati ) con risoluzione di 2.5 km ed intervallo temporale ogni ora. L'accesso ai dati poi è semplificato dal protocollo OpenDAP che ti permette di scaricare solo la parte di territorio che ti interessa.


Per NCL è un pò ostico ma ho dovuto sbatterci il naso per i windgram .La versione di TJ che scarichi ( la 1 nel mio sito ) non era un gran che.

Ora con OMD pensavo di usare qualcosa di più semplice ( mtplotlib e Basemap Matplotlib Toolkit ) che sto proprio ora compilando per raspberry.

Ottenere i blipmap dai dati OMD è semplicissimo e sono 10 istruzioni in python .. ora per i windgram sarà un pò più complicato.

Si installare e giocare con RASP è il modo migliore di capirci un pò qualcosa. Sarebbe poi interessante vedere se è possibile compilare NCL e il suo wrapper in python sul Raspberry ... www.ncl.ucar.edu/

RaspGM per ora lo puoi saltare . E' praticamente la parte di Paul Scorer per visualizzare le blipmap su Google Map. rasp.inn.leedsmet.ac.uk/

I GEOG invece li devi scaricare altrimenti non riuscirai a fare nuovi domini .. vedi bene .. nella guida c'è ma forse nella parte di creazione domini ... per far girare l'esempio SW_SOUTHAFFRICA non servono in quanto il dominio è già bello che pronto.
Sono praticamente i dati territoriali : DTM,Albedo, Riflettanza ,uso del suolo .... Il GUI di creazione domini legge quei file e ne fà un subset locale ( quando lanci la localizzazione .... ma questo lo vedrai dopo ) per ora cerca di far girare l'esempio sia in single stage che in doppio stage.

Consiglio : Installa il tutto nella tua home sotto la cartella DRJACK .. è tutto più semplice.

vedi il file log è quello che finisce in .printout

dopo aver lanciato : ./run.rasp SW_SOUTHAFRICA da un'altra finestra fai un :

tail -f rasp.SW_SOUTHAFRICA.printout

Controlla il nome non mi ricordo bene.

Se esce subito vuol dire che c'è qualcosa che non non và nello scarico dei dati GFS via ftp.

se non scrive nulla in printout allora controlla il file che.finisce in stderr

Trovi un pò di documentazione qui :

openmeteodata.org/index.php/projects/rasp-france/data/

ho fatto dei test la settimana scorsa e mi sono un pò bloccato con le API php quando ho notato che le richieste per una determinata posizioni restituivano i dati per una posizione diversa !!!

per esempio la richiesta ( per Temperatura e pioggia ) :

ti da i dati per la la location :

Non ho indagato oltre perché sono passato poi ad analizzare l'accesso ai dati tramite OpenDAP.

Ho già una prima versione del sw ma prima di pubblicarla devo separare la configurazione dal codice ... i risultati li vedi nel menu Blipmap a sinistra.

I plottaggi con il Raspberry e MatPlotLib sono abbastanza pesanti quindi la mia idea è quella di fare 4 mappe sintetiche che raccolgano un po tutti i parametri principali per il volo, Le prime 2 le vedi già : Convergenza e Top BL.

Per cosa ti interessa la cosa ??? Volo o Meteo ???

TT

Si quindi puoi usare le API php ( vedi il riferimento che ti ho postato ) ma bisogna analizzare il problema di localizzazione che ti dicevo.

Questo è un esempietto in python su come leggere i dati :

TT

Pero tu hai richiesto :

41.7004957,12.6501506

e lui ha risposto con :

[location] => Albano Laziale, Rome, IT
[lat] => 41.697448730469
[lon] => 12.645589828491

che non è la stessa cosa !!!!

Prova a vedere la distanza tra quei 2 punti . Inutile avere una risoluzione di 2.5km se poi risponde con una posizione oltre un margine di 2.5km !!! Fammi sapere

si hai ragione .. ho approfondito. Anche nel mio caso la distanza è solo di 1.3 km

Mi sembra che il modello a 12 km dovrebbe arrivare alle 72 ore mentre in quello a 2.5 , malgrado la documentazione, arriva a 48 ma le previsioni del giorno successivo non escono di norma prima delle 2-3 del pomeriggio.

Tieni presente poi che l'orario del file dovrebbe essere GMT

Una prima versione del sw per l'uso di OMD è disponibile nel mio repository. Praticamente è una prima versione che uso per plottare queste mappe : www.vololiberomontecucco.it/vlmc/index.p...w=wrapper&Itemid=265

Se mi date due info veloci su come installarle e come usarle, le provo anche io volentieri !

La fase di test è oramai conclusa e le blitmap a 2.5 km funzionano oramai da più di una settimana su www.vololiberomontecucco.it/vlmc/index.p...w=wrapper&Itemid=265

L'uso è molto semplice e simile a quello di swpi. Al primo avvio del programma python ( omdpi.py ) viene generato un file di configurazione che va editato ( nessun configuratone web con i propri parametri.

La generazione delle mappe è particolarmente gravosa con il Raspberry quindi l'area deve essere ristretta magari ad un comune o provincia. Nel mio caso è la nostra area di gara che è di circa 100x100km.
In questo caso il raspberry impiega circa 5 minuti a generare una mappa.

Dove trovo il sorgente da scaricare tony?

Il mio solito repository

github.com/ToninoTarsi/OpenMeteoData-PI