← 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-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/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
  </variables>
 
76
    <pass_tests>  
 
77
    <test name="check end time" language="python">
 
78
assert endtime &gt; 0.35
 
79
    </test>
 
80
    <test name="check end timestep" language="python">
 
81
assert endtimestep &gt; 0.001
 
82
    </test>
 
83
    <test name="check water volume fraction was conserved" language="python">
 
84
print 'conservation error = ', abs(vfracintwater0-vfracintwater)
 
85
assert abs(vfracintwater0-vfracintwater) &lt; max(maxdivergence, mindivergence)
 
86
    </test>
 
87
    <test name="check droplet volume fraction was conserved" language="python">
 
88
print 'conservation error = ', abs(vfracintdrop0-vfracintdrop)
 
89
assert abs(vfracintdrop0-vfracintdrop) &lt; max(maxdivergence, mindivergence)
 
90
    </test>
 
91
    <test name="check that reasonable tolerance (1.0e-08) is being used to test conservation" language="python">
 
92
print 'divergence tolerance = ', max(maxdivergence, mindivergence)
 
93
assert max(maxdivergence, mindivergence) &lt; 1.0e-08
 
94
    </test>
 
95
    <test name="check minimum water volume fraction is 0.0" language="python">
 
96
assert abs(mnvfracwater) &lt; 1e-10
 
97
    </test>
 
98
    <test name="check maximum water volume fraction is 1.0" language="python">
 
99
assert abs(mxvfracwater-1.0) &lt; 1e-10
 
100
    </test>
 
101
    <test name="check minimum droplet volume fraction is 0.0" language="python">
 
102
assert abs(mnvfracdrop) &lt; 1e-10
 
103
    </test>
 
104
    <test name="check maximum droplet volume fraction is 1.0" language="python">
 
105
assert abs(mxvfracdrop-1.0) &lt; 1e-10
 
106
    </test>
 
107
    <test name="check minimum air volume fraction is 0.0" language="python">
 
108
assert abs(mnvfracair) &lt; 1e-10
 
109
    </test>
 
110
    <test name="check maximum air volume fraction is 1.0" language="python">
 
111
assert abs(mxvfracair-1.0) &lt; 1e-10
 
112
    </test>
 
113
    <test name="check minimum sum volume fraction is 0.0" language="python">
 
114
assert abs(mnvfracsum) &lt; 1e-10
 
115
    </test>
 
116
    <test name="check maximum sum volume fraction is 1.0" language="python">
 
117
assert abs(mxvfracsum-1.0) &lt; 1e-10
 
118
    </test>
 
119
  </pass_tests>
 
120
  <warn_tests>
 
121
  </warn_tests>
 
122
</testproblem>