Merge pull request #478 from firedrakeproject/fix_moist_thermal_example

Fix moist thermal example
firedrakeproject · Feb 29, 2024 · 6ac3f4b · 6ac3f4b
2 parents 605064a + d18f0ec
commit 6ac3f4b
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/examples/shallow_water/moist_thermal_williamson5.py b/examples/shallow_water/moist_thermal_williamson5.py
@@ -34,8 +34,8 @@
 mu2 = 0.98
 q0 = 135  # chosen to give an initial max vapour of approx 0.02
 beta2 = 10
-qprecip = 10e-4
-gamma_r = 10e-3
+qprecip = 1e-4
+gamma_r = 1e-3
 # topography parameters
 R0 = pi/9.
 R0sq = R0**2
@@ -104,8 +104,10 @@ def gamma_v(x_in):
 InstantRain(eqns, qprecip, vapour_name="cloud_water", rain_name="rain",
             gamma_r=gamma_r)
 
+transport_methods = [DGUpwind(eqns, field_name) for field_name in eqns.field_names]
+
 # Timestepper
-stepper = Timestepper(eqns, RK4(domain), io)
+stepper = Timestepper(eqns, RK4(domain), io, spatial_methods=transport_methods)
 
 # ----------------------------------------------------------------- #
 # Initial conditions