From bf66646e6aad48743747872db6f2590998d856f8 Mon Sep 17 00:00:00 2001
From: yoctoyotta1024 <92624833+yoctoyotta1024@users.noreply.github.com>
Date: Mon, 28 Jul 2025 11:23:45 +0200
Subject: [PATCH] fix: bug on interpolation of qv and _th when apprx_drhod_dz
 is false

---
 .../PyMPDATA_examples/Shipway_and_Hill_2012/settings.py   | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/examples/PyMPDATA_examples/Shipway_and_Hill_2012/settings.py b/examples/PyMPDATA_examples/Shipway_and_Hill_2012/settings.py
index 1490cfcb..1f86ad92 100644
--- a/examples/PyMPDATA_examples/Shipway_and_Hill_2012/settings.py
+++ b/examples/PyMPDATA_examples/Shipway_and_Hill_2012/settings.py
@@ -37,8 +37,12 @@ def __init__(
         self.z_max = z_max
         self.t_max = t_max
 
-        self.qv = interp1d((0, 740, 3260), (0.015, 0.0138, 0.0024))
-        self._th = interp1d((0, 740, 3260), (297.9, 297.9, 312.66))
+        self.qv = interp1d(
+            (0, 740, 3260), (0.015, 0.0138, 0.0024), fill_value="extrapolate"
+        )
+        self._th = interp1d(
+            (0, 740, 3260), (297.9, 297.9, 312.66), fill_value="extrapolate"
+        )
 
         # note: not in the paper,
         # https://github.com/BShipway/KiD/tree/master/src/physconst.f90#L43