#include "grid/multigrid.h"
#include "compressible.h"

double tend = 3.;

int main() {

  foreach_dimension() {
    w.n[right] = neumann(0);
    w.n[left]  = neumann(0);
  }
  
  size(20);
  X0 = -L0/2.;
  N = 128;
  run(); 
}

event init (t = 0)
{ 
  
  double cson = sqrt(gammao);

  foreach() {
    double p = 1.+ 1.e-3*exp(-x*x);
    rho[] = 1. + (p - 1.)/sq(cson);
    foreach_dimension()
      w.x[] = 0.;
    E[] = p/(gammao - 1.);
  }
}

event endprint (t = tend) {
  
  foreach () {
    double xref = fabs(x) - tend*sqrt(gammao);
    printf("%g %g \n", xref, (E[]*(gammao-1)-1.)/1.e-3);
  }
  printf("\n");

}