Skip to content

GitLab

Switch branch/tag
Find file Normal viewHistoryPermalink
momentum_test.py 1.9 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
62
63 
from veros.core import momentum
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,
    coord_degree=True,
    enable_cyclic_x=True,
    enable_conserve_energy=True,
    enable_bottom_friction_var=True,
    enable_hor_friction_cos_scaling=True,
    enable_implicit_vert_friction=True,
    enable_explicit_vert_friction=True,
    enable_TEM_friction=True,
    enable_hor_friction=True,
    enable_biharmonic_friction=True,
    enable_ray_friction=True,
    enable_bottom_friction=True,
    enable_quadratic_bottom_friction=True,
    enable_momentum_sources=True,
)


def test_momentum_advection(pyom2_lib):
    vs_state, pyom_obj = get_random_state(pyom2_lib, extra_settings=TEST_SETTINGS)
    vs_state.variables.update(momentum.momentum_advection(vs_state))
    pyom_obj.momentum_advection()

    # not a part of momentum_advection in PyOM
    m = pyom_obj.main_module
    m.du[..., m.tau - 1] += m.du_adv
    m.dv[..., m.tau - 1] += m.dv_adv

    compare_state(vs_state, pyom_obj)


def test_vertical_velocity(pyom2_lib):
    vs_state, pyom_obj = get_random_state(pyom2_lib, extra_settings=TEST_SETTINGS)
    vs_state.variables.update(momentum.vertical_velocity(vs_state))
    pyom_obj.vertical_velocity()
    compare_state(vs_state, pyom_obj)


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

    # results are only identical if initial guess is already cyclic
    from veros.core import utilities

    vs = vs_state.variables
    m = pyom_obj.main_module
    m.psi[...] = utilities.enforce_boundaries(m.psi, vs_state.settings.enable_cyclic_x)
    vs.psi = utilities.enforce_boundaries(vs.psi, vs_state.settings.enable_cyclic_x)

    vs_state.variables.update(momentum.momentum(vs_state))
    pyom_obj.momentum()
    compare_state(vs_state, pyom_obj)