6 12 GfsPoisson GfsBox GfsGEdge {} {
  PhysicalParams { L = 2.*M_PI/4. }
  MetricCubed M LEVEL
  Time { iend = 1 }
  Refine LEVEL

  ProjectionParams { tolerance = 1e-12 }
  ApproxProjectionParams { tolerance = 1e-12 }

  Init {} { Div = exp(-2.*5.*5.*(1.-cos((y + 90.)/180.*M_PI))) }

  OutputLocation { start = end } prof.dat profile

  OutputSimulation { start = end } end.gfs
  OutputProjectionStats { istep = 1 } stderr

  EventScript { start = end } {
      gnuplot <<EOF
      set term pos enhanced eps color solid 20 lw 3 
      set out 'profile.eps'

      set key bottom right
      set xl "Latitude"
      set yl "{/Symbol F}"
      set xr [-90:90]

      plot './prof.dat' u 3:5 every 5 ps 2 t "Gerris",\
           'prof.ref'u 1:2 w l t "(Boyd and Zhou, 2009)"
EOF
  }
  
}
GfsBox {}
GfsBox {}
GfsBox {}
GfsBox {}
GfsBox {}
GfsBox {}
1 2 right
2 3 top
3 4 right
4 5 top
5 6 right
6 1 top
1 3 top left
3 5 top left
5 1 top left
2 6 bottom right
4 2 bottom right
6 4 bottom right