Find the LCL / 'define' problem

[Wolfgang Linder]

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