#9 indexing bugfix to ensure correct PWV for each beam

deshima-dev · Mar 26, 2021 · bcb695e · bcb695e
1 parent f03f20e
commit bcb695e
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/pyproject.toml b/pyproject.toml
@@ -1,9 +1,9 @@
 [tool.poetry]
 name = "tiempo_deshima"
-version = "0.1.16"
+version = "0.1.17"
 description = "time-dependent end-to-end model for post-process optimization"
 authors = ["Esmee Huijten"]
-maintainers = ["Stefanie Brackenhoff <[email protected]>", "Anne-Kee Doing <[email protected]>"]
+maintainers = ["Stefanie Brackenhoff <[email protected]>", "Anne-Kee Doing <[email protected]>"]
 license = "MIT"
 readme = "README.md"
 homepage = "https://github.com/deshima-dev/tiempo_deshima"

diff --git a/tiempo_deshima/DESHIMA/desim/minidesim.py b/tiempo_deshima/DESHIMA/desim/minidesim.py
@@ -359,7 +359,7 @@ def eta_atm_func(F, pwv, EL=60., eta_atm_df = pd.Series([]), F_highres = pd.Seri
     pwv = np.squeeze(pwv)
     eta_atm = np.zeros((len(F), len(pwv)))
     for i in range(len(pwv)):
-        eta_atm = np.abs(eta_atm_func_zenith(pwv[i], F)) ** (1./np.sin(EL*np.pi/180.))
+        eta_atm[:,i] = np.abs(np.squeeze(eta_atm_func_zenith(pwv[i], F))) ** (1./np.sin(EL*np.pi/180.)) #squeeze removes additional dimension of size 1, such that only 2 indices are needed
     return np.squeeze(eta_atm)