1                   method: 1-upwind, 2-leapfrog, 3-lax-wendroff
1500                   ntmax (max no of time steps)
 12.0                  tmax: max time for integration
 -3.0                  xmin
 15.0                  xmax
   .001                alpha: diffusion coefficient 
  1.0                  u (velocity)
  0.5                  c = u*dt/dx
  0.05                 dt
  0                    cordt (dt from: 0-c or 1-dt directly
  0.0                  delta (param for mass operator)
  0.0                  que (param for upwind operator)
 40                    nout: output every nout steps
  1.0                  f0: normalization of function (t=f0*fnew)
   1                   init: 1-trunc sine, 2-sine, 3-step function
  1.5                  lambda: half wavelength for initial condition
   1                   parbound: 1-periodic, 2-dirichlet, 3-neumann
  0.0                  fb(1): boundary value or derivative at xmin 
  0.0                  fb(2): boundary value or derivative at xmax 
 401                   maxex: no of terms in exact solution<=ntex