1 0 GfsOcean GfsBox GfsGEdge {} {
  # dt = 1000 s
  Time { iend = 1580 dtmax = 0.10285 }
  Refine 6
  # Lx = Ly = 1000 km
  # H0 = 1000 m
  # g = 0.01 m/s^2
  PhysicalParams { g = 9.4534734306584e-4 }
  ApproxProjectionParams { tolerance = 1e-6 }
  Solid (z + 1.)
  Init {} {
    # e-folding radius = 100 km
    # umax = 0.5 m/s = sqrt(200)*exp(-1/2)
    U = (5.667583815e-4*200.*y*exp (-100.*(x*x + y*y)))
    V = (- 5.667583815e-4*200.*x*exp (-100.*(x*x + y*y)))
    P = (5.667583815e-4*exp (-100.*(x*x + y*y)))
  }
  # f0 = 1.0285e-4 s-1
  # beta = 1.607e-11 m-1s-1
  SourceCoriolis (1. + 0.156246961595*(y + 0.5))

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