[ Home ] [ Su ]
 

Realizzare diagrammi con OpenGrADS

La rappresentazione delle variabili meteorologiche mediante diagrammi possiede l'indubbio vantaggio di rendere più facilmente visibile il comportamento delle medesime variabili in relazione, ad esempio, all'altezza rispetto alla superficie terrestre.

L'esempio che segue esamina il comportamento della temperatura dell'aria (in °C) e della temperatura del punto di rugiada (in °C) con l'altezza.

La temperatura del punto di rugiada viene calcolata attraverso la formula:

 

td = tc-((14.55+0.114*tc)*(1-0.01*rhprs) + pow((2.5+0.007*tc)*(1-0.01*rhprs),3) + (15.9+0.117*tc)*pow((1-0.01*rhprs),14))

dove:

  • tc è la temperatura dell'aria in °C

  • rhprs è la variabile che contiene i dati di umidità relativa

 

*nomefile diagramma.gs
'set display color white'
'c'
'set t 1'

'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 grads off'
'set lat 40.76'
'set lon 16.96'
'set lev 1000 100'
'set zlog on'
'tc = tmpprs-273.16'
'td = tc-((14.55+0.114*tc)*(1-0.01*rhprs) + pow((2.5+0.007*tc)*(1-0.01*rhprs),3) + (15.9+0.117*tc)*pow((1-0.01*rhprs),14))'
'set ccolor 2'
'd tc'
'set ccolor 3'
'd td'
'draw title Curva di stato e rugiada\Valid: 'gg_set' 'giorno_v''mese_v''anno_v' at 'ora_v'Z'
'draw string 1.2 0.3 Dati GFS High Resolution - Elaborazione: V.Villasmunta'

'printim diagramma.png x800 y600'
 

Nell'esempio che segue, tracciamo in un diagramma a parte, la curva dell'umidità relativa rispetto alla quota (espressa sempre in pressioni in scala logaritmica):

*nomefile diagramma.gs
'set display color white'
'c'
'set t 1'

'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 grads off'
'set lat 40.76'
'set lon 16.96'
'set lev 1000 100'
'set zlog on'
'tc = tmpprs-273.16'
'td = tc-((14.55+0.114*tc)*(1-0.01*rhprs) + pow((2.5+0.007*tc)*(1-0.01*rhprs),3) + (15.9+0.117*tc)*pow((1-0.01*rhprs),14))'
'set ccolor 2'
'd tc'
'set ccolor 3'
'd td'
'draw title Curva di stato e rugiada\Valid: 'gg_set' 'giorno_v''mese_v''anno_v' at 'ora_v'Z'
'draw string 1.2 0.3 Dati GFS High Resolution - Elaborazione: V.Villasmunta'

'printim diagramma.png x800 y600'
'c'
'set grads off'
'set lat 40.76'
'set lon 16.96'
'set lev 1000 100'
'set zlog on'
'set ccolor 3'
'd rhprs'
'draw title Ur%\Valid: 'gg_set' 'giorno_v''mese_v''anno_v' at 'ora_v'Z'
'draw string 1.2 0.3 Dati GFS High Resolution - Elaborazione: V.Villasmunta'
'printim diagramma_ur.png x800 y600'
 
Questa pagina è stata realizzata da Vittorio Villasmunta

Ultimo aggiornamento: 14/05/16