Home ] Su ]
 

Vento ageostrofico con OpenGrADS

 

Vedi anche:

'g=9.8'costante per l'accelerazione di gravitÓ
'dtr=pi/180'per convertire i gradi in radianti
'a=6.37122e6' 
'omega=7.2921e-5' 
'define f=2*omega*sin(lat*dtr)'computo della rotazione locale (parametro di Coriolis).
'dy=cdiff(lat,y)*dtr*a' 
'dx=cdiff(lon,x)*dtr*a*cos(lat*dtr)' 
'dhgtx=cdiff(hgtprs,x)' 
'dhgty=cdiff(hgtprs,y)' 
'define ug=-1*(g/f)*(dhgty/dy)' 
'define vg=(g/f)*(dhgtx/dx)' 
'define ua=ugrdprs-ug' 
'define va=vgrdprs-vg' 

Nello script seguente, calcoliamo il vento ageostrofico a 850 hPa (distanza tra isoipse: 60 mgp) e lo confrontiamo con il vento reale:

*calcolo del vento ageostrofico
*nomefile va.gs

'set display color white'
'c'
'rgbset2'
'map europe'
'set map 1 1 4'
'set lev 850'

'dtr=pi/180'
'g=9.8'
'a=6.37122e6'
'omega=7.2921e-5'

'define f=2*omega*sin(lat*dtr)'
'dy=cdiff(lat,y)*dtr*a'
'dx=cdiff(lon,x)*dtr*a*cos(lat*dtr)'
'dhgtx=cdiff(hgtprs,x)'
'dhgty=cdiff(hgtprs,y)'
'define ug=-1*(g/f)*(dhgty/dy)'
'define vg=(g/f)*(dhgtx/dx)'
'define ua=ugrdprs-ug'
'define va=vgrdprs-vg'

'set gxout barb'
'set ccolor 2'
'd skip(ua,5,5);va'

'set gxout barb'
'set ccolor 3'
'd skip(ugrdprs,5,5);vgrdprs'

'set gxout contour'
'set cint 60'
'set ccolor 4'
'd hgtprs'

'draw title Vento ageostrofico (rosso), vento reale (verde) [m/s]'

'printim vento_ageostrofico.gif x800 y600'


Questa pagina Ŕ stata realizzata da Vittorio Villasmunta

Ultimo aggiornamento: 14/05/16