Skip to content

Commit

Permalink
kgos for incompressile and compressible boussinesq test
Browse files Browse the repository at this point in the history
  • Loading branch information
jshipton committed Mar 22, 2024
1 parent 9fc53d7 commit 3e4308d
Show file tree
Hide file tree
Showing 3 changed files with 5 additions and 6 deletions.
Binary file not shown.
Binary file not shown.
11 changes: 5 additions & 6 deletions integration-tests/equations/test_boussinesq.py
Original file line number Diff line number Diff line change
Expand Up @@ -25,7 +25,7 @@ def run_boussinesq(tmpdir, compressible):
ncols = 10 # number of columns
Lx = 1000.0
Lz = 1000.0
mesh_name = 'incompressible_mesh'
mesh_name = 'boussinesq_mesh'
m = PeriodicIntervalMesh(ncols, Lx)
mesh = ExtrudedMesh(m, layers=nlayers, layer_height=Lz/nlayers, name=mesh_name)
domain = Domain(mesh, dt, "CG", 1)
Expand Down Expand Up @@ -64,8 +64,7 @@ def run_boussinesq(tmpdir, compressible):
# Time stepper
stepper = SemiImplicitQuasiNewton(eqn, io, transported_fields,
transport_methods,
linear_solver=linear_solver,
num_outer=4, num_inner=1)
linear_solver=linear_solver)

# ------------------------------------------------------------------------ #
# Initial conditions
Expand Down Expand Up @@ -100,9 +99,9 @@ def run_boussinesq(tmpdir, compressible):

# State for checking checkpoints
if compressible:
checkpoint_name = 'compressible_chkpt.h5'
checkpoint_name = 'compressible_boussinesq_chkpt.h5'
else:
checkpoint_name = 'incompressible_chkpt.h5'
checkpoint_name = 'incompressible_boussinesq_chkpt.h5'
new_path = join(abspath(dirname(__file__)), '..', f'data/{checkpoint_name}')
check_output = OutputParameters(dirname=output_dirname,
checkpoint_pickup_filename=new_path,
Expand All @@ -117,7 +116,7 @@ def run_boussinesq(tmpdir, compressible):
return stepper, check_stepper


@pytest.mark.parametrize("compressible", [False])
@pytest.mark.parametrize("compressible", [True, False])
def test_boussinesq(tmpdir, compressible):

dirname = str(tmpdir)
Expand Down

0 comments on commit 3e4308d

Please sign in to comment.