-
Notifications
You must be signed in to change notification settings - Fork 1
alt0
Tentativo di sintetizzare una strategia implementativa per l'altezza dello zero termico (#5.
Per COSMO è un semplice rendering di contouring di:
arki-query --data 'reftime:=today 0:00;product:GRIB1,80,201,84;level:GRIB1,4' http://arkiope.metarpa:8090/dataset/cosmo_5M_ita > hzero_cosmo.grib
ECMWF usa l'altezza dal suolo invece che dal livello del mare, e quindi ai dati di ECMWF va sommata l'altezza dall'orografia.
Altezza dello zero termico dal suolo:
arki-query --data 'reftime:=today 0:00;product:GRIB1,98,228,24;level:GRIB1,1' http://arkiope.metarpa:8090/dataset/ifs_ita010 > hzero_ecmwf.grib
Orografia:
arki-query --data 'reftime:=today 0:00;product:GRIB1,98,128,129;level:GRIB1,1' http://arkiope.metarpa:8090/dataset/ifs_ita010 > orografia_ita.grib
Le due variabili non sono presenti in tutti i dataset (ad esempio in ifs_atl025 non è presente).
Si è optato per plottare l'altezza dello zero termico solo se il pacco dati contiene sia l'altezza da terra che l'orografia.
Si potrebbe mettere in cima alla ricetta un commento che spiega come ottenere l'orografia e come renderla disponibile ad arkimaps.
Rimane il problema di dove metterli: questo introduce il bisogno di una directory di dati locale, la cui posizione andrebbe configurata in qualche modo.
Si può usare un elenco di directory hardcoded in arkimaps, magari con override dall'environment, o una directory configurabile in un file.
Nella seconda ipotesi, questo introduce l'uso di un file di configurazione di arkimaps, che implica decidere in quali luoghi tenerlo e cercarlo.
Per scaricare le orografie da arkimet, c'è il problema di dove trovarle: servirebbe una query del tipo "dei risultati che hai, raggruppa per reftime e dammi l'ultimo", in modo da poter fare query tipo "dammi l'ultima orografia disponibile prima del giorno X". Se serve possiamo aprire un ticket in arkimet per implementarla.
TODO: passare questa sezione da una raccolta di possibilità a una strategia implementativa specifica.
Se abbiamo il permesso di redistribuire le orografie, possiamo includerle in arkimaps per una selezione di aree e risoluzioni, però questo rende il prodotto calcolabile solo per le aree e risoluzioni per le quali arkimaps fornisce orografie, e gonfia parecchio la dimensione di arkimaps (solo per il dataset ita le orografie 2Mb).
TODO: Scartabile
Una volta che abbiamo una serie di orografie disponibili, serve scegliere quale si applica ai dati. In teoria, possiamo confrontare i metadati di griglia di GRIB, e selezionare tra le orografie quella che ha i metadati di griglia identici.
TODO: Possibile strategia futura in caso di necessità