40                mx          = cells in x direction
 40                my          = cells in y direction
 3                 mxnest      = max number of grid levels
 2 4 2 2           inrat       = refinement ratios (1..mxnest)


 5                 nout        = number of output times to print results
2                  outstyle    = style of specifying output times
0.1 0.2 0.3 0.4 0.5 tout  = output times

0.0005d0           dtv(1)      = initial dt (used in all steps if method(1)=0)
1.0d99             dtv(2)      = max allowable dt 
1.0d0              cflv(1)     = max allowable Courant number
0.9d0              cflv(2)     = desired Courant number
500                nv(1)       = max number of time steps per call to claw2

1                  method(1)   = 1 for variable dt
2                  method(2)   = order 
2                  method(3)   = transverse order
3                  method(4)   = verbosity of output
1                  method(5)   = source term splitting
0                  method(6)   = mcapa
2                  method(7)   = maux (should agree with parameter in driver)
 "center"       aux(i,j,1) = impedance
 "center"       aux(i,j,2) = sound speed

3  		meqn     (# variables - e.g. 4 for Euler)
2 		mwaves for clawpack  (number of waves in soln to Riemann prob)
4 4 4             (mthlim(mw), mw=1,mwaves)    selects limiter

0.d0               t0          = initial time
0.0d0              xlower      = left edge of computational domain
1.0d0              xupper      = right edge of computational domain
0.0d0              ylower      = bottom edge of computational domain
1.0d0              yupper      = top edge of computational domain

2                  mbc         = number of ghost cells at each boundary
1                  mthbc(1)    = type of boundary conditions at left
1                  mthbc(2)    = type of boundary conditions at right
1                  mthbc(3)    = type of boundary conditions at bottom
1                  mthbc(4)    = type of boundary conditions at top

F               do a restart (input from rest2d)
1000            iousr

0.05 	        tol 	 (tolerance for Richardson estimator)
0.05            tolsp	 (spatial tolerance for refinement)
2               kcheck   (how often to est. error - related to buffer size)
3               ibuff    (buffer zone size - add # of flagged pts to add) 
.70             cutoff   (efficiency cutoff for grid generator)


F               PRINT option  (if T prints soln. to out2d)
F		NCAR graphics
T               Matlab graphics


F		dprint -  verbose domain flags
F		eprint -  error estimation - output the flags
F		edebug -  even more error est. output
F		gprint - grid generation output (bisection, clustering,...)
F		nprint - proper nesting output
F		pprint - projectiong of tagged pts. output
F		rprint - verbose regridding - output new grid summary
F		sprint - space (memory) output
T		tprint - time step reporting - for each level
F		uprint - updating/upbnding reporting