bugfix off of #235

geoschem · Jul 2, 2024 · 06aa4c1 · 06aa4c1
1 parent a94e574
commit 06aa4c1
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/inversion_scripts/invert.py b/src/inversion_scripts/invert.py
@@ -151,10 +151,10 @@ def do_inversion(
             ]
         )
         # Define observational errors (diagonal entries of S_o matrix)
-        obs_error = np.power(obs_error, 2)
+        obs_error = np.power(obs_err, 2)
         gP = s_superO_p**2 / s_superO_1**2
         # scale error variance by gP
-        obs_error = gP * obs_err
+        obs_error = gP * obs_error
 
         # check to make sure obs_err isn't negative, set 1 as default value
         obs_error = [obs if obs > 0 else 1 for obs in obs_error]