xan/xmet
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Fix LCL, LFC (I think)

Browse source
This commit is contained in:
XANTRONIX 2025-03-15 01:32:57 -04:00
parent 4538b3bc84
commit 64fe879d40
1 changed files with 6 additions and 10 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

16
lib/xmet/sounding.py
View file

@ -251,10 +251,9 @@ class SoundingParams():
surface.pressure) surface.pressure)
temp_line = sounding.follow_temp() temp_line = sounding.follow_temp()
dewpoint_line = sounding.follow_dewpoint()
virtual_temp_line = sounding.follow_virtual_temp() virtual_temp_line = sounding.follow_virtual_temp()
dry_adiabat = follow_dry_adiabat(surface.temp, surface.pressure) dry_adiabat = follow_dry_adiabat(tv, surface.pressure)
mixing_ratio_line = follow_saturated_mixing_ratio(surface.dewpoint, mixing_ratio_line = follow_saturated_mixing_ratio(surface.dewpoint,
surface.pressure) surface.pressure)
@ -262,18 +261,15 @@ class SoundingParams():
lcl = dry_adiabat.intersect(mixing_ratio_line, lcl = dry_adiabat.intersect(mixing_ratio_line,
SeriesIntersection.LESSER) SeriesIntersection.LESSER)
tv2 = virtual_temp(temp_line.nearest(lcl[1]), moist_adiabat = follow_moist_adiabat(*lcl)
dewpoint_line.nearest(lcl[1]),
lcl[1])
moist_adiabat = follow_moist_adiabat(tv2-1, lcl[1]) lfc = virtual_temp_line.intersect(moist_adiabat,
SeriesIntersection.LESSER,
lfc = moist_adiabat.intersect(virtual_temp_line, lcl[1]-1)
SeriesIntersection.GREATER)
el = temp_line.intersect(moist_adiabat, el = temp_line.intersect(moist_adiabat,
SeriesIntersection.GREATER, SeriesIntersection.GREATER,
lfc[1]) lfc[1]-1)
self.lcl = lcl self.lcl = lcl
self.lfc = lfc self.lfc = lfc