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.

Note

There is no constraint on the coordinates of the input locations but only locations which are within the simulation domain will be written in the output file.

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