Skip to content

GitLab

Switch branch/tag
Find file Normal viewHistoryPermalink
idemix_test.py 1.41 KB
Newer Older
mashun1's avatar
veros  
mashun1 committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49 
import pytest

from veros.core import idemix
from veros.pyom_compat import get_random_state

from test_base import compare_state


TEST_SETTINGS = dict(
    nx=70,
    ny=60,
    nz=50,
    dt_tracer=3600,
    dt_mom=3600,
    enable_idemix=True,
    enable_idemix_hor_diffusion=True,
    enable_idemix_superbee_advection=True,
    enable_idemix_upwind_advection=True,
    enable_eke=True,
    enable_store_cabbeling_heat=True,
    enable_eke_diss_bottom=True,
    enable_eke_diss_surfbot=True,
    enable_store_bottom_friction_tke=True,
    enable_TEM_friction=True,
)

PROBLEM_SETS = {
    "eke": dict(enable_eke=True),
    "no-eke": dict(enable_eke=False),
    "no-eke_diss_bottom": dict(enable_eke_diss_bottom=False),
    "no-eke_diss_surfbot": dict(enable_eke_diss_surfbot=False),
}


def test_set_idemix_parameter(pyom2_lib):
    vs_state, pyom_obj = get_random_state(pyom2_lib, extra_settings=TEST_SETTINGS)
    vs_state.variables.update(idemix.set_idemix_parameter(vs_state))
    pyom_obj.set_idemix_parameter()
    compare_state(vs_state, pyom_obj)


@pytest.mark.parametrize("problem_set", PROBLEM_SETS)
def test_integrate_idemix(pyom2_lib, problem_set):
    settings = {**TEST_SETTINGS, **PROBLEM_SETS[problem_set]}
    vs_state, pyom_obj = get_random_state(pyom2_lib, extra_settings=settings)

    vs_state.variables.update(idemix.integrate_idemix(vs_state))
    pyom_obj.integrate_idemix()
    compare_state(vs_state, pyom_obj)