From bcd4b9ae4a8f704927be069a52cd4793be63c29c Mon Sep 17 00:00:00 2001
From: Mathieu Roule <72818492+MathieuRoule@users.noreply.github.com>
Date: Fri, 14 Jun 2024 21:55:14 +0200
Subject: [PATCH] Update WMat.jl

Fix wrong call to OrbitalElements functions
---
 src/WMat.jl | 12 ++++++------
 1 file changed, 6 insertions(+), 6 deletions(-)

diff --git a/src/WMat.jl b/src/WMat.jl
index 2dfbc34..45803ae 100644
--- a/src/WMat.jl
+++ b/src/WMat.jl
@@ -61,19 +61,19 @@ function angle_fouriertransform!(
     @assert length(result) == output_length "Result array length does not 
         match the function's output length."
 
+    # Integration step
+    Kwp = (params.ADAPTIVEKW) ? ceil(Int64,params.Kw/(0.1+(1-e))) : params.Kw
+    Orbitalparams = params.Orbitalparams
+    
     if L == 0.0
         # need angular momentum
-        _, L = EL_from_ae(a, e, model, params)
+        _, L = EL_from_ae(a, e, model, Orbitalparams)
     end
     if Ω1 == 0.0 || Ω2 == 0.0
         # need frequencies
-        Ω1, Ω2 = frequencies_from_ae(a, e, model, params)
+        Ω1, Ω2 = frequencies_from_ae(a, e, model, Orbitalparams)
     end
 
-    # Integration step
-    Kwp = (params.ADAPTIVEKW) ? ceil(Int64,params.Kw/(0.1+(1-e))) : params.Kw
-    Orbitalparams = params.Orbitalparams
-
     # Caution : Reverse integration (lower error at apocenter than pericenter)
     # -> Result to multiply by -1
     dw = -2/Kwp