Home ] Su ]
 

Fog Stability Index

Nella pagina "Indice di stabilità per la formazione della nebbia da irraggiamento", abbiamo trattato alcuni aspetti intorno a questo indice. In questa sezione, relativa ad OpenGrADS, vedremo come ottenere sia un grafico dell'indice su una data località (ovviamente con tutte le limitazioni imposte dal modello numerico), sia una carta riferita ad un'area geografica.

I valori che l'indice può assumere sono illustrati nella seguente tabella:

VALORE DELL'INDICECOLOREPOSSIBILITA' CHE SI FORMI LA NEBBIA DA IRRAGGIAMENTO
fogsi< 31rossoforte
31<=  fogsi < 55giallomoderato
fogsi >= 55 verdedebole

 

* vittorio villasmunta
* fog stability index LIBV
* RH media 975-950
* nomefile: FOGSI_LIBV.gs

'reset'
'set display color white'

'c'
'rgbset2'

'set grads off'

'set lat 40.76'
'set lon 16.96'

'set t 1 12'
'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)

'define t2m=tmp2m-273.16'
'define dewp2m = tmp2m-((14.55+0.114*tmp2m)*(1-0.01*rh2m)+pow((2.5+0.007*tmp2m)
*(1-0.01*rh2m),3)+(15.9+0.117*tmp2m)*pow((1-0.01*rh2m),14))'
'define d2m=dewp2m-273.16'
'define t850=tmpprs(lev=850)-273.16'
'define u850=ugrdprs(lev=850)*1.9438'
'define v850=vgrdprs(lev=850)*1.9438'
'define w850=mag(u850,v850)'

'define FSI=4*t2m-2*(t850+d2m)+w850'

'set gxout contour'
'set ccolor 1'
'set csmooth on'
'd FSI'

* RH media
'define RH=ave(rhprs,lev=975,lev=950)'

'define dewp = tmp2m-((14.55+0.114*tmp2m)*(1-0.01*rh)+pow((2.5+0.007*tmp2m)*(1-
0.01*rh),3)+(15.9+0.117*tmp2m)*pow((1-0.01*rh),14))'
'define dm2m=dewp2m-273.16'
'define FSIm=4*t2m-2*(t850+dm2m)+w850'

'set gxout contour'
'set ccolor 3'
'set csmooth on'
'd FSIm'

'draw title Fog Stability Index\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 fogsi_libv.png x1100 y850'

 

* Vittorio Villasmunta
* Fog Stability Index
* nomefile: FOGSI.gs
* 20111030, 20150118
 

'reset'
'set display color white'
'c'
palette()
index=1

'set grads off'
'set grid off'
'map world'
'set mpdset hires'
'set lat 38.5 42'
'set lon 15 19'

while(index<=33)
'c'

'set t 'index
'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)

'define t2m=tmp2m-273.16'
'define dewp2m = tmp2m-((14.55+0.114*tmp2m)*(1-0.01*rh2m)+pow((2.5+0.007*tmp2m)*(1-0.01*rh2m),3)+(15.9+0.117*tmp2m)*pow((1-0.01*rh2m),14))'
'define d2m=dewp2m-273.16'
'define t850=tmpprs(lev=850)-273.16'
'define u850=ugrdprs(lev=850)*1.9438 '
'define v850=vgrdprs(lev=850)*1.9438'
'define w850=mag(u850,v850)'
'define FSI=4*t2m-2*(t850+d2m)+w850'

'set clevs 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110'
'set ccols 0 16 17 18 19 20 21 22 23 24 25 26 31 32 33 34 35 36 37 38 39 40'
'set gxout shaded'
'set csmooth on'
'd FSI'
'cbarn'

'set gxout contour'
'set ccolor 2'
'set csmooth on'
'd FSI'

'draw title Fog Stability Index\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'
'draw mark 3 5.76 5.26 .1'
'draw string 5.89 5.26 LIBV'

'printim fogsi_'index'.png x800 y600'

index=index+1
endwhile

function palette()
'set rgb 16 0 0 255'
'set rgb 17 20 20 255'
'set rgb 18 50 50 255'
'set rgb 19 75 75 255'
'set rgb 20 100 100 255'
'set rgb 21 125 125 255'
'set rgb 22 150 150 255'
'set rgb 23 175 175 255'
'set rgb 24 200 200 255'
'set rgb 25 225 225 255'
'set rgb 26 255 255 255'

* These are the RED shades
'set rgb 31 255 225 225'
'set rgb 32 255 200 200'
'set rgb 33 255 175 175'
'set rgb 34 255 150 150'
'set rgb 35 255 125 125'
'set rgb 36 255 100 100'
'set rgb 37 255 75 75'
'set rgb 38 255 50 50'
'set rgb 39 255 25 25'
'set rgb 40 255 0 0'

return

E' possibile attribuire colori che si riferiscono direttamente alla tabella di cui sopra (rosso, giallo e verde):

* Vittorio Villasmunta
* Fog Stability Index
* nomefile: FOGSI3c.gs
* 20111030, 20150118
* If fogsi < 31 Then red
* If fogsi >= 31 And fogsi < 55 Then yellow
* If fogsi >= 55 Then green

'reset'
'set display color white'
'c'
palette()
index=1

'set grads off'
'set grid off'
'map world'
'set mpdset hires'
'set lat 38.5 42'
'set lon 15 19'

while(index<=5)
'c'

'set t 'index
'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)

'define t2m=tmp2m-273.16'
'define dewp2m = tmp2m-((14.55+0.114*tmp2m)*(1-0.01*rh2m)+pow((2.5+0.007*tmp2m)
*(1-0.01*rh2m),3)+(15.9+0.117*tmp2m)*pow((1-0.01*rh2m),14))'
'define d2m=dewp2m-273.16'
'define t850=tmpprs(lev=850)-273.16'
'define u850=ugrdprs(lev=850)*1.9438 '
'define v850=vgrdprs(lev=850)*1.9438'
'define w850=mag(u850,v850)'
'define FSI=4*t2m-2*(t850+d2m)+w850'

'set clevs 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110'
'set ccols 0 16 17 18 19 20 21 22 23 24 25 26 31 32 33 34 35 36 37 38 39 40'
'set gxout shaded'
'set csmooth on'
'd FSI'
'cbarn'

'set gxout contour'
'set ccolor 2'
'set csmooth on'
'd FSI'

'draw title Fog Stability Index\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'
'draw mark 3 5.76 5.26 .1'
'draw string 5.89 5.26 LIBV'

'printim fogsi_'index'.png x800 y600'

index=index+1
endwhile

function palette()
'set rgb 16 255 0 0'
'set rgb 17 255 50 0'
'set rgb 18 255 75 0'
'set rgb 19 255 100 0'
'set rgb 20 255 150 0'
'set rgb 21 255 255 0'
'set rgb 22 245 255 0'
'set rgb 23 235 255 0'
'set rgb 24 225 255 0'
'set rgb 25 210 255 0'
'set rgb 26 200 255 0'
'set rgb 31 190 255 0'
'set rgb 32 180 255 0'
'set rgb 33 170 255 0'
'set rgb 34 160 255 0'
'set rgb 35 150 255 0'
'set rgb 36 140 255 0'
'set rgb 37 130 255 0'
'set rgb 38 120 255 0'
'set rgb 39 110 255 0'
'set rgb 40 100 255 0'

return

 


Questa pagina è stata realizzata da Vittorio Villasmunta

Ultimo aggiornamento: 14/05/16