From 7e59d94119dec8896ceed07b77b218526e4d3056 Mon Sep 17 00:00:00 2001
From: Adam Cheng <52572642+adamchengtkc@users.noreply.github.com>
Date: Thu, 16 May 2024 18:42:46 +0200
Subject: [PATCH] fix edge case when beta is too low

---
 warmth/forward_modelling.py | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/warmth/forward_modelling.py b/warmth/forward_modelling.py
index 1d72ea2..b88c041 100644
--- a/warmth/forward_modelling.py
+++ b/warmth/forward_modelling.py
@@ -125,6 +125,9 @@ def _check_beta(self, wd_diff: float, beta_current: float, beta_all: np.ndarray[
         """
         beta_found = False
         beta_all = np.append(beta_all, beta_current)
+        if Wd_diff_all.size > 0:
+            if wd_diff >= Wd_diff_all[0]:
+                wd_diff +=1
         Wd_diff_all = np.append(Wd_diff_all, wd_diff)
         if wd_diff > 0:
             beta_found = True