#include "navier-stokes/centered.h"
#define U_GEO 1.
#include "force_geo.h"
scalar * tracers;      //We need to declare the tracers

uz[bottom] = dirichlet (0.);
u.t[bottom] = dirichlet (0.);

FILE * fp;

int main() {
  fp = fopen ("vy", "w");
  periodic (left);
  f_cor = 1.;
  mu = unityf;
  L0 = 25.;
  for (TOLERANCE = 1e-2; TOLERANCE >= 1e-6; TOLERANCE /= 10.)
    run();
}

event init (t = 0) {
  foreach()
    u.x[] = U_GEO;
}

event adapt (i++)
  adapt_wavelet ({u.x, uz}, (double[]){0.005, 0.005}, 8);

#include "profile5c.h"
scalar u_test[];
event end (t < HUGE; i += 10) {
  if (change (u.x, u_test) < 2e-3) {
    fprintf (fp, "%g %g\n", TOLERANCE, normf(u.y).avg);
    if (fabs (TOLERANCE - 1e-3) < 1e-4)
      profile ({u.x, uz}, fname = "prof");
    return 1;
  }
}