1 0 GfsRiver GfsBox GfsGEdge { x = 0.5 } {
PhysicalParams { L = 10 g = 9.81 }
Refine 5
RefineSolid 7
Solid (sphere(5.,0.,0.,2.0))
Init {} {
Zb = x/10.
P = 12.1/(100. - 4.*M_PI)
}
Time { end = 200 }
SourceCoriolis 0 1.0e-01
OutputScalarNorm { start = end } u { v = U }
OutputErrorNorm { start = end } ep { v = P } {
s = MAX(0, 0.5 - x/10.)
unbiased = 1
relative = 1
}
GModule gfsview
OutputView { start = end } still.eps { format = EPS } still.gfv
EventScript { start = end } {
status=0
if awk '{if ($9 > 1e-5) { print "u: " $9 > "/dev/stderr"; exit (1); }}' < u ; then :
else
status=$GFS_STOP;
fi
if awk '{if ($9 > 4e-3) { print "ep: " $9 > "/dev/stderr"; exit (1); }}' < ep ; then :
else
status=$GFS_STOP;
fi
exit $status
}
}
GfsBox {
left = Boundary
right = Boundary
top = Boundary
bottom = Boundary
}