src/test/bore1.c

    Simple Saint-Venant Riemann problem

    This test is used by Le Métayer et al, 2010 (section 6.1.1) to provide the base solution for comparison with the Green-Naghdi equations.

    See also

    #include "grid/cartesian1D.h"
#include "saint-venant.h"

int main()
{
  X0 = -300.;
  L0 = 600.;
  G = 10.;
  N = 2048;
  run();
}

event init (i = 0)
{
  foreach()
    h[] = x < 0. ? 1.8 : 1.;
}

event output (t = 48) {
  foreach()
    fprintf (stdout, "%g %g %g\n", x, h[], u.x[]);
  fprintf (stdout, "\n");
}