Showalter Index

[Davide Sacchetti]

I use this (R is relative humidity):
Bye bye
Davide

   Rd=287
   Cp=1004
   k1=0.0091379024
   k2=6106.396
   k3=223.1986
   k4=0.0182758048
   k5=-0.37329638
   k6=41.178204
   k7=0.0015945203
   k10=0.1419709
   k11=20.746475
   k12=6672.5713
   k13=0.000036336427
   k15=26.66082
   P0=500
   level=850
   'set lev 'level
   'define epres=exp('k15'-'k1'*t-'k2'/t)*'R'/100'
   'define pom='k15'-log(epres)'
   'define dewtemp=(pom-sqrt(pom*pom-'k3'))/'k4
   'define tcp=dewtemp-(t-dewtemp)*('k5'+'k6'/t+'k7'*t)'
   'define pcp=lev*pow((tcp/t),('Rd'/'Cp'))'
   'define ww=epres/(lev-epres)'
   'define pd=(pcp-'P0'+abs(pcp-'P0'))/(2*abs(pcp-'P0'))'
   if(pd=1)
    'define
showalter=t(lev='P0')-t(lev='level')*pow(('P0'/'level'),('Rd'/'Cp'))'
   else
    'define showalter=t(lev='P0')-tcp
+(pcp-'P0')*('k10'+'k11'*ww-'k12'*ww/tcp-'k13'*pcp)'
   endif




On Thu, 2008-09-18 at 15:27 +0000, sansaad wrote:

>
> Hi everybody,
>
> Anyone has a script to calculate Showalter Index?
>
> Thanks in advance,
>
> Sandra

--
Sacchetti Davide
ARPAL UO3 Centro Meteo Idrologico Regione Liguria - Dir. Scientifica
V.le Brigare Partigiane 2 16121 Genova (I)
tel: +39 010 6437535
mail: davide.sacchetti at arpal.org   web: www.meteoliguria.it
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080919/cba4f5b0/attachment.html