← Back to branch summary

~fluidity-core/fluidity/sea-ice-branch

« back to all changes in this revision

Viewing changes to tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg/3material-droplet-p1dgp2-test-cty-cv-rhop1dg.xml

  • Committer: Simon Mouradian
  • Date: 2012-03-30 11:11:47 UTC
  • mfrom: (3520.32.239 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20120330111147-y4n4y4th52x2dqcf
merge from lp:fluidity

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/3material-droplet-p1dgp2-test-cty-cv-rhop1dg/3material-droplet-p1dgp2-test-cty-cv-rhop1dg.xml
1
 
<?xml version="1.0" encoding="UTF-8" ?>
2
 
 
3
 
<testproblem>
4
 
  <name>p1dgp2-test-cty-cv-rhop1dg multi-material flow test</name>
5
 
  <owner userid="btollit"/>
6
 
  <tags>flml</tags>
7
 
    <problem_definition length="medium" nprocs="1">
8
 
    <command_line>fluidity -v1 -l 3material-droplet-p1dgp2-test-cty-cv-rhop1dg.flml </command_line>
9
 
  </problem_definition>  
10
 
  <variables>  
11
 
    <variable name="endtime" language="python">
12
 
from fluidity_tools import stat_parser as stat
13
 
endtime=stat("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")["ElapsedTime"]["value"][-1]
14
 
    </variable>
15
 
    <variable name="endtimestep" language="python">
16
 
from fluidity_tools import stat_parser as stat
17
 
endtimestep=stat("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")["dt"]["value"][-1]
18
 
    </variable>
19
 
    <variable name="mindivergence" language="python">
20
 
from fluidity_tools import stat_parser as stat
21
 
mindivergence=max(abs(stat("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")["Water"]["ControlVolumeDivergence"]["min"]))
22
 
    </variable>
23
 
    <variable name="maxdivergence" language="python">
24
 
from fluidity_tools import stat_parser as stat
25
 
maxdivergence=max(stat("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")["Water"]["ControlVolumeDivergence"]["max"])
26
 
    </variable>
27
 
     <variable name="vfracintwater0" language="python"> 
28
 
import fluidity_tools
29
 
vfracintwater0 = fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Water']['MaterialVolumeFraction']['cv_integral'][0]
30
 
    </variable>    
31
 
     <variable name="vfracintwater" language="python"> 
32
 
import fluidity_tools
33
 
vfracintwater = fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Water']['MaterialVolumeFraction']['cv_integral'][-1]
34
 
    </variable>    
35
 
     <variable name="vfracintdrop0" language="python"> 
36
 
import fluidity_tools
37
 
vfracintdrop0 = fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Droplet']['MaterialVolumeFraction']['cv_integral'][0]
38
 
    </variable>    
39
 
     <variable name="vfracintdrop" language="python"> 
40
 
import fluidity_tools
41
 
vfracintdrop = fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Droplet']['MaterialVolumeFraction']['cv_integral'][-1]
42
 
    </variable>    
43
 
     <variable name="mxvfracwater" language="python"> 
44
 
import fluidity_tools
45
 
mxvfracwater = min(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Water']['MaterialVolumeFraction']['max'])
46
 
    </variable>    
47
 
     <variable name="mnvfracwater" language="python"> 
48
 
import fluidity_tools
49
 
mnvfracwater = min(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Water']['MaterialVolumeFraction']['min'])
50
 
    </variable>    
51
 
    <variable name="mxvfracdrop" language="python">
52
 
import fluidity_tools
53
 
mxvfracdrop = max(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Droplet']['MaterialVolumeFraction']['max'])
54
 
    </variable>
55
 
     <variable name="mnvfracdrop" language="python">
56
 
import fluidity_tools
57
 
mnvfracdrop = min(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Droplet']['MaterialVolumeFraction']['min'])
58
 
    </variable>
59
 
     <variable name="mxvfracair" language="python"> 
60
 
import fluidity_tools
61
 
mxvfracair = max(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Air']['MaterialVolumeFraction']['max'])
62
 
    </variable>    
63
 
     <variable name="mnvfracair" language="python"> 
64
 
import fluidity_tools
65
 
mnvfracair = min(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Air']['MaterialVolumeFraction']['min'])
66
 
    </variable>    
67
 
    <variable name="mxvfracsum" language="python">
68
 
import fluidity_tools
69
 
mxvfracsum = max(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Water']['SumMaterialVolumeFractions']['max'])
70
 
    </variable>
71
 
     <variable name="mnvfracsum" language="python">
72
 
import fluidity_tools
73
 
mnvfracsum = min(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Water']['SumMaterialVolumeFractions']['min'])
74
 
    </variable>
75
 
     <variable name="max_reynolds" language="python">
76
 
import fluidity_tools
77
 
max_reynolds = max(fluidity_tools.stat_parser("3material-droplet-p1dgp2-test-cty-cv-rhop1dg.stat")['Water']['GridReynoldsNumber']['max'][1:])
78
 
    </variable>    
79
 
  </variables>
80
 
    <pass_tests>  
81
 
    <test name="check end time" language="python">
82
 
assert endtime &gt; 0.35
83
 
    </test>
84
 
    <test name="check end timestep" language="python">
85
 
assert endtimestep &gt; 0.001
86
 
    </test>
87
 
    <test name="check water volume fraction was conserved" language="python">
88
 
print 'conservation error = ', abs(vfracintwater0-vfracintwater)
89
 
assert abs(vfracintwater0-vfracintwater) &lt; max(maxdivergence, mindivergence)
90
 
    </test>
91
 
    <test name="check droplet volume fraction was conserved" language="python">
92
 
print 'conservation error = ', abs(vfracintdrop0-vfracintdrop)
93
 
assert abs(vfracintdrop0-vfracintdrop) &lt; max(maxdivergence, mindivergence)
94
 
    </test>
95
 
    <test name="check that reasonable tolerance (1.0e-08) is being used to test conservation" language="python">
96
 
print 'divergence tolerance = ', max(maxdivergence, mindivergence)
97
 
assert max(maxdivergence, mindivergence) &lt; 1.0e-08
98
 
    </test>
99
 
    <test name="check minimum water volume fraction is 0.0" language="python">
100
 
assert abs(mnvfracwater) &lt; 1e-10
101
 
    </test>
102
 
    <test name="check maximum water volume fraction is 1.0" language="python">
103
 
assert abs(mxvfracwater-1.0) &lt; 1e-10
104
 
    </test>
105
 
    <test name="check minimum droplet volume fraction is 0.0" language="python">
106
 
assert abs(mnvfracdrop) &lt; 1e-10
107
 
    </test>
108
 
    <test name="check maximum droplet volume fraction is 1.0" language="python">
109
 
assert abs(mxvfracdrop-1.0) &lt; 1e-10
110
 
    </test>
111
 
    <test name="check minimum air volume fraction is 0.0" language="python">
112
 
assert abs(mnvfracair) &lt; 1e-10
113
 
    </test>
114
 
    <test name="check maximum air volume fraction is 1.0" language="python">
115
 
assert abs(mxvfracair-1.0) &lt; 1e-10
116
 
    </test>
117
 
    <test name="check minimum sum volume fraction is 0.0" language="python">
118
 
assert abs(mnvfracsum) &lt; 1e-10
119
 
    </test>
120
 
    <test name="check maximum sum volume fraction is 1.0" language="python">
121
 
assert abs(mxvfracsum-1.0) &lt; 1e-10
122
 
    </test>
123
 
    <test name="check maximum grid reynolds number is 4599.44185868 to tolerance 2.5e-02" language="python">
124
 
assert abs(max_reynolds - 4599.44185868)/4599.44185868 &lt; 2.5e-02
125
 
    </test>    
126
 
  </pass_tests>
127
 
  <warn_tests>
128
 
  </warn_tests>
129
 
</testproblem>