1.0e0      advection velocity u
 200.e0      parameter determining width of initial Gaussian pulse
 1.0e0      rate lambda in  q_t + uq_x = -lambda q