GfsOutputScalarSum

From Gerris

GfsOutputScalarSum is used to write the volume-weighted 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 ]

Examples

• Dam break on complex topography

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

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

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

• Mass conservation

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

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

• Mass conservation with solid boundary

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

    OutputScalarSum { istep = 1 } srt1 { v = y*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
    }

• Geostrophic adjustment on a beta-plane

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

• 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 { istep = 10 } vol-LEVEL { v = P }

    OutputScalarSum { istep = 10 } U-LEVEL { v = U }