From cda2ab9eaff80b98204d6b8756c3ca7e97727369 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Mikkel=20Friis-M=C3=B8ller?= <mikf@dtu.dk>
Date: Fri, 2 Jul 2021 10:40:46 +0000
Subject: [PATCH] Tilt defl
---
py_wake/deflection_models/jimenez.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/py_wake/deflection_models/jimenez.py b/py_wake/deflection_models/jimenez.py
index a8f4b525d..1bb475f41 100644
--- a/py_wake/deflection_models/jimenez.py
+++ b/py_wake/deflection_models/jimenez.py
@@ -25,8 +25,8 @@ class JimenezWakeDeflection(DeflectionModel):
nominator_ijxl = (1 + (self.beta / D_src_il)[:, na, na, :] * dw_ijxl)**2
alpha = denominator_ilk[:, na, na] / nominator_ijxl[..., na]
deflection_ijlk = np.trapz(np.sin(alpha), dw_ijxl[..., na], axis=2)
- self.hcw_ijlk = hcw_ijl[..., na] + deflection_ijlk * np.cos(theta_deflection_ilk[..., na])
- self.dh_ijlk = dh_ijl[..., na] + deflection_ijlk * np.sin(theta_deflection_ilk[..., na])
+ self.hcw_ijlk = hcw_ijl[..., na] + deflection_ijlk * np.cos(theta_deflection_ilk[:, na])
+ self.dh_ijlk = dh_ijl[..., na] + deflection_ijlk * np.sin(theta_deflection_ilk[:, na])
return dw_ijl[..., na], self.hcw_ijlk, self.dh_ijlk
--
GitLab