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

<?xml version="1.0" encoding="UTF-8" ?>

 

<testproblem>

  <name>tracer_inflow_2d_dg_cg</name>

  <owner userid="dham"/>

  <tags>flml</tags>

  <problem_definition length="short" nprocs="1">

    <command_line>fluidity  tracer.flml </command_line>

    <!-- Two dimensional problem in which a pulse of tracer us introduced via a dirichlet boundary condition on the tracer. The tracer then flows out of the domain. This version tests the DG code path on a continuous velocity space.-->

  </problem_definition>

  <variables>

    <variable name="solvers_converged" language="python">

import os

files = os.listdir("./")

solvers_converged = not "matrixdump" in files and not "matrixdump.info" in files

    </variable>

    <variable name="tracer_integral_max" language="python">

from fluidity_tools import stat_parser

s = stat_parser("tracer.stat")

tracer_integral_max = max(s["Water"]["Tracer"]["integral"])

    </variable>

    <variable name="tracer_integral_final" language="python">

from fluidity_tools import stat_parser

s = stat_parser("tracer.stat")

tracer_integral_final = s["Water"]["Tracer"]["integral"][-1]

    </variable>

  </variables>

  <pass_tests>

    <test name="Solvers converged" language="python">

      assert(solvers_converged)

    </test>

    <test name="maximum tracer integral is 20/pi" language="python">

from math import pi

assert abs(tracer_integral_max - 20/pi) &lt; 1.5e-2

    </test>

    <test name="final tracer integral is 0" language="python">

assert abs(tracer_integral_final) &lt; 1.0e-2

    </test>

  </pass_tests>

  <warn_tests>

  </warn_tests>

</testproblem>