GfsOutputScalarStats

From Gerris

GfsOutputScalarStats is used to write the volume-weighted statistics over the whole domain of a given scalar.

The statistics are written using the following formatting:

DESCRIPTION time: T min: MIN avg: AVG | STDEV max: MAX

with:

DESCRIPTION

a description of the scalar field (without any spaces),

T

the physical time,

MIN

the minimum value,

AVG

the volume-weighted average value,

STDEV

the volume-weighted standard deviation,

MAX

the maximum value.

The syntax in parameter files is:

[ GfsOutputScalar ]

Examples

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

    OutputScalarStats { istep = 1 } r {
	v = (T > 1e-2 && T < 1. - 1e-2 ? 
	    (sqrt((Y + 0.5)*(Y + 0.5) + (Z + 0.5)*(Z + 0.5))/RADIUS - 1.)/EPSILON : 0)
    }

    OutputScalarStats { istep = 1 } k { v = K }

    OutputScalarStats { istep = 1 } kmax { v = Kmax }

• Forced isotropic turbulence in a triply-periodic box

  OutputScalarStats { istep = 1 } log { v = Unbar }

  OutputScalarStats { istep = 1 } log { v = U }

  OutputScalarStats { istep = 1 } Reynolds.dat { 
      v = 2./3.*FluctKinEn/VOLUME/MU*sqrt(15*MU/(Dissipation/VOLUME))
  }

  OutputScalarStats { istep = 1 } Dissipation.dat { v = Dissipation/VOLUME }

  OutputScalarStats { istep = 1 } Energy.dat { v = FluctKinEn/VOLUME }

  OutputScalarStats { istep = 1 } Vorticity.dat { v = Vorticity }  

• Tsunami runup onto a complex three-dimensional beach

    OutputScalarStats { istep = 1 } p { v = (Zb > 0. ? P : P + Zb) }

• The 2004 Indian Ocean tsunami

    OutputScalarStats { istep = 1 } p { v = P }

• "Garden sprinkler effect" in wave model

    OutputScalarStats { step = 12 } hs-MINLEVEL-NTHETA { v = Hs }

• Cyclone-generated wave field

    OutputScalarStats { step = 0.25 } hs { v = Hs }

    OutputScalarStats { step = 0.25 } vr { v = sqrt(U10*U10 + V10*V10) }

• Conservation of diffusive tracer

    OutputScalarStats { istep = 1 } t { v = T }

    OutputScalarStats { istep = 1 } te { v = T }

• Coriolis formulation in 3-D

  OutputScalarStats { step = 3000 } { awk '{print $3, $7}' > u.dat } { v = U }

  OutputScalarStats { step = 3000 } { awk '{print $3, $7}' > v.dat } { v = V }

  OutputScalarStats { step = 3000 } { awk '{print $3, $7}' > w.dat } { v = W }

  OutputScalarStats { istep = 10 } { awk '{print $3,$7}' > error.dat } {
      v = sqrt((U - Usol0)*(U - Usol0) + (V - Vsol0)*(V - Vsol0) + (W - Wsol0)*(W - Wsol0))
  }

• Wind-driven lake

    OutputScalarStats { istep = 1 } p { v = P }

• Circular droplet in equilibrium

  OutputScalarStats { istep = 1 } {
    awk '{ print MU*$3/(0.8*0.8), $5, $7, $9, $11 }' > K-LAPLACE-LEVEL
  } { v = (K - 2.50771) }

• Axisymmetric spherical droplet in equilibrium

  OutputScalarStats { istep = 1 } {
    awk '{ print MU*$3/(0.8*0.8), $5, $7, $9, $11; fflush (stdout); }' > K-LAPLACE-LEVEL
  } { v = (K - 2.50771) }

• Scalings for Plateau--Rayleigh pinchoff

   OutputScalarStats { istep = 1 } p { v = P }

   OutputScalarStats { istep = 1 } y { v = Y format = "%10.6e" }

   OutputScalarStats { istep = 1 } k { v = K }

• Sessile drop

    OutputScalarStats { istep = 10 } k { v = (T > 0.05 && T < 0.95 ? K : NODATA) }

• Dielectric-dieletric planar balance

    OutputScalarStats { start = end } p { v = P }

    OutputScalarStats { start = end } ey { v = Ey }

• Balance with solid boundaries

    OutputScalarStats { start = end } p { v = P }

    OutputScalarStats { start = end } ey { v = Ey }