Skip to content

GitLab

Switch branch/tag
Find file Normal viewHistoryPermalink
diffusion_test.py 1.73 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
50
51
52
53
54
55
56
57
58
59
60
61 
from veros.core import diffusion
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_cyclic_x=True,
    enable_conserve_energy=True,
    enable_hor_friction_cos_scaling=True,
    enable_tempsalt_sources=True,
    K_hbi=1,
    K_h=1,
    hor_friction_cosPower=2,
)


def prepare_inputs(vs_state, pyom_obj):
    # implementations are only identical if non-water values are 0
    vs = vs_state.variables
    for var in (
        "P_diss_sources",
        "P_diss_hmix",
    ):
        getattr(pyom_obj.main_module, var.lower())[...] *= vs.maskT
        with vs.unlock():
            setattr(vs, var, vs.get(var) * vs.maskT)

    return vs_state, pyom_obj


def test_tempsalt_biharmonic(pyom2_lib):
    vs_state, pyom_obj = get_random_state(pyom2_lib, extra_settings=TEST_SETTINGS)
    vs_state, pyom_obj = prepare_inputs(vs_state, pyom_obj)

    vs_state.variables.update(diffusion.tempsalt_biharmonic(vs_state))
    pyom_obj.tempsalt_biharmonic()
    compare_state(vs_state, pyom_obj)


def test_tempsalt_diffusion(pyom2_lib):
    vs_state, pyom_obj = get_random_state(pyom2_lib, extra_settings=TEST_SETTINGS)
    vs_state, pyom_obj = prepare_inputs(vs_state, pyom_obj)

    vs_state.variables.update(diffusion.tempsalt_diffusion(vs_state))
    pyom_obj.tempsalt_diffusion()
    compare_state(vs_state, pyom_obj)


def test_tempsalt_sources(pyom2_lib):
    vs_state, pyom_obj = get_random_state(pyom2_lib, extra_settings=TEST_SETTINGS)
    vs_state, pyom_obj = prepare_inputs(vs_state, pyom_obj)

    vs_state.variables.update(diffusion.tempsalt_sources(vs_state))
    pyom_obj.tempsalt_sources()
    compare_state(vs_state, pyom_obj)