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 w 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 }