GfsOutputScalarSum

From Gerris

GfsOutputScalarSum is used to write the sum over the whole domain of a given scalar.

The sum is written using the following formatting:

Description time: 1.2 sum: 2.578

where Description is a description of the scalar field.

The syntax in parameter files is:

[ GfsOutputScalar ]

By default the sum is weighted by the volume of each cell (i.e. it is the volume integral of the given scalar). This can be changed using the weight option of GfsOutputScalar.

Examples

• Collapse of a column of grains

    OutputScalarSum { istep = 10 } t-LEVEL { v = T }

    OutputScalarSum { istep = 10 } {
	awk '{
          print $3,$5/(H0*R0);
          fflush (stdout);
        }' > xg-H0-LEVEL
    } { v = T*x }

    OutputScalarSum { istep = 10 } {
	awk '{
          print $3,$5/(H0*R0);
          fflush (stdout);
        }' > yg-H0-LEVEL
    } { v = T*y }

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

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

• Dam break on complex topography

    OutputScalarSum { istep = 10 } ke { v = (P > 0. ? U*U/P : 0.) }

    OutputScalarSum { istep = 10 } vol { v = P }

• Tsunami runup onto a complex three-dimensional beach

    OutputScalarSum { istep = 1 } ke { v = (P > 0. ? Velocity2*P : 0.) }

• Convergence of the Godunov advection scheme

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

• Time-reversed VOF advection in a shear flow

    OutputScalarSum { istep = 1 } { awk '{ print $3,$5-8.743441e-01 }' > t } { v = T }

• Time-reversed advection of a VOF concentration

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

    OutputScalarSum { istep = 1 } t1 { v = C }

    OutputScalarSum { istep = 1 } t2 { v = G }

• Time-reversed advection with curvature-based refinement

    OutputScalarSum { istep = 1 } { awk '{ print $3,$5-0.8743665 }' > t } { v = T }

• Rotation of a straight interface

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

• Comparison between explicit and implicit diffusion schemes on concentration tracer

    OutputScalarSum {istep = 10}  T1vol { v = T1 }

    OutputScalarSum {istep = 10}  Cvol { v = C }

• Conservation of diffusive tracer

    OutputScalarSum { istep = 1 } st { v = T }

    OutputScalarSum { istep = 1 } ste { v = T }

• Estimation of the numerical viscosity

  OutputScalarSum { istep = 1 } kineticLEVEL { v = Velocity2 }

  OutputScalarSum { istep = 1 } stdout { v = Velocity2 }

• Estimation of the numerical viscosity with refined box

  OutputScalarSum { istep = 1 } kineticLEVEL { v = Velocity2 }

  OutputScalarSum { istep = 1 } stdout { v = Velocity2 }

• Numerical viscosity for the skew-symmetric scheme

  OutputScalarSum { istep = 1 } kineticLEVEL { v = Velocity2 }

• Numerical viscosity for the skew-symmetric scheme with refined box

  OutputScalarSum { istep = 1 } kinetic-LEVEL { v = Velocity2 }

• Numerical viscosity for vorticity/streamfunction formulation

  OutputScalarSum { istep = 1 } kineticLEVEL { v = Velocity2 }

  OutputScalarSum { istep = 1 } stdout { v = Velocity2 }

• Mass conservation

    OutputScalarSum { istep = 1 end = 0.8 } srt { v = T }

    OutputScalarSum { istep = 1 end = 0.8 } srt1 { v = T1 }

• Mass conservation with solid boundary

    OutputScalarSum { istep = 1 } srt { v = T }

    OutputScalarSum { istep = 1 } srt1 { v = T1 }

• Momentum conservation for large density ratios

    OutputScalarSum { istep = 1 } k { v = Velocity2*rho(T1) }

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

• Translation of an hexagon in a uniform flow

  OutputScalarSum { istep = 1 } {
      awk '{ printf ("%e %e\n", $3, $5 - 1.953125) }' > tracersum-ORDER
  } { v = T }

• Shape oscillation of an inviscid droplet

    OutputScalarSum { istep = 1 } k-LEVEL {
        v = RHO(T1)*Velocity2
    }

• Scalings for Plateau--Rayleigh pinchoff

   OutputScalarSum { istep = 1 } ke { v = Velocity2*T }

• Sessile drop

    OutputScalarSum { start = end } vol { v = T }

• Geostrophic adjustment on a beta-plane

  OutputScalarSum { istep = 150 } { 
      awk '{print $3/1.0285e-4/3600./24. " " $5/9.683940e-11}' > energy 
  } { v = (Velocity2 + P*P/9.4534734306584e-4) }

• Non-linear geostrophic adjustment

  OutputScalarSum { istep = 1 } energy-OMEGA { v = Velocity2 }

• Gravity waves in a realistic ocean basin

    OutputScalarSum { istep = 10 } k { v = Velocity2 }

• Oscillations in a parabolic container

    OutputScalarSum { istep = 10 } ke-LEVEL { v = (P > 0. ? U*U/P : 0.) }

    OutputScalarSum { step = 50 } vol-LEVEL { v = P }

    OutputScalarSum { step = 50 } U-LEVEL { v = U }

• Parabolic container with embedded solid

    OutputScalarSum { istep = 10 } ke-LEVEL { v = (P > 0. ? U*U/P : 0.) }

    OutputScalarSum { step = 50 } vol-LEVEL { v = P }

    OutputScalarSum { step = 50 } U-LEVEL { v = U*cos (ANGLE) + V*sin(ANGLE) }

• Lake-at-rest balance in an inclined domain with bipolar metric

    OutputScalarSum { start = 0 } vol { v = Zb }

• Circular dam break on a sphere

    OutputScalarSum { istep = 1 } sp { v = P }

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

    OutputScalarSum { istep = 1 } sp { v = P }

• Advection of a cosine bell around the sphere

  OutputScalarSum { istep = 1 } t-LEVEL-ALPHA { v = T }

  OutputScalarSum { istep = 1 } area-LEVEL-ALPHA { v = 1 }

• Rossby--Haurwitz wave

    OutputScalarSum { istep = 10 } ec-LEVEL { v = Velocity2 }

    OutputScalarSum { istep = 10 } zeta-LEVEL { v = Vorticity }

    OutputScalarSum { istep = 10 } p-LEVEL { v = P }

• Rossby--Haurwitz wave with a free surface

    OutputScalarSum { istep = 10 } ec-LEVEL { v = Velocity2 }

    OutputScalarSum { istep = 10 } zeta-LEVEL { v = Vorticity }

    OutputScalarSum { istep = 10 } p-LEVEL { v = P }

• Rossby--Haurwitz wave with Saint-Venant

    OutputScalarSum { istep = 10 } ec-LEVEL { v = Velocity2 }

    OutputScalarSum { istep = 10 } p-LEVEL { v = P }

• Charge relaxation in an axisymmetric insulated conducting column

    OutputScalarSum { step = 1 } {
	awk 'BEGIN { R0 = 0.1 ; rhoinic = 0.5 ; L =1.0 ; Q = 0.5*R0*R0*L*rhoinic }
	     { print $3,$5,100*sqrt((1.0 - $5/Q)*(1.0 - $5/Q)) }' > rhoe-LEVEL 
    } { v = Rhoe }

    OutputScalarSum { step = 1 } {
	awk 'BEGIN { R0 = 0.1 ; rhoinic = 0.5 ; L =1.0 ; Q = 0.5*R0*R0*L*rhoinic }
	     { print $3,$5,100*sqrt((1.0 - $5/Q)*(1.0 - $5/Q)) }' > C-LEVEL
    } { v = C }

• Charge relaxation in a planar cross-section

    OutputScalarSum { istep = 1 } { 
       awk 'BEGIN { rhoinic = 0.5 ; R0 = 0.1 ; Q = rhoinic*R0*R0*3.141592654 }
            { print $3,$5,100*sqrt((1.-$5/Q)*(1.-$5/Q)); fflush(stdout); }' > rhoe-LEVEL
    } { v = Rhoe}

• Equilibrium of a droplet suspended in an electric field

    OutputScalarSum { istep = 10 } rhoe { v = Rhoe }