Add unit and type tests (#117)

* add unit tests for ec fluxes

* add test unit for tec, etec and lmars fluxe

* include test_unit.jl

* add unit tests for shallow water

* fix formatting

* increase coverage

* format

* change pressure computation

* add first type stable tests

* add some more type stable tests

* add all tests for potential temperature

* fix formatting

* add type stable test for moist euler

* fix exa
mples

* fix typo

* add ShalloWater 3D tests

* remove @show

* update test_type.jl

* update test_type.jl

* update runtests.jl

* increase coverage

* apply code review suggestions

* fix elixirs

* apply review suggestions

---------

Co-authored-by: Hendrik Ranocha <ranocha@users.noreply.github.com>

Author: MarcoArtiano

examples/elixir_euler_potential_temperature_baroclinic_instability.jl

@@ -193,7 +193,9 @@ end
 
     return SVector(du1, du2, du3, du4, du5, zero(eltype(u)))
 end
-equations = CompressibleEulerPotentialTemperatureEquationsWithGravity3D()
+equations = CompressibleEulerPotentialTemperatureEquationsWithGravity3D(c_p = 1004,
+                                                                        c_v = 717,
+                                                                        gravity = 9.81)
 
 initial_condition = initial_condition_baroclinic_instability
 

examples/elixir_euler_potential_temperature_ec.jl

@@ -9,7 +9,7 @@ function initial_condition_density_wave(x, t,
     return prim2cons(SVector(rho, v1, p), equations)
 end
 
-equations = CompressibleEulerPotentialTemperatureEquations1D()
+equations = CompressibleEulerPotentialTemperatureEquations1D(c_p = 1004, c_v = 717)
 
 polydeg = 3
 basis = LobattoLegendreBasis(polydeg)

examples/elixir_euler_potential_temperature_inertia_gravity_waves.jl

@@ -36,7 +36,9 @@ function initial_condition_gravity_waves(x, t,
     return prim2cons(SVector(rho, v1, v2, p, g * x[2]), equations)
 end
 
-equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D()
+equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D(c_p = 1004,
+                                                                        c_v = 717,
+                                                                        gravity = 9.81)
 
 # We have an isothermal background state with T0 = 250 K. 
 # The reference speed of sound can be computed as:

examples/elixir_euler_potential_temperature_linear_hydrostatic.jl

@@ -95,7 +95,9 @@ end
     return prim2cons(SVector(rho, v1, v2, p, g * x[2]), equations)
 end
 
-equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D()
+equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D(c_p = 1004,
+                                                                        c_v = 717,
+                                                                        gravity = 9.81)
 alpha = 0.035
 xr_B = 60000.0
 linear_hydrostatic_setup = HydrostaticSetup(alpha, xr_B, equations)

examples/elixir_euler_potential_temperature_linear_nonhydrostatic.jl

@@ -91,7 +91,9 @@ end
     return prim2cons(SVector(rho, v1, v2, p, g * x[2]), equations)
 end
 
-equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D()
+equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D(c_p = 1004,
+                                                                        c_v = 717,
+                                                                        gravity = 9.81)
 alpha = 0.03
 xr_B = 40000.0
 

examples/elixir_euler_potential_temperature_robert_bubble.jl

@@ -51,7 +51,7 @@ end
     return SVector(zero(eltype(u)), zero(eltype(u)), -g * rho, zero(eltype(u)))
 end
 
-equations = CompressibleEulerPotentialTemperatureEquations2D()
+equations = CompressibleEulerPotentialTemperatureEquations2D(c_p = 1004, c_v = 717)
 
 surface_flux = FluxLMARS(340)
 volume_flux = flux_tec

examples/elixir_euler_potential_temperature_schaer_mountain.jl

@@ -91,7 +91,9 @@ end
 
 ###############################################################################
 # semidiscretization of the compressible Euler equations
-equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D()
+equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D(c_p = 1004,
+                                                                        c_v = 717,
+                                                                        gravity = 9.81)
 alpha = 0.03
 xr_B = 20000
 schär_setup = SchärSetup(alpha, xr_B)

examples/elixir_euler_potential_temperature_taylor_green_vortex.jl

@@ -14,7 +14,7 @@ function initial_condition_taylor_green_vortex(x, t,
     return prim2cons(SVector(rho, v1, v2, v3, p), equations)
 end
 
-equations = CompressibleEulerPotentialTemperatureEquations3D()
+equations = CompressibleEulerPotentialTemperatureEquations3D(c_p = 1004, c_v = 717)
 polydeg = 3
 basis = LobattoLegendreBasis(polydeg)
 surface_flux = flux_etec

examples/elixir_euler_potential_temperature_well_balanced.jl

@@ -32,7 +32,9 @@ function initial_condition_adiabatic(x, t,
     return prim2cons(SVector(rho, v1, p, g * x[1]), equations)
 end
 
-equations = CompressibleEulerPotentialTemperatureEquationsWithGravity1D()
+equations = CompressibleEulerPotentialTemperatureEquationsWithGravity1D(c_p = 1004,
+                                                                        c_v = 717,
+                                                                        gravity = 9.81)
 polydeg = 2
 basis = LobattoLegendreBasis(polydeg)
 cs = 340.0

examples/elixir_euler_potential_temperature_well_balanced_curvilinear.jl

@@ -33,7 +33,9 @@ function initial_condition_adiabatic(x, t,
     return prim2cons(SVector(rho, v1, v2, p, g * x[2]), equations)
 end
 
-equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D()
+equations = CompressibleEulerPotentialTemperatureEquationsWithGravity2D(c_p = 1004,
+                                                                        c_v = 717,
+                                                                        gravity = 9.81)
 polydeg = 2
 basis = LobattoLegendreBasis(polydeg)
 cs = 340