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

Stephane Popinet popinet at users.sf.net
Fri May 15 02:55:54 UTC 2009


The following commit has been merged in the upstream branch:
commit d35ad2458debf20fc7f4bb3661cbf38513398009
Author: Stephane Popinet <popinet at users.sf.net>
Date:   Fri Nov 14 11:53:37 2008 +1100

    Added "none" option for gradient of AdvectionParams
    
    This means that the advection schemes can be made first-order in space.
    
    darcs-hash:20081114005337-d4795-0ee829ebe97039fc9f12a112c98fc0f2b2f3946b.gz

diff --git a/src/advection.c b/src/advection.c
index fdeb426..d765f51 100644
--- a/src/advection.c
+++ b/src/advection.c
@@ -912,7 +912,9 @@ void gfs_advection_params_write (GfsAdvectionParams * par, FILE * fp)
 	   par->cfl,
 	   par->gradient == gfs_center_gradient ? 
 	   "gfs_center_gradient" :
-	   "gfs_center_van_leer_gradient",
+	   par->gradient == gfs_center_van_leer_gradient ?
+	   "gfs_center_van_leer_gradient" :
+	   "none",
 	   par->flux == gfs_face_advection_flux ?
 	   "gfs_face_advection_flux" :
 	   par->flux == gfs_face_velocity_advection_flux ?
@@ -946,6 +948,11 @@ void gfs_advection_params_init (GfsAdvectionParams * par)
   par->gc = FALSE;
 }
 
+static gdouble none (FttCell * cell, FttComponent c, guint v)
+{
+  return 0.;
+}
+
 void gfs_advection_params_read (GfsAdvectionParams * par, GtsFile * fp)
 {
   GtsFileVariable var[] = {
@@ -980,6 +987,8 @@ void gfs_advection_params_read (GfsAdvectionParams * par, GtsFile * fp)
       par->gradient = gfs_center_gradient;
     else if (!strcmp (gradient, "gfs_center_van_leer_gradient"))
       par->gradient = gfs_center_van_leer_gradient;
+    else if (!strcmp (gradient, "none"))
+      par->gradient = none;
     else if (fp->type != GTS_ERROR)
       gts_file_variable_error (fp, var, "gradient",
 			       "unknown gradient parameter `%s'", gradient);

-- 
Gerris Flow Solver



More information about the debian-science-commits mailing list