Find the LCL / 'define' problem
Wolfgang Linder
linder at METEOTEST.CH
Mon May 9 08:54:00 EDT 2005
Hello Marco,
you can interpolate your variables to pressure levels (zinterp) and
then use this:
**************************************************************
**** Berechnung der thermodynamischen Indices
**** siehe /home/meteo/bin/gp_dwd.pl
**************************************************************
**************************************************************
**** LCL
**************************************************************
*** Temperatur, Druck und Theta auf LCL ***
mb=850
'set lev 'mb
'td=1/(1/tmp-((1.85*log(rh/100))/10000))'
'a = 1/(td-56)'
'b = log(tmp/td)/800'
'tlcl = 56+1/(a+b)'
'thet = tmp*pow(1000/850,0.286)'
'Lcp=(2.501-0.00237*tmp)*1000000/1004.67'
'Tq850=(tmp+Lcp*m)*pow((1000/'mb'),0.286))-273.15'
'plcl = exp(log(1000)-(log(thet/tlcl)/0.286))'
Best regards
Wolfgang
More information about the gradsusr
mailing list