Compare commits

..

2 commits

View file

@ -260,10 +260,13 @@ class SoundingParams():
moist_adiabat = follow_moist_adiabat(*lcl)
lfc = virtual_temp_line.intersect(moist_adiabat,
SeriesIntersection.LESSER,
lcl[1]-1)
SeriesIntersection.LESSER,
lcl[1]-1)
moist_adiabat = follow_moist_adiabat(*lfc)
tv = virtual_temp_line.nearest(lcl[1])
moist_adiabat = follow_moist_adiabat((lcl[0] + tv) / 2,
lcl[1])
el = temp_line.intersect(moist_adiabat,
SeriesIntersection.GREATER,