4 3 GfsSimulation GfsBox GfsGEdge {} {
  Time { end = 1 dtmax = 5e-3 }
  Refine 7

  # The tracer T is used to track both phases
  VariableTracerVOF {} T
 
  # The initial sinusoidal interface (translated by 0.5 along the y-axis)
  InitFraction {} T (0.05*cos (2.*M_PI*x) + y) { ty = 0.5 }

  AdaptVorticity { istep = 1 } { maxlevel = 7 cmax = 2e-2 }
  AdaptGradient { istep = 1 } { maxlevel = 7 cmax = 1e-2 } T

  # The dynamic viscosity for both phases
  SourceViscosity {} 0.00313

  # This defines the inverse of the density of the fluids as a
  # function of T
  PhysicalParams { alpha = 1./(T*1.225 + (1. - T)*0.1694) }

  # We also need gravity
  Source {} V -9.81

  OutputTime { istep = 10 } stderr
  OutputBalance { istep = 10 } stderr
  OutputProjectionStats { istep = 10 } stderr
  OutputDiffusionStats { istep = 10 } stderr
  OutputPPM { istep = 2 } { ppm2mpeg > vort.mpg} {
    min = -30 max = 30 v = Vorticity
  }
  OutputPPM { istep = 2 } { ppm2mpeg > t.mpg } {
    min = 0 max = 1 v = T
  }
  OutputPPM { start = end } { convert -colors 256 ppm:- vort.eps } {
    min = -30 max = 30 v = Vorticity
  }
  OutputPPM { start = end } { convert -colors 256 ppm:- t.eps } {
    min = 0 max = 1 v = T
  }
  OutputTiming { start = end } stderr
  OutputSimulation { step = 0.1 } stdout
  EventScript { start = 0 } { echo "Save t-0.eps { format = EPS }" }
  EventScript { start = 0.7 step = 0.1 } { echo "Save t-$GfsTime.eps { format = EPS }" }
}
GfsBox {}
GfsBox {}
GfsBox {}
GfsBox {}
1 2 top
2 3 top
1 4 bottom