GfsOutputScalarStats

From Gerris

Jump to: navigation, search

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 }
    

Personal tools
communication