From 75a3b0ed1b31cacd38609b95052b755029799626 Mon Sep 17 00:00:00 2001 From: Ansgar Wehrhahn Date: Thu, 3 Dec 2020 15:31:48 +0100 Subject: [PATCH] Bugfix in Atmosphere interpolation --- src/pysme/atmosphere/interpolation.py | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/src/pysme/atmosphere/interpolation.py b/src/pysme/atmosphere/interpolation.py index c5ad38d3..1a03efd8 100644 --- a/src/pysme/atmosphere/interpolation.py +++ b/src/pysme/atmosphere/interpolation.py @@ -784,8 +784,8 @@ def interpolate_corner_models( p = monh for (i, j) in np.ndindex(2, 2): - m0 = atmo_grid[icor[0, i, j]] - m1 = atmo_grid[icor[1, i, j]] + m0 = atmo_grid[icor[0, j, i]] + m1 = atmo_grid[icor[1, j, i]] p0 = getattr(m0, param) p1 = getattr(m1, param) pfrac = (p - p0) / (p1 - p0) if p0 != p1 else 0 @@ -796,11 +796,11 @@ def interpolate_corner_models( # Interpolate 4 models at the desired [M/H] to create 2 models at desired # [M/H] and log(g). atmo2 = [None, None] - param = "teff" - p = teff + param = "logg" + p = logg for (k,) in np.ndindex(2): - m0 = atmo[0, k] - m1 = atmo[1, k] + m0 = atmo[k, 0] + m1 = atmo[k, 1] p0 = getattr(m0, param) p1 = getattr(m1, param) pfrac = (p - p0) / (p1 - p0) if p0 != p1 else 0