From df1659e13b7094a217df283799e7d3d70d3f06ee Mon Sep 17 00:00:00 2001 From: Tom Bendall Date: Wed, 1 May 2024 14:45:05 +0100 Subject: [PATCH] switch to vector advection form and make stepper 4x1 for compressible Eady --- examples/compressible/compressible_eady.py | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/examples/compressible/compressible_eady.py b/examples/compressible/compressible_eady.py index 1ccccd601..3a6752f2f 100644 --- 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