Vento ageostrofico con OpenGrADSVedi 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' | 
|
|
