GfsOutputLocation

From Gerris

GfsOutputLocation writes the values of all permanent variables at a set of given locations. The first time GfsOutputLocation is called it writes a comment describing the fields being written.

The values are linearly interpolated to the given location using gfs_interpolate().

The syntax in parameter files is as follows:

[ GfsOutput ] X Y Z

where X, Y and Z are the coordinates of the location.

or:

[ GfsOutput ] FILE

where FILE is the name of a text file containing the space-separated X, Y and Z coordinates of a set of locations.

or:

[ GfsOutput ] { X1 Y1 Z1 ... XN YN ZN }

where the Xi, Yi and Zi are the coordinates of the set of locations.

Examples

• Viscous flow past a sphere

    OutputLocation { step = 0.1 } {
	awk 'BEGIN { t = 2.; oldl = -1.; oldt = 0.; } {
          if ($1 != t) { t = $1; x1 = $2; u1 = $7; }
          else {
            x2 = $2; u2 = $7;
            if (u1 <= 0. && u2 > 0.) {
              l = (u1*x2 - u2*x1)/(u1 - u2) - A0;
              dl = (l - oldl)/(t - oldt);
              print t, l, dl;
              fflush (stdout);
              oldl = l;
              oldt = t;
            }
            x1 = x2; u1 = u2;
          }
        }' > l-LEVEL-RE
    } axis

• Lid-driven cavity at Re=1000

  OutputLocation { start = end } xprof xprofile

  OutputLocation { start = end } yprof yprofile

• Lid-driven cavity at Re=1000 (explicit scheme)

  OutputLocation { start = end } xprof ../xprofile

  OutputLocation { start = end } yprof ../yprofile

• Creeping Couette flow of Generalised Newtonian fluids

  OutputLocation { start = end } { awk '{if ($1 != "#") print $2,$8;}' > prof-MODEL } profile