Vento geostrofico con OpenGrADSE' un vento inerziale, di equilibrio dinamico, nel quale la forza deviante equilibra la forza di gradiente. Esso corre parallelamente alle isobare, lasciando le basse pressioni alla sua sinistra. Vedi anche: | 'g=9.8' | costante per l'accelerazione di gravità | | 'dtr=3.1415/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 | | '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)' | |
Nello script seguente, calcoliamo il vento geostrofico a 850 hPa (distanza tra isoipse: 30 mgp): *calcolo del vento geostrofico
*nomefile vg.gs
'set display color white'
'c'
'map europe'
'set map 1 1 4'
'set lev 850'
'dtr=3.1415/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)'
'set gxout stream'
'set ccolor 3'
'd ug;vg'
'set gxout contour'
'set cint 30'
'set ccolor 4'
'd hgtprs'
'draw title Vento geostrofico'
'printim vento_geostrofico.gif x800 y600' | 
|
|
