def calcolaventogeostrofico(lat,GradientePressione,R):
omega=2*math.pi/(24*3600)#rotazione terrestre
f=2*omega*math.sin(lat)#parametro di coriolis
densitaaria=1.2#kgm/mc vale per T=20gradi
R=R*1853.0#curvatura isobare espressa in m
GradientePressione=GradientePressione*100.0/1853.0#Gradiente espresso in Pa/m
segno=GradientePressione/math.copysign(GradientePressione,1)
a=segno/R
b=-f
c=math.copysign(GradientePressione,1)/densitaaria
discriminante=(b**2-4*a*c)
if discriminante >0:
tws=(-b-discriminante**0.5)/(2*a)
tws=tws*3600.0/1853.0#converte tws in knts
else:
tws=0.0#caso del centro di alta pressione
return tws
评论列表
文章目录