[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