<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">
<META NAME="GENERATOR" CONTENT="GtkHTML/3.16.1">
</HEAD>
<BODY>
I use this (R is relative humidity):<BR>
Bye bye<BR>
Davide<BR>
<BR>
Rd=287<BR>
Cp=1004<BR>
k1=0.0091379024<BR>
k2=6106.396<BR>
k3=223.1986<BR>
k4=0.0182758048<BR>
k5=-0.37329638<BR>
k6=41.178204<BR>
k7=0.0015945203<BR>
k10=0.1419709<BR>
k11=20.746475<BR>
k12=6672.5713<BR>
k13=0.000036336427<BR>
k15=26.66082<BR>
P0=500<BR>
level=850<BR>
'set lev 'level<BR>
'define epres=exp('k15'-'k1'*t-'k2'/t)*'R'/100'<BR>
'define pom='k15'-log(epres)'<BR>
'define dewtemp=(pom-sqrt(pom*pom-'k3'))/'k4<BR>
'define tcp=dewtemp-(t-dewtemp)*('k5'+'k6'/t+'k7'*t)'<BR>
'define pcp=lev*pow((tcp/t),('Rd'/'Cp'))'<BR>
'define ww=epres/(lev-epres)'<BR>
'define pd=(pcp-'P0'+abs(pcp-'P0'))/(2*abs(pcp-'P0'))'<BR>
if(pd=1)<BR>
'define showalter=t(lev='P0')-t(lev='level')*pow(('P0'/'level'),('Rd'/'Cp'))'<BR>
else<BR>
'define showalter=t(lev='P0')-tcp+(pcp-'P0')*('k10'+'k11'*ww-'k12'*ww/tcp-'k13'*pcp)'<BR>
endif<BR>
<BR>
<BR>
<BR>
<BR>
On Thu, 2008-09-18 at 15:27 +0000, sansaad wrote:<BR>
<BLOCKQUOTE TYPE=CITE>
<BR>
<FONT COLOR="#000000">Hi everybody,</FONT><BR>
<BR>
<FONT COLOR="#000000">Anyone has a script to calculate Showalter Index?</FONT><BR>
<BR>
<FONT COLOR="#000000">Thanks in advance,</FONT><BR>
<BR>
<FONT COLOR="#000000">Sandra</FONT><BR>
</BLOCKQUOTE>
<TABLE CELLSPACING="0" CELLPADDING="0" WIDTH="100%">
<TR>
<TD>
<PRE>
--
Sacchetti Davide
ARPAL UO3 Centro Meteo Idrologico Regione Liguria - Dir. Scientifica
V.le Brigare Partigiane 2 16121 Genova (I)
tel: +39 010 6437535
mail: <A HREF="mailto:davide.sacchetti@arpal.org">davide.sacchetti@arpal.org</A> web: <A HREF="http://www.meteoliguria.it">www.meteoliguria.it</A>
</PRE>
</TD>
</TR>
</TABLE>
</BODY>
</HTML>