From f958d558866a884b3474c66ee528db19a054a43c Mon Sep 17 00:00:00 2001
From: Andrew Simms <simms.andrew@gmail.com>
Date: Mon, 1 Jul 2024 08:24:56 -0600
Subject: [PATCH] Test: Gamma parameter of jonswap_spectrum

---
 mhkit/tests/Wave_TestResourceSpectrum.m | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/mhkit/tests/Wave_TestResourceSpectrum.m b/mhkit/tests/Wave_TestResourceSpectrum.m
index 423180d72..83d9bee2c 100644
--- a/mhkit/tests/Wave_TestResourceSpectrum.m
+++ b/mhkit/tests/Wave_TestResourceSpectrum.m
@@ -113,6 +113,21 @@ function test_jonswap_spectrum(testCase)
             assertLessThan(testCase,errorHm0, 0.01);
             assertLessThan(testCase,errorTp0, 0.01);
         end
+
+        function test_jonswap_spectrum_gamma(testCase)
+            Obj.f = 0.1/(2*pi):0.01/(2*pi):3.5/(2*pi);
+            Obj.Tp = 8;
+            Obj.Hs = 2.5;
+            Obj.gamma = 2.0
+
+            S = jonswap_spectrum(Obj.f, Obj.Tp, Obj.Hs, Obj.gamma);
+            Hm0 = significant_wave_height(S);
+            Tp0 = peak_period(S);
+            errorHm0 = abs(Obj.Tp - Tp0)/Obj.Tp;
+            errorTp0 = abs(Obj.Hs - Hm0)/Obj.Hs;
+            assertLessThan(testCase,errorHm0, 0.01);
+            assertLessThan(testCase,errorTp0, 0.01);
+        end
         
         function test_plot_spectrum(testCase)
             Obj.f = 0.1/(2*pi):0.01/(2*pi):3.5/(2*pi);