Vertical Total IndexIl Vertical Total Index (spesso abbreviato in VTOT) prende in esame le temperature delle topografie assolute di 850 e 500 hPa. In particolare, tiene conto della differenza tra le due temperature per trarre conclusioni sulla capacità dello strato considerato di favorire fenomeni convettivi. La formula, semplicissima, è la seguente: con le temperature espresse in gradi Celsius. Valori di VTOT superiori a 28 indicherebbero un forte potenziale per i temporali. Ottenere con OpenGrADS una cartina che rappresenti questo indice è, di conseguenza, abbastanza elementare. Vediamo come. * Vertical Total Index (VTOT)
* 20150124
* nomefile: vtot.gs
'reset'
'set display color white'
'c'
'rgbset2'
'set t 5'
'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 grid off'
'define VTOT=(tmpprs(lev=850)-tmpprs(lev=500))'
'map europe'
'set gxout shaded'
'set cint 4'
'd VTOT'
'cbarn'
'set gxout contour'
'set cint 4'
'set ccolor 2'
'd VTOT'
'printim vtot_eu.png x800 y600'
'c'
'set map 1 1 10'
'map world'
'set mpdset hires'
'set lon 6 20'
'set lat 36 48'
*'set lat 37 43'
*'set lon 11 20'
'set gxout shaded'
'set cint 2'
'd vtot'
'cbarn'
'set gxout contour'
'set cint 2'
'set ccolor 2'
'd vtot'
'draw title Vertical Total Index\Valid: 'gg_set' 'giorno_v''mese_v''anno_v' at
'ora_v'Z'
'draw string 2.0 0.1 Dati GFS High Resolution - Elaborazione: V.Villasmunta'
'printim vtot_it.png x800 y600' | 
| Aggiungendo opportunamente la riga: 'set cmin 28' possiamo evidenziare solo le aree con valori significativi, | 
|
Altre risorse in rete: |
