← 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-3d/3material-droplet-p1dgp2-test-cty-cv-rhop1dg-3d.xml

  • Committer: Simon Mouradian
  • Date: 2012-05-30 09:08:23 UTC
  • mfrom: (3520.32.274 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20120530090823-pq0fgemjpqn72j8d
merge from lp:fluidity

Show diffs side-by-side

added added

removed removed

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