switch to vector advection form and make stepper 4x1 for compressible…

… Eady
firedrakeproject · May 1, 2024 · df1659e · df1659e
1 parent 27d4074
commit df1659e
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/examples/compressible/compressible_eady.py b/examples/compressible/compressible_eady.py
@@ -45,7 +45,8 @@
 # Equation
 Omega = as_vector([0., 0., f*0.5])
 parameters = CompressibleEadyParameters(N=sqrt(2.5e-5), H=H, f=f, Pi0=Constant(1.0))
-eqns = CompressibleEadyEquations(domain, parameters, Omega=Omega)
+eqns = CompressibleEadyEquations(domain, parameters, Omega=Omega,
+                                 u_transport_option='vector_advection_form')
 
 # I/O
 output = OutputParameters(dirname=dirname,
@@ -84,7 +85,8 @@
 # Time stepper
 stepper = SemiImplicitQuasiNewton(eqns, io, transport_schemes,
                                   transport_methods,
-                                  linear_solver=linear_solver)
+                                  linear_solver=linear_solver,
+                                  num_outer=4, num_inner=1)
 
 # ---------------------------------------------------------------------------- #
 # Initial conditions