← Back to branch summary

~fluidity-core/fluidity/sediment_restructure

« back to all changes in this revision

Viewing changes to tests/test_advection_diffusion/cdisk_diff.adml

merge with branch including remap_surface_to_field but then removed this change and approached it a different way which seems to work

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/test_advection_diffusion/cdisk_diff.adml
1
 
<?xml version='1.0' encoding='utf-8'?>
2
 
<fluidity_options>
3
 
  <simulation_name>
4
 
    <string_value lines="1">cdisk_ad</string_value>
5
 
  </simulation_name>
6
 
  <geometry>
7
 
    <dimension>
8
 
      <integer_value rank="0">2</integer_value>
9
 
    </dimension>
10
 
    <mesh name="CoordinateMesh">
11
 
      <from_file file_name="cdisk">
12
 
        <format name="triangle">
13
 
          <string_value>triangle</string_value>
14
 
        </format>
15
 
        <stat>
16
 
          <include_in_stat/>
17
 
        </stat>
18
 
      </from_file>
19
 
    </mesh>
20
 
    <mesh name="VelocityMesh">
21
 
      <from_mesh>
22
 
        <mesh name="CoordinateMesh"/>
23
 
        <stat>
24
 
          <exclude_from_stat/>
25
 
        </stat>
26
 
      </from_mesh>
27
 
    </mesh>
28
 
    <quadrature>
29
 
      <degree>
30
 
        <integer_value rank="0">4</integer_value>
31
 
      </degree>
32
 
    </quadrature>
33
 
  </geometry>
34
 
  <io>
35
 
    <dump_format>
36
 
      <string_value>vtk</string_value>
37
 
    </dump_format>
38
 
    <dump_period_in_timesteps>
39
 
      <constant>
40
 
        <integer_value rank="0">1</integer_value>
41
 
      </constant>
42
 
    </dump_period_in_timesteps>
43
 
    <output_mesh name="VelocityMesh"/>
44
 
  </io>
45
 
  <timestepping>
46
 
    <current_time>
47
 
      <real_value rank="0">0.0</real_value>
48
 
    </current_time>
49
 
    <timestep>
50
 
      <real_value rank="0">0.0125663706</real_value>
51
 
      <comment>2 * pi/50 (50 steps per revolution)</comment>
52
 
    </timestep>
53
 
    <finish_time>
54
 
      <real_value rank="0">1000</real_value>
55
 
    </finish_time>
56
 
    <final_timestep>
57
 
      <integer_value rank="0">50</integer_value>
58
 
      <comment>Go around once.</comment>
59
 
    </final_timestep>
60
 
  </timestepping>
61
 
  <material_phase name="Fluid">
62
 
    <scalar_field name="Tracer" rank="0">
63
 
      <prognostic>
64
 
        <mesh name="VelocityMesh"/>
65
 
        <solver>
66
 
          <iterative_method name="cg"/>
67
 
          <preconditioner name="sor"/>
68
 
          <relative_error>
69
 
            <real_value rank="0">1.0e-7</real_value>
70
 
          </relative_error>
71
 
          <max_iterations>
72
 
            <integer_value rank="0">1000</integer_value>
73
 
          </max_iterations>
74
 
          <never_ignore_solver_failures/>
75
 
          <diagnostics>
76
 
            <monitors/>
77
 
          </diagnostics>
78
 
        </solver>
79
 
        <initial_condition name="WholeMesh">
80
 
          <python>
81
 
            <string_value lines="20" type="python">def val(X,t):
82
 
    from numpy import matrix
83
 
    from math import sqrt
84
 
    dx= (matrix(X)-matrix((-0.5,0)))
85
 
    r=sqrt(dx*dx.T)
86
 
    if (r&lt;0.25):
87
 
        return 1.0
88
 
    else:
89
 
        return 0</string_value>
90
 
          </python>
91
 
        </initial_condition>
92
 
        <tensor_field name="Diffusivity" rank="2">
93
 
          <prescribed>
94
 
            <value name="WholeMesh">
95
 
              <isotropic>
96
 
                <constant>
97
 
                  <real_value rank="0">0.001</real_value>
98
 
                </constant>
99
 
              </isotropic>
100
 
            </value>
101
 
          </prescribed>
102
 
        </tensor_field>
103
 
        <output/>
104
 
        <stat/>
105
 
      </prognostic>
106
 
    </scalar_field>
107
 
  </material_phase>
108
 
</fluidity_options>