GfsOutputSimulation

From Gerris

GfsOutputSimulation will write a description of the current state of the simulation which contains both standard simulation parameters, layout of the cell hierarchy and associated variable values.

The syntax in parameter files is:

[ GfsOutput ] {
  depth = 5 
  variables = U,V,C 
  binary = 1 
  solid = 0 
  format = text
  precision = %.8f
}

with:

depth

an optional parameter setting the maximum depth which will be considered (the default is infinity).

variables

a comma-separated list of the names of the variables which will be written. The default is all the simulation permanent variables.

binary

an optional parameter specifying if a binary format should be used. Binary files allow faster reading of simulation files and are the default. Setting binary to 0 will use a more portable ASCII format.

solid

an optional parameter specifying if the surface describing solid boundaries should be included in the simulation file. The default is to include the solid surface.

format

an optional parameter specifying which file format to use:

• gfs: the standard Gerris format. This is the default.
• text: an ASCII file containing only cell coordinates and variable values. The first line of the text-formatted file contains a description of the different fields (column by column).
• VTK: the VTK file format.
• Tecplot: the Tecplot file format.

precision

an optional C formatting string used to convert floating-point numbers to their ASCII representation (only relevant for ASCII file formats). The default is %g.

Examples

• Parallel simulation on four processors

  OutputSimulation { start = end } end.gfs

• Rayleigh-Taylor instability

  OutputSimulation { step = 0.1 } stdout

• Boussinesq flow generated by a heated cylinder

  OutputSimulation { istep = 4 } stdout

• Coalescence of a pair of Gaussian vortices (Gerris logo)

    OutputSimulation { istep = 10 } stdout

• Collapse of a column of grains

    OutputSimulation { istep = 10 } stdout

    OutputSimulation { start = 0 } snapshot-%g.gfs

    OutputSimulation { start = 1.65132 } snapshot-%g.gfs

    OutputSimulation { start = 2.3769 } snapshot-%g.gfs

    OutputSimulation { start = 3.10248 } snapshot-%g.gfs

    OutputSimulation { start = 3.80304 } snapshot-%g.gfs

    OutputSimulation { start = 5.70456 } snapshot-%g.gfs

• Starting vortex of a NACA 2414 aerofoil

  OutputSimulation { step = FRAMEPERIOD } stdout

• Turbulent air flow around RV Tangaroa

  OutputSimulation { istep = 4 } stdout

  OutputSimulation { start = end } simulation-sum {
      variables = SU,SV,SW,SU2,SV2,SW2
  }

• Savart--Plateau--Rayleigh instability of a water column

    OutputSimulation { istep = 1000 } plateau-%ld.gfs

• Atomisation of a pulsed liquid jet

    OutputSimulation { istep = 100 } snapshot-%ld.gfs { }

• Air-water flow around a Series 60 cargo ship

    OutputSimulation { istep = 5 end = 4 } stdout

    OutputSimulation { start = 1 step = 1 } sim-%g.gfs

• Forced isotropic turbulence in a triply-periodic box

  OutputSimulation { istep = 100 } snapshot-%ld.gfs

  OutputSimulation { start = end } end.gfs

• Wingtip vortices behind a rectangular NACA 2414 wing

  OutputSimulation { step = FRAMEPERIOD } stdout

• Lunar tides in Cook Strait, New Zealand

    OutputSimulation { istep = 20 } stdout { solid = 0 }

    OutputSimulation { start = end } end.gfs

• Dam break on complex topography

    OutputSimulation { step = 0.3 } {
	gfsplot "
          load 'dam.plot'
          set title sprintf('t = %4.1f', (t))
          plot [-4.:4.]'-' u (x):(Zb):(H) w filledcu lc 3, \
                       '-' u (x):(Zb) w l lw 4 lc 1 lt 1
          set term pngcairo size 800,600
          set output sprintf('sim-%04.1f.png', (t))
          replot
          set term wxt noraise
        "
    } { format = text }

• Small amplitude solitary wave interacting with a parabolic hump

    OutputSimulation { istep = 10 } stdout

    OutputSimulation { start = 0.6 step = 0.3 } sim-%g.gfs

• Shock reflection by a circular cylinder

    OutputSimulation { istep = 10 } stdout

• Tsunami runup onto a complex three-dimensional beach

    OutputSimulation { step = 1 } stdout

    OutputSimulation { step = 1 } sim-%g.gfs

• The 2004 Indian Ocean tsunami

    OutputSimulation { istep = 100 } snapshot-%ld.gfs

    OutputSimulation { istep = 20 } stdout

    OutputSimulation { step = 900 } sim-%g.gfs

• "Garden sprinkler effect" in wave model

    OutputSimulation { step = 12 } sim-MINLEVEL-NTHETA-%g.gfs

    OutputSimulation { start = end } end-MINLEVEL-NTHETA.gfs    

• Cyclone-generated wave field

    OutputSimulation { istep = 10 } stdout

    OutputSimulation { step = 4 } sim-%g.gfs

• Convergence of the Poisson solver

  OutputSimulation { start = end } end-SOLVER.gfs

• Convergence with a refined circle

  OutputSimulation { start = end } end-SOLVER.gfs

• Dirichlet boundary condition

    OutputSimulation { start = end } end-SOLVER.gfs

• Convergence of the Poisson solver with solid boundaries

  OutputSimulation { start = end } sim-LEVEL-SOLVER { variables = P }

• Star-shaped solid boundary with refinement

    OutputSimulation { start = end } sim-LEVEL-SOLVER { variables = P }

• Star-shaped solid boundary

    OutputSimulation { start = end } sim-LEVEL-SOLVER { variables = P }

• Thin wall at box boundary

  OutputSimulation { start = end } sim-LEVEL-SOLVER { variables = P }

• Time-reversed advection of a VOF concentration

    OutputSimulation { start = 2.5 } half-LEVEL.gfs    

    OutputSimulation { start = end } end-LEVEL.gfs    

• Comparison between the explicit and implicit diffusion schemes

    OutputSimulation { start = 0.01 } end.gfs

    OutputSimulation { start = 2.5e-3 step = 2.5e-3 } {
	awk '{
           if ($1 != "#")
             print sqrt($1*$1+$2*$2),$8;
        }' > prof
    } { format = text }

• Comparison between explicit and implicit diffusion schemes on concentration tracer

    OutputSimulation { start = 0.01 } end.gfs

    OutputSimulation {start = 0.005} {
	awk '{
           if ($1 != "#" && $8 > 0.)
             print sqrt($1*$1+$2*$2),$12,$13;
        }' > prof
    } { format = text}

• Convergence for the three-way vortex merging problem

  OutputSimulation { start = 0.05 } stdout

  OutputSimulation { start = 0.15 } stdout

  OutputSimulation { start = 0.25 } stdout

  OutputSimulation { start = 0.25 } SIM

• Flow created by a cylindrical volume source

    OutputSimulation { start = end } end-LEVEL.gfs

• Potential flow around a sphere

    OutputSimulation { start = end } sim-LEVEL.gfs

• Viscous flow past a sphere

    OutputSimulation { start = end } end-LEVEL-RE.gfs

• Mass conservation

    OutputSimulation { step = 0.2 } stdout

• Boundary layer on a rotating disk

    OutputSimulation { start = end } end.txt { format = text variables = U,W }

    OutputSimulation { start = end } end.gfs

• Lid-driven cavity at Re=1000

  OutputSimulation { start = end } end.gfs

• Lid-driven cavity with a non-uniform metric

  OutputSimulation { start = end } end.gfs

• Lid-driven cavity on an anisotropic mesh

  OutputSimulation { start = end } end.gfs

• Bagnold flow of a granular material

    OutputSimulation { start = end } end-LEVEL.txt { format = text }

• Poiseuille flow with multilayer Saint-Venant

    OutputSimulation { start = end } {
	awk '{
  	  nl = NL;
          if ($1 == "#") {
            for (i = 2; i <= NF; i++) {
              split($i,a,":")
              if (a[2] == "U0")
                start = a[1];
            }
          }
          else {
            dz = $4/nl;
            emax = 0.;
  	    for (i = 0; i < nl; i++) {
              z = dz*(0.5+i)
              u = A/(2.*NU)*(1./4 - (0.5 - z)*(0.5 - z))
              eu = u - $(start+i)/dz
              if (eu < 0.)
                eu = -eu;
              if (eu > emax)
                emax = eu;
            }
            print nl,emax
          }
        }'
    } { format = text }

• Coriolis formulation in 3-D

  OutputSimulation { start = end } end.gfs

• Wind-driven lake

    OutputSimulation { start = end } end.gfs

• Multi-layer Saint-Venant solver

    OutputSimulation { start = end } {
	awk '{
  	  nl = NL;
          if ($1 == "#") {
            for (i = 2; i <= NF; i++) {
              split($i,a,":")
              if (a[2] == "U0")
                start = a[1];
            }
          }
          else if ($1 == RATIO/2**(LEVEL + 1)) {
            dz = $4/nl;
  	    for (i = 0; i < nl; i++) {
              z = dz*(0.5+i)
              print $1,$2,z,$(start+i)/dz
            }
          }
        }' > uprof-RATIO-RE-NL
    } { format = text }

• Wind-driven stratified lake

    OutputSimulation { start = end } end-NL.txt { format = text }

    OutputSimulation { start = end } {
	awk '{
  	  nl = NL;
          if ($1 == "#") {
            for (i = 2; i <= NF; i++) {
              split($i,a,":")
              if (a[2] == "U0")
                start = a[1];
            }
          }
          else if ($1 == RATIO/2**(LEVEL + 1)) {
            dz = $4/nl;
  	    for (i = 0; i < nl; i++) {
              z = dz*(0.5+i)
              print $1,$2,z,$(start+i)/dz
            }
          }
        }' > uprof-NL
    } { format = text }

• Convergence of a potential flow solution

    OutputSimulation { start = end } sim-LEVEL {
        variables = U,V,P
    }

• Flow through a divergent channel

    OutputSimulation { start = end } sim-LEVEL {
        variables = U,V,P
    }

• Translation of an hexagon in a uniform flow

  OutputSimulation { start = end } end-ORDER.gfs

• B\'enard--von K\'arm\'an vortex street behind a cylinder translating in a fluid at rest

  OutputSimulation { start = end } end.gfs 

• Circular droplet in equilibrium

  OutputSimulation { start = end } stdout { depth = LEVEL }

• Axisymmetric spherical droplet in equilibrium

  OutputSimulation { start = end } stdout { depth = LEVEL }

• Height-Function on parallel subdomains

    OutputSimulation { start = end } stdout

• Sessile drop

    OutputSimulation { start = end } end-ANGLE.gfs

• Non-linear geostrophic adjustment

  OutputSimulation { start = end } end-OMEGA.gfs

• Coastally-trapped waves

  OutputSimulation { start = end } sim-LEVEL

• Coastally-trapped waves with adaptive refinement

  OutputSimulation { start = end } sim-LEVEL

• Oscillations in a parabolic container

    OutputSimulation { start = 1500 } sim-LEVEL-%g.txt { format = text }

    OutputSimulation { start = end } end-LEVEL.txt { format = text }

• Parabolic container with embedded solid

    OutputSimulation { start = 1500 } { 
	awk '
        function atan(x)
        {
          return atan2(x,1.);
        }
        function pow(x,y)
        {
          return x**y;
        }
        {
          if ($1 == "#")
            print $0;
          else {
            printf ("%g", R($1,$2));
            for (i = 2; i <= NF; i++)
              printf (" %s", $i);
            printf ("\n");
          }
        }' > sim-LEVEL-1500.txt 
    } { format = text }

    OutputSimulation { start = end } end-LEVEL.txt { format = text }

• Transcritical flow over a bump

    OutputSimulation { start = end } {
	gfsplot "
          set term postscript eps color solid lw 2 16
          set output 'h-LEVEL.eps'
          set xlabel 'x'
          set ylabel 'z'
          plot [5:15]'swashes' u 1:4 t 'topography' w l, \
                     'swashes' u 1:6 t 'free surface (analytical)' w l, \
                     '-' u (x):(H) w p t 'free surface (numerical)'
        "
    } { format = text }

    OutputSimulation { start = end } end-LEVEL.txt { format = text }

• Transcritical flow with multiple layers

    OutputSimulation { start = end } {
	awk '{ if ($1 != "#") print $1,$8; }' > prof-LEVEL-NL
    } { format = text }

    OutputSimulation { start = end } end-LEVEL-NL { format = text }

• Tsunami runup onto a plane beach

    OutputSimulation { start = 0 } sim-LEVEL-%g.txt { format = text }

    OutputSimulation { start = 160 } sim-LEVEL-%g.txt { format = text }

    OutputSimulation { start = 175 } sim-LEVEL-%g.txt { format = text }

    OutputSimulation { start = 220 } sim-LEVEL-%g.txt { format = text }

• Terrain reconstruction

    OutputSimulation { start = end } {
	awk '{
           if ($1 != "#") {
             r = sqrt ($1*$1 + $2*$2);
             print r,$4,$5,$6,$7;
           }
        }' | sort -n -k 1,2 > end-LEVEL.txt
    } { format = text variables = T,H0,Et,Eh }

    OutputSimulation { start = end } end-LEVEL.gfs

• Circular dam break on a sphere

    OutputSimulation { step = 0.3 } sim-%g.txt { variables = U,V,P format = text }

• Circular dam break on a rotating sphere

    OutputSimulation { step = 0.4 } sim-%g.txt { variables = U,V,P format = text }

• Circular dam break on a ``cubed sphere''

    OutputSimulation { step = 0.3 } sim-%g.txt { variables = U,V,P format = text }

• Advection of a cosine bell around the sphere

  OutputSimulation { start = end } end-LEVEL-ALPHA.gfs

• Poisson problem with a pure spherical harmonics solution

  OutputSimulation { start = end } end-SOLVER-LEVEL.gfs

• Spherical harmonics with longitude-latitude coordinates

  OutputSimulation { start = end } end-SOLVER-LEVEL.gfs

• Poisson equation on a sphere with Gaussian forcing

  OutputSimulation { start = end } end.gfs

• Gaussian forcing using longitude-latitude coordinates

  OutputSimulation { start = end } end.gfs

• Creeping Couette flow between cylinders

    OutputSimulation { start = end } end-LEVEL.gfs

    OutputSimulation { start = end } end-LEVEL.txt { format = text }

• Creeping Couette flow between eccentric cylinders

  OutputSimulation { start = end } end-LEVEL.gfs

• Flow between eccentric cylinders using bipolar coordinates

  OutputSimulation { start = end } end-LEVEL.gfs

• Flow between eccentric cylinders on a stretched grid

  OutputSimulation { start = end } end-LEVEL.gfs

• Rossby--Haurwitz wave

    OutputSimulation { start = end } end-LEVEL.gfs

• Rossby--Haurwitz wave with a free surface

    OutputSimulation { start = end } end-LEVEL.gfs

• Rossby--Haurwitz wave with Saint-Venant

    OutputSimulation { start = end } end-LEVEL.gfs

• Dielectric-dieletric planar balance

    OutputSimulation { start = end } result.txt { format = text }

• Balance with solid boundaries

    OutputSimulation { start = end } result.txt { format = text }

• Relaxation of a charge bump

    OutputSimulation { start = 0 } {
	awk '{ if ($1 != "#") print sqrt($1*$1+$2*$2),$11 > "t_0-SCHEME"; }'
    } { format = text }

    OutputSimulation { start = 2 } {
	awk '{ if ($1 != "#") print sqrt($1*$1+$2*$2),$11 > "t_2-SCHEME"; }'
    } { format = text }

    OutputSimulation { start = 4 } {
	awk '{ if ($1 != "#") print sqrt($1*$1+$2*$2),$11 > "t_4-SCHEME"; }'
    } { format = text }

    OutputSimulation { start = 6 } {
	awk '{ if ($1 != "#") print sqrt($1*$1+$2*$2),$11 > "t_6-SCHEME"; }'
    } { format = text }

• Charge relaxation in an axisymmetric insulated conducting column

     OutputSimulation { start = end } {
	awk '{ if ($1 != "#") print $2,sqrt($4*$4+$5*$5); }' > prof-LEVEL
    } { 
	format = text 
	variables = Ex,Ey 
    }

    OutputSimulation { start = end } result-LEVEL.gfs 

• Charge relaxation in a planar cross-section

    OutputSimulation { start = end } {
	awk '{ if ($1 != "#") 
                 print sqrt($2*$2+$1*$1),sqrt($4*$4+$5*$5); }' > prof-LEVEL
    } { 
	format = text 
	variables = Ex,Ey 
    }

    OutputSimulation { start = end } result-LEVEL.gfs 

• Equilibrium of a droplet suspended in an electric field

    OutputSimulation { start = end } result.gfs

• Simple example of groundwater flow following Darcy's law

	OutputSimulation solution.gfs

• Groundwater flow with piecewise constant permeability

	OutputSimulation {} solution.gfs