1.0e0      advection velocity u
 100.e0      parameter determining width of wave packet
 80.e0       frequency of wave packet
 6500.0e0       phiM