6
6
<tags>flml gls</tags>
7
7
<problem_definition length="medium" nprocs="1">
9
fluidity gls-Kato_Phillips-mixed_layer_depth-k_e-CA.flml
10
fluidity gls-Kato_Phillips-mixed_layer_depth-k_kl-KC94.flml
9
fluidity gls-Kato_Phillips-mixed_layer_depth-gen-GL.flml
10
fluidity gls-Kato_Phillips-mixed_layer_depth-k_w-CB.flml
12
12
</problem_definition>
15
<variable name="FinalTimeKECA" language="python">
17
FinalTimeKECA = fluidity_tools.stat_parser("Kato_Phillips-mld-k_e-CA.stat")['ElapsedTime']['value'][-1]
19
<variable name="FinalTimeKKLKC" language="python">
21
FinalTimeKKLKC = fluidity_tools.stat_parser("Kato_Phillips-mld-k_kl-KC.stat")['ElapsedTime']['value'][-1]
15
<variable name="FinalTimeGenGL" language="python">
17
FinalTimeGenGL = fluidity_tools.stat_parser("Kato_Phillips-mld-gen-GL.stat")['ElapsedTime']['value'][-1]
19
<variable name="FinalTimeKWCB" language="python">
21
FinalTimeKWCB = fluidity_tools.stat_parser("Kato_Phillips-mld-k_w-CB.stat")['ElapsedTime']['value'][-1]
24
<variable name="KECA_MLD" language="python">
26
sys.path.append("gls-Kato_Phillips-mixed_layer_depth")
29
filelist = glob.glob("Kato_Phillips-mld-k_e-CA*.vtu")
30
mld_calc.sort_nicely(filelist)
31
KECA_MLD = mld_calc.MLD(filelist)
33
<variable name="KKLKC_MLD" language="python">
35
sys.path.append("gls-Kato_Phillips-mixed_layer_depth")
38
filelist = glob.glob("Kato_Phillips-mld-k_kl-KC_*.vtu")
39
mld_calc.sort_nicely(filelist)
40
KKLKC_MLD = mld_calc.MLD(filelist)
24
<variable name="GenGL_MLD" language="python">
26
sys.path.append("gls-Kato_Phillips-mixed_layer_depth")
29
filelist = glob.glob("Kato_Phillips-mld-gen-GL_*.vtu")
30
mld_calc.sort_nicely(filelist)
31
GenGL_MLD = mld_calc.MLD(filelist)
33
<variable name="KWCB_MLD" language="python">
35
sys.path.append("gls-Kato_Phillips-mixed_layer_depth")
38
filelist = glob.glob("Kato_Phillips-mld-k_w-CB_*.vtu")
39
mld_calc.sort_nicely(filelist)
40
KWCB_MLD = mld_calc.MLD(filelist)
46
<test name="Final time equals 10 hours for k-e CA" language="python">
47
assert FinalTimeKECA > 17900
46
<test name="Final time equals 10 hours for Gen GL" language="python">
47
assert FinalTimeGenGL > 35900
49
<test name="Final time equals 10 hours for k-kl KC" language="python">
50
assert FinalTimeKKLKC > 17900
49
<test name="Final time equals 10 hours for k-w CB" language="python">
50
assert FinalTimeKWCB > 35900
53
<test name ="MLD for K-E KC vs analytical" language="python">
54
from numpy import array, all
55
assert (all( abs(array(KECA_MLD[1]) - array(KECA_MLD[2]) ) < 1.5))
57
<test name ="MLD for K-KL vs analytical" language="python">
58
from numpy import array, all
59
assert (all( abs(array(KKLKC_MLD[1]) - array(KKLKC_MLD[2]) ) < 1.5))
53
<test name ="MLD for Gen GL vs analytical" language="python">
54
from numpy import array, all
55
assert (all( abs(array(GenGL_MLD[1]) - array(GenGL_MLD[2]) ) < 2.0))
57
<test name ="MLD for K-W CB vs analytical" language="python">
58
from numpy import array, all
59
print abs(array(KWCB_MLD[1]) - array(KWCB_MLD[2]) )
60
assert (all( abs(array(KWCB_MLD[1]) - array(KWCB_MLD[2]) ) < 2.0))