 maxstep  = 20,              # number of steps to run OR
 stoptime = 100.0,           # max simulation time
 intplot  = 5,               # who often to write plot file
 maxplots = 1,               # ignore, just set to 1
 nxcell   = 40,              # number of cells in X direction
 nycell   = 20,              # number of cells in Y direction
 xlft     = 0.0,             # physical position of left edge of domain
 xrgt     = 20.0,	     # physical position of right edge of domain
 ybot     = 0.0,	     # physical position of bottom edge of domain
 ytop     = 10.0,	     # physical position of top edge of domain
 cfl      = 0.9,             # cfl number, 0.9 is good
 gamma    = 1.4,             # polytropic gas constant
 difmag   = 0.1,             # coef of artificial visc.
 timefact = 0.5,             # initial scaling of timestep
 lstrip   = .false.,          # stripmine if true
 ldebug   = .true.,          # ??? i forgot
 lnewton  = .false.,         # print out stats of riemann solver
 ncheck   = 100,             # how often to write checkpoint file
 nupdate  = 1,               # how often to inform user of progress
 niter    = 2,               # number of iterations in Riemann solver

 # ---------------------------------------------------------------
 # problem specific inputs
 shockdir   = 1.0,           # direction of shock, 1=x, 2=y
 shockpos   = 5.0,           # position of shock front
 rmach      = 2.0,           # shock mach number
 xcld       = 10.0,          # (xcld,ycld) = position of circular
 ycld       = 5.0,           #               cloud of gas
 radius     = 3.0,           # radius of cloud
 denfact    = 2.0,           # density ratio between cloud and ambient
 xperiod    = .false.,       # true if periodic in x direction
 yperiod    = .false.,       # true if periodic in y direction
 godorder   = 1		     # Godunov order
