From b5b5d4efbab080522273bd23d9941c5c4ab133ae Mon Sep 17 00:00:00 2001
From: XANTRONIX Industrial <xan@xantronix.com>
Date: Wed, 19 Mar 2025 11:35:06 -0400
Subject: [PATCH] Declare LAPSE_RATE_MAGIC_FACTOR
---
lib/xmet/thermo.py | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/lib/xmet/thermo.py b/lib/xmet/thermo.py
index e14dfb6..c4fe55f 100644
--- a/lib/xmet/thermo.py
+++ b/lib/xmet/thermo.py
@@ -5,7 +5,8 @@ from typing import Callable
from xmet.series import Series
from xmet.units import rad, deg, kelvin, celsius
-LAPSE_RATE_DRY = 9.8 / 1000 # degrees C per km
+LAPSE_RATE_DRY = 9.8 / 1000 # degrees C per km
+LAPSE_RATE_MAGIC_FACTOR = 1.09 # adjustment for accuracy
PRESSURE_MAX = 1050 # millibar
PRESSURE_MSL = 1013.25
@@ -130,7 +131,8 @@ def loft_parcel(start_temp: float,
except OverflowError:
break
- temp = lapse(temp, height - last_height, 1.09*rate)
+ temp = lapse(temp, height - last_height,
+ LAPSE_RATE_MAGIC_FACTOR*rate)
yield temp, pressure