!:NumericalSchemeParameters:! one value if LevelMethods==F else (>=MaxLevel) values
!==================================================================================
0                  method(1)   = (reserved)
2                  method(2)   = convergence order
2                  method(3)   = transverse convergence order
0                  method(4)   = verbosity of wavebear output
0                  method(5)   = source term splitting
0                  method(6)   = 0 split q differences, 1 split flux differences
0                  method(7)   = radius of slab around current 1D array of cells

4                  mwaves      = number of waves in each Riemann solution
3 3 3 3            mthlim(mw)  = limiter for each wave  (mw=1,mwaves)
!==================================================================================