From 69cc2ffb24264f45bfc7bc5d9b3ae6b1f055e44e Mon Sep 17 00:00:00 2001
From: Florian Atteneder <florian.atteneder@uni-jena.de>
Date: Mon, 31 Jan 2022 13:09:21 +0000
Subject: [PATCH] fix wrong GW frequency cut off for FFI method
---
watpy/wave/wave.py | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/watpy/wave/wave.py b/watpy/wave/wave.py
index 9ad911d..9511ebc 100644
--- a/watpy/wave/wave.py
+++ b/watpy/wave/wave.py
@@ -618,7 +618,7 @@ class wave(object):
"""
if self.prop['var']=='Psi4':
if fcut < 0. :
- fcut = 2 * self.prop['init.frequency'] / max(1,abs(self.prop['mmode']))
+ fcut = max(1,abs(self.prop['mmode'])) * self.prop['init.frequency'] / 2
dt = self.time[1] - self.time[0]
return win * fixed_freq_int_2( win * self.p4, fcut, dt = dt)
else:
--
GitLab