From 326bcdb2df1a38c7f244d3d2ed7dc18ec7fcac2b Mon Sep 17 00:00:00 2001
From: chrisbrahms <c.brahms@hw.ac.uk>
Date: Wed, 9 Aug 2023 14:43:49 +0100
Subject: [PATCH] add prop_gnlse case

---
 src/Luna.jl | 15 +++++++++++++++
 1 file changed, 15 insertions(+)

diff --git a/src/Luna.jl b/src/Luna.jl
index 792da1a8..c452cdce 100644
--- a/src/Luna.jl
+++ b/src/Luna.jl
@@ -381,6 +381,21 @@ Logging.with_logger(Logging.NullLogger()) do
                     τfwhm=10e-15, λlims=(150e-9, 4e-6), trange=1e-12, saveN=11,
                     modes=4)
     p = Tools.capillary_params(120e-6, 10e-15, 800e-9, 125e-6, :He, P=1.0)
+
+    # gnlse_sol.jl example but with N=1 and 100th of the fibre length
+    γ = 0.1
+    β2 = -1e-26
+    τ0 = 280e-15
+    τfwhm = (2*log(1 + sqrt(2)))*τ0
+    fr = 0.18
+    P0 = abs(β2)/((1-fr)*γ*τ0^2)
+    flength = π*τ0^2/abs(β2)/100
+    βs =  [0.0, 0.0, β2]
+    λ0 = 835e-9
+    λlims = [450e-9, 8000e-9]
+    trange = 4e-12
+    output = prop_gnlse(γ, flength, βs; λ0, τfwhm, power=P0, pulseshape=:sech, λlims, trange,
+                        raman=false, shock=false, fr, shotnoise=false, saveN=11)
 end
 
 end # module