1.e0     umax = speed limit in flux function
0.4e0    q0 = initial uniform density
0.008e0  D = strength of source term