SSPRK2/3 stability tests

Author: milankl

test/runtests.jl

@@ -38,3 +38,39 @@ end
     @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
     @test all(Prog.u .!= 0.0f0)
 end
+
+@testset "SSPRK2 tests" begin
+    Prog = RunModel(time_scheme="SSPRK2",RKs=2,cfl=0.7,nstep_advcor=1)
+    @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
+    @test all(Prog.u .!= 0.0f0)
+
+    Prog = RunModel(time_scheme="SSPRK2",RKs=3,cfl=1.2,nstep_advcor=1)
+    @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
+    @test all(Prog.u .!= 0.0f0)
+
+    Prog = RunModel(time_scheme="SSPRK2",RKs=4,cfl=1.4,nstep_advcor=1)
+    @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
+    @test all(Prog.u .!= 0.0f0)
+
+    Prog = RunModel(time_scheme="SSPRK2",RKs=2,cfl=0.7,nstep_advcor=0)
+    @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
+    @test all(Prog.u .!= 0.0f0)
+
+    Prog = RunModel(time_scheme="SSPRK2",RKs=3,cfl=1.2,nstep_advcor=0)
+    @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
+    @test all(Prog.u .!= 0.0f0)
+
+    Prog = RunModel(time_scheme="SSPRK2",RKs=4,cfl=1.2,nstep_advcor=0)
+    @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
+    @test all(Prog.u .!= 0.0f0)
+end
+
+@testset "SSPRK3 tests" begin
+    Prog = RunModel(time_scheme="SSPRK3",cfl=1.2,nstep_advcor=1)
+    @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
+    @test all(Prog.u .!= 0.0f0)
+
+    Prog = RunModel(time_scheme="SSPRK3",cfl=1.2,nstep_advcor=0)
+    @test all(abs.(Prog.η) .< 10)    # sea surface height shouldn't exceed +-10m
+    @test all(Prog.u .!= 0.0f0)
+end