Mathematica code:

S[n_, t_] := Sin[n*2 Pi/50 + t]

ListAnimate[
 Show[
   Table[
    Plot[
     100 - n + (14Sin[t]*S[n,t] + .02)*Exp[-(x - 4Sin[t]*S[n, 0])^2/Abs[S[n,t]]], 
     {x, -10, 10},
    PlotStyle -> Directive[Black, Thick], PlotRange -> {{-7, 7}, {0, 100.5}},
    Filling -> Axis, FillingStyle -> White, Axes -> False, AspectRatio -> Full, 
    ImageSize -> {500, 750}],
  {n, -10, 100, 1}]],
{t, .001, 2 Pi + .001, (2 Pi + .001)/40}, AnimationRunning->False]