[SCM] Gerris Flow Solver branch, upstream, updated. b3aa46814a06c9cb2912790b23916ffb44f1f203

[Stephane Popinet]

The following commit has been merged in the upstream branch:
commit 28e1b95d321700dd5e6e64f7f92002bd371dd130
Author: Stephane Popinet <popinet at users.sf.net>
Date:   Tue Mar 4 11:35:16 2008 +1100

    New quadratic hydrostatic pressure test case
    
    darcs-hash:20080304003516-d4795-8c266e3c3e8745c785e0809e67c941c6490acc74.gz

diff --git a/test/hydrostatic/quadratic/quadratic.gfs b/test/hydrostatic/quadratic/quadratic.gfs
new file mode 100644
index 0000000..782b59a
--- /dev/null
+++ b/test/hydrostatic/quadratic/quadratic.gfs
@@ -0,0 +1,46 @@
+# Title: Hydrostatic balance with quadratic pressure profile
+#
+# Description:
+#
+# Same test as before but for a quadratic pressure profile.
+#
+# Author: St\'ephane Popinet
+# Command: gerris2D quadratic.gfs
+# Version: 1.2.0
+# Required files: quadratic.gfs
+1 0 GfsSimulation GfsBox GfsGEdge {} {
+    Refine 3
+
+    # This test case only works for constant refinement
+    #    Refine (x*x + y*y < 0.2*0.2 ? 4 : 3) 
+
+    # Note: it is important to use 'cy' rather than 'y' in the formula
+    # below so that the hydrostatic density distribution is correct
+    # even for 'cut cells'
+    Init {} { rho = (cy + 0.5) }
+
+    Source V -rho
+    SourceViscosity 1e-2
+    Solid (ellipse(0.,0.,0.24,0.24))
+    Time { iend = 10 }
+    ApproxProjectionParams { tolerance = 1e-12 }
+    ProjectionParams { tolerance = 1e-12 }
+
+    OutputScalarNorm { istep = 1 } v { v = V }
+    # Checks that the pressure profile is close to the exact solution
+    OutputErrorNorm { istep = 1 } p { v = P } {
+	s = (cy*cy/2. + 0.5*cy) 
+	unbiased = 1 
+    }
+    EventScript { start = end } { 
+        if awk '{if ($9 > 1e-12) exit (1);}' < v ; then :
+        else
+            return $GFS_STOP;
+        fi        
+	if awk '{if ($9 > 1e-12) exit (1);}' < p ; then :
+        else
+            return $GFS_STOP;
+        fi
+    } 
+}
+GfsBox {}
diff --git a/test/template.tex b/test/template.tex
index fc4526b..b79d7d3 100644
--- a/test/template.tex
+++ b/test/template.tex
@@ -60,6 +60,7 @@ current stable branch of the version-controlled source code.
 \input{couette/couette.tex}
 \input{kinetic/kinetic.tex}
 \input{hydrostatic/hydrostatic.tex}
+\input{hydrostatic/quadratic/quadratic.tex}
 
 \section{Solid boundaries}
 

-- 
Gerris Flow Solver