← Back to branch summary

~fluidity-core/fluidity/refactor-netcdf

« back to all changes in this revision

Viewing changes to tests/sigma_layer_sphere_parallel/sphere_sigma.flml

  • Committer: Jon Hill
  • Date: 2013-02-16 09:01:40 UTC
  • mfrom: (3981.7.159 fluidity)
  • Revision ID: jon.hill@imperial.ac.uk-20130216090140-bplzxqzdk1eik4za
Megre from trunk, fixing several conflicts

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/sigma_layer_sphere_parallel/sphere_sigma.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">sigma_layers</string_value>
 
5
  </simulation_name>
 
6
  <problem_type>
 
7
    <string_value lines="1">oceans</string_value>
 
8
  </problem_type>
 
9
  <geometry>
 
10
    <dimension>
 
11
      <integer_value rank="0">3</integer_value>
 
12
    </dimension>
 
13
    <mesh name="CoordinateMesh">
 
14
      <from_mesh>
 
15
        <mesh name="BaseMesh"/>
 
16
        <stat>
 
17
          <include_in_stat/>
 
18
        </stat>
 
19
      </from_mesh>
 
20
    </mesh>
 
21
    <mesh name="PressureMesh">
 
22
      <from_mesh>
 
23
        <mesh name="BaseMesh"/>
 
24
        <mesh_shape>
 
25
          <polynomial_degree>
 
26
            <integer_value rank="0">2</integer_value>
 
27
          </polynomial_degree>
 
28
        </mesh_shape>
 
29
        <stat>
 
30
          <exclude_from_stat/>
 
31
        </stat>
 
32
      </from_mesh>
 
33
    </mesh>
 
34
    <mesh name="InputMesh">
 
35
      <from_file file_name="mesh/np_rectangle">
 
36
        <format name="triangle"/>
 
37
        <stat>
 
38
          <exclude_from_stat/>
 
39
        </stat>
 
40
      </from_file>
 
41
    </mesh>
 
42
    <mesh name="BaseMesh">
 
43
      <from_mesh>
 
44
        <mesh name="InputMesh"/>
 
45
        <extrude>
 
46
          <regions name="WholeMesh">
 
47
            <bottom_depth>
 
48
              <python>
 
49
                <string_value type="code" lines="20" language="python">def val(X,t):
 
50
  z = X[2]
 
51
  # make the depth a linear function of Z (which shows little
 
52
  # change over the NP, but X and Y do)
 
53
  # This should vary from 4000-ish to 1000-ish
 
54
  return (z-6.331010e+06) / 10
 
55
  </string_value>
 
56
              </python>
 
57
            </bottom_depth>
 
58
            <sizing_function>
 
59
              <sigma_layers>
 
60
                <standard>
 
61
                  <integer_value rank="0">2</integer_value>
 
62
                </standard>
 
63
              </sigma_layers>
 
64
            </sizing_function>
 
65
            <top_surface_id>
 
66
              <integer_value rank="0">1</integer_value>
 
67
            </top_surface_id>
 
68
            <bottom_surface_id>
 
69
              <integer_value rank="0">2</integer_value>
 
70
            </bottom_surface_id>
 
71
          </regions>
 
72
        </extrude>
 
73
        <stat>
 
74
          <include_in_stat/>
 
75
        </stat>
 
76
      </from_mesh>
 
77
    </mesh>
 
78
    <quadrature>
 
79
      <degree>
 
80
        <integer_value rank="0">4</integer_value>
 
81
      </degree>
 
82
      <surface_degree>
 
83
        <integer_value rank="0">3</integer_value>
 
84
      </surface_degree>
 
85
    </quadrature>
 
86
    <spherical_earth>
 
87
      <linear_mapping/>
 
88
    </spherical_earth>
 
89
    <ocean_boundaries>
 
90
      <top_surface_ids>
 
91
        <integer_value shape="1" rank="1">1</integer_value>
 
92
      </top_surface_ids>
 
93
      <bottom_surface_ids>
 
94
        <integer_value shape="1" rank="1">2</integer_value>
 
95
      </bottom_surface_ids>
 
96
      <scalar_field name="DistanceToTop" rank="0">
 
97
        <diagnostic>
 
98
          <algorithm name="Internal" material_phase_support="multiple"/>
 
99
          <mesh name="CoordinateMesh"/>
 
100
          <output/>
 
101
          <stat/>
 
102
          <convergence>
 
103
            <include_in_convergence/>
 
104
          </convergence>
 
105
          <detectors>
 
106
            <include_in_detectors/>
 
107
          </detectors>
 
108
          <steady_state>
 
109
            <include_in_steady_state/>
 
110
          </steady_state>
 
111
        </diagnostic>
 
112
      </scalar_field>
 
113
      <scalar_field name="DistanceToBottom" rank="0">
 
114
        <diagnostic>
 
115
          <algorithm name="Internal" material_phase_support="multiple"/>
 
116
          <mesh name="CoordinateMesh"/>
 
117
          <output/>
 
118
          <stat/>
 
119
          <convergence>
 
120
            <include_in_convergence/>
 
121
          </convergence>
 
122
          <detectors>
 
123
            <include_in_detectors/>
 
124
          </detectors>
 
125
          <steady_state>
 
126
            <include_in_steady_state/>
 
127
          </steady_state>
 
128
        </diagnostic>
 
129
      </scalar_field>
 
130
    </ocean_boundaries>
 
131
  </geometry>
 
132
  <io>
 
133
    <dump_format>
 
134
      <string_value>vtk</string_value>
 
135
    </dump_format>
 
136
    <dump_period>
 
137
      <constant>
 
138
        <real_value rank="0">864000</real_value>
 
139
        <comment>10*86400 - every 10 days</comment>
 
140
      </constant>
 
141
    </dump_period>
 
142
    <output_mesh name="CoordinateMesh"/>
 
143
    <stat/>
 
144
  </io>
 
145
  <timestepping>
 
146
    <current_time>
 
147
      <real_value rank="0">0.0</real_value>
 
148
      <time_units date="seconds since 1987-01-05 00:00.0"/>
 
149
    </current_time>
 
150
    <timestep>
 
151
      <real_value rank="0">360.0</real_value>
 
152
    </timestep>
 
153
    <finish_time>
 
154
      <real_value rank="0">720</real_value>
 
155
    </finish_time>
 
156
    <nonlinear_iterations>
 
157
      <integer_value rank="0">2</integer_value>
 
158
    </nonlinear_iterations>
 
159
  </timestepping>
 
160
  <physical_parameters>
 
161
    <gravity>
 
162
      <magnitude>
 
163
        <real_value rank="0">9.81</real_value>
 
164
      </magnitude>
 
165
      <vector_field name="GravityDirection" rank="1">
 
166
        <prescribed>
 
167
          <mesh name="CoordinateMesh"/>
 
168
          <value name="WholeMesh">
 
169
            <python>
 
170
              <string_value lines="20" type="code" language="python">def val(X, t):
 
171
 
 
172
  a = X[0]
 
173
  b = X[1]
 
174
  c = X[2]
 
175
  
 
176
  x_component = -(a/((a**2+b**2+c**2)**0.5)) 
 
177
  y_component = -(b/((a**2+b**2+c**2)**0.5))
 
178
  z_component = -(c/((a**2+b**2+c**2)**0.5))
 
179
  
 
180
  return [x_component, y_component, z_component]</string_value>
 
181
            </python>
 
182
          </value>
 
183
          <output>
 
184
            <exclude_from_vtu/>
 
185
          </output>
 
186
          <stat>
 
187
            <include_in_stat/>
 
188
          </stat>
 
189
          <detectors>
 
190
            <exclude_from_detectors/>
 
191
          </detectors>
 
192
        </prescribed>
 
193
      </vector_field>
 
194
    </gravity>
 
195
  </physical_parameters>
 
196
  <material_phase name="Fields">
 
197
    <vector_field name="Velocity" rank="1">
 
198
      <prescribed>
 
199
        <mesh name="CoordinateMesh"/>
 
200
        <value name="WholeMesh">
 
201
          <constant>
 
202
            <real_value shape="3" dim1="dim" rank="1">0 0 0</real_value>
 
203
          </constant>
 
204
        </value>
 
205
        <output/>
 
206
        <stat>
 
207
          <include_in_stat/>
 
208
        </stat>
 
209
        <detectors>
 
210
          <exclude_from_detectors/>
 
211
        </detectors>
 
212
      </prescribed>
 
213
    </vector_field>
 
214
  </material_phase>
 
215
</fluidity_options>