lib/xmet/skew_t.py

@@ -2,11 +2,14 @@ import math
 import cairo
 
 from typing    import Callable
+from itertools import product
 
-from xmet.sounding import Sounding
-from xmet.thermo   import loft_parcel, moist_lapse_rate, \
+from xmet.list     import nearest
+from xmet.sounding import Sounding, SoundingSample
+from xmet.thermo   import pressure_height, loft_parcel, moist_lapse_rate, \
                           LAPSE_RATE_DRY, PRESSURE_MAX, PRESSURE_MIN, \
-                          PRESSURE_STEP
+                          PRESSURE_STEP, mixing_ratio, saturated_mixing_ratio, \
+                          mixing_ratio_temp
 
 PRESSURE_LOG_MAX   = math.log(PRESSURE_MAX)
 PRESSURE_LOG_MIN   = math.log(PRESSURE_MIN)

lib/xmet/thermo.py

@@ -106,6 +106,8 @@ def loft_parcel(start_temp: float,
     dewpoint, pressure value, will be dispatched to obtain the lapse rate
     for each height.
     """
+    sx_last      = None
+    sy_last      = None
     last_height  = None
 
     temp     = start_temp