[cig-commits] commit: Make hydrostatic use step() instead of if()

Mercurial hg at geodynamics.org
Sun May 15 11:02:03 PDT 2011


changeset:   348:5ca06f9a13d7
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Sun May 15 11:00:34 2011 -0700
files:       input/cookbook/hydrostatic.xml
description:
Make hydrostatic use step() instead of if()


diff -r f308d0928bb9 -r 5ca06f9a13d7 input/cookbook/hydrostatic.xml
--- a/input/cookbook/hydrostatic.xml	Sun May 15 11:00:14 2011 -0700
+++ b/input/cookbook/hydrostatic.xml	Sun May 15 11:00:34 2011 -0700
@@ -350,8 +350,8 @@
     </struct>
     <struct name="hydrostaticTerm">
       <param name="Type">HydrostaticTerm</param>
-      <param name="densityEquation">if(y &lt; .35, 1,0)</param>
-      <param name="pressureEquation">if(y &lt; .35, .35-y,0)</param>
+      <param name="densityEquation">step(.35-y)</param>
+      <param name="pressureEquation">step(.35-y)*(.35-y)</param>
     </struct>
     <struct name="stressBC">
       <param name="Type">StressBC</param>



More information about the CIG-COMMITS mailing list