6 nRootGrids = number of root level grids (if >1 repeat all previous input data) F Restart = restart flag 1 MaxLevel = allowed number of grid refinement levels 2 2 2 2 2 2 2 grid refinement ratios 2.0e-3 2.0e-3 2.0e-3 2.0e-3 2.0e-3 2.0e-3 2.0e-3 field variable tolerances on each level 1.0e-2 1.0e-2 1.0e-2 1.0e-2 1.0e-2 1.0e-2 1.0e-2 spatial variable tolerances on each level 0 0 0 0 0 0 0 xbear run log options 2 iBuffer = size of buffer zone around area flagged for refinement 4 MinimumGridPoints along one dimension 0.85 DesiredFillRatio in determining a new subgrid 5 MaxCFLRetry = maximum number of times to reduce dt so that cfl<=cflv(1) before giving up 4 Output style = 1 AMRCLAW, 2 TECPLOT, 3 HDF, 4 UserHDF, 5 UserOut 40 nout = number of output times to print results 0 InterpOpt = interpolation method 0=minmod, 1=constant, 2=centered, 3=left gradient, 4=right gradient 1 TimeStepMethod 0=fixed dt, 1=variable dt 0.0d0 t0 = initial time 8.0d0 tfinal = final time 0 Reserved option (unused at present) 0 Reserved option (unused at present) 0 Reserved option (unused at present) 0 Reserved option (unused at present) 0 Reserved option (unused at present) 2 nDim !!! Data for gas in left part of shock tube !!! 2 MethodOrder = order of precision. used in error estimator 80 mx = cells in x direction 80 my = cells in y direction 1 80 mGlobal(1) = Grid extents along dimension 1 in global indexing system 1 80 mGlobal(2) = Grid extents along dimension 2 in global indexing system 0.0d0 xlower = left edge of computational domain 2.0d0 xupper = right edge of computational domain -1.0d0 ylower = bottom edge of computational domain 1.0d0 yupper = top edge of computational domain 2 mbc = number of ghost cells at each boundary 3 mthbc(1) = type of boundary conditions at left 999 mthbc(2) = type of boundary conditions at right 13 mthbc(3) = type of boundary conditions at bottom 15 mthbc(4) = type of boundary conditions at top 0.5d-3 dtv(1) = initial dt (used in all steps if method(1)=0) 1.0d99 dtv(2) = max allowable dt 0.90d0 cflv(1) = max allowable Courant number 0.80d0 cflv(2) = desired Courant number 4 1 NrVars,nEquationSet = number of primary fields, equation set for these fields 7 maux = number of auxilliary fields 0 method(1) = (unused) 2 method(2) = order 2 method(3) = transverse order 0 method(4) = (unused) 0 method(5) = source term splitting 0 method(6) = (unused) 7 method(7) = maux (should agree with parameter in driver) 4 mwaves = number of waves in each Riemann solution 4 4 4 4 mthlim(mw) = limiter for each wave (mw=1,mwaves) 2 nDim !!! Data for gas in right part of shock tube !!! 2 MethodOrder = order of precision. used in error estimator 80 mx = cells in x direction 80 my = cells in y direction 81 160 mGlobal(1) = Grid extents along dimension 1 in global indexing system 1 80 mGlobal(2) = Grid extents along dimension 2 in global indexing system 2.0d0 xlower = left edge of computational domain 4.0d0 xupper = right edge of computational domain -1.0d0 ylower = bottom edge of computational domain 1.0d0 yupper = top edge of computational domain 2 mbc = number of ghost cells at each boundary 999 mthbc(1) = type of boundary conditions at left 3 mthbc(2) = type of boundary conditions at right 14 mthbc(3) = type of boundary conditions at bottom 16 mthbc(4) = type of boundary conditions at top 0.5d-3 dtv(1) = initial dt (used in all steps if method(1)=0) 1.0d99 dtv(2) = max allowable dt 0.90d0 cflv(1) = max allowable Courant number 0.80d0 cflv(2) = desired Courant number 4 2 NrVars,nEquationSet = number of primary fields, equation set for these fields 7 maux = number of auxilliary fields 0 method(1) = (unused) 2 method(2) = order 2 method(3) = transverse order 0 method(4) = (unused) 0 method(5) = source term splitting 0 method(6) = (unused) 7 method(7) = maux (should agree with parameter in driver) 4 mwaves = number of waves in each Riemann solution 4 4 4 4 mthlim(mw) = limiter for each wave (mw=1,mwaves) 1 nDim !!! Data for string AB !!! 2 MethodOrder = order of precision. used in error estimator 80 mx = cells in x direction 1001 1080 mGlobal(1) = Grid extents along dimension 1 in global indexing system 0.0d0 xlower = left edge of computational domain 2.0d0 xupper = right edge of computational domain 2 mbc = number of ghost cells at each boundary 11 mthbc(1) = type of boundary conditions at left 12 mthbc(2) = type of boundary conditions at right 0.5d-3 dtv(1) = initial dt (used in all steps if method(1)=0) 1.0d99 dtv(2) = max allowable dt 0.90d0 cflv(1) = max allowable Courant number 0.80d0 cflv(2) = desired Courant number 3 3 NrVars,nEquationSet = number of primary fields, equation set for these fields 0 maux = number of auxilliary fields 0 method(1) = (unused) 2 method(2) = order 2 method(3) = transverse order 0 method(4) = (unused) 1 method(5) = source term splitting 0 method(6) = (unused) 0 method(7) = maux (should agree with parameter in driver) 2 mwaves = number of waves in each Riemann solution 4 4 mthlim(mw) = limiter for each wave (mw=1,mwaves) 1 nDim !!! Data for string BC !!! 2 MethodOrder = order of precision. used in error estimator 80 mx = cells in x direction 1081 1160 mGlobal(1) = Grid extents along dimension 1 in global indexing system 2.d0 xlower = left edge of computational domain 4.0d0 xupper = right edge of computational domain 2 mbc = number of ghost cells at each boundary 11 mthbc(1) = type of boundary conditions at left 12 mthbc(2) = type of boundary conditions at right 0.5d-3 dtv(1) = initial dt (used in all steps if method(1)=0) 1.0d99 dtv(2) = max allowable dt 0.90d0 cflv(1) = max allowable Courant number 0.80d0 cflv(2) = desired Courant number 3 4 NrVars,nEquationSet = number of primary fields, equation set for these fields 0 maux = number of auxilliary fields 0 method(1) = (unused) 2 method(2) = order 2 method(3) = transverse order 0 method(4) = (unused) 1 method(5) = source term splitting 0 method(6) = (unused) 0 method(7) = maux (should agree with parameter in driver) 2 mwaves = number of waves in each Riemann solution 4 4 mthlim(mw) = limiter for each wave (mw=1,mwaves) 1 nDim !!! Data for string DE !!! 2 MethodOrder = order of precision. used in error estimator 80 mx = cells in x direction 1161 1240 mGlobal(1) = Grid extents along dimension 1 in global indexing system 0.0d0 xlower = left edge of computational domain 2.0d0 xupper = right edge of computational domain 2 mbc = number of ghost cells at each boundary 11 mthbc(1) = type of boundary conditions at left 12 mthbc(2) = type of boundary conditions at right 0.5d-3 dtv(1) = initial dt (used in all steps if method(1)=0) 1.0d99 dtv(2) = max allowable dt 0.90d0 cflv(1) = max allowable Courant number 0.80d0 cflv(2) = desired Courant number 3 5 NrVars,nEquationSet = number of primary fields, equation set for these fields 0 maux = number of auxilliary fields 0 method(1) = (unused) 2 method(2) = order 2 method(3) = transverse order 0 method(4) = (unused) 1 method(5) = source term splitting 0 method(6) = (unused) 0 method(7) = maux (should agree with parameter in driver) 2 mwaves = number of waves in each Riemann solution 4 4 mthlim(mw) = limiter for each wave (mw=1,mwaves) 1 nDim !!! Data for string BC !!! 2 MethodOrder = order of precision. used in error estimator 80 mx = cells in x direction 1241 1320 mGlobal(1) = Grid extents along dimension 1 in global indexing system 2.0d0 xlower = left edge of computational domain 4.0d0 xupper = right edge of computational domain 2 mbc = number of ghost cells at each boundary 11 mthbc(1) = type of boundary conditions at left 12 mthbc(2) = type of boundary conditions at right 0.5d-3 dtv(1) = initial dt (used in all steps if method(1)=0) 1.0d99 dtv(2) = max allowable dt 0.90d0 cflv(1) = max allowable Courant number 0.80d0 cflv(2) = desired Courant number 3 6 NrVars,nEquationSet = number of primary fields, equation set for these fields 0 maux = number of auxilliary fields 0 method(1) = (unused) 2 method(2) = order 2 method(3) = transverse order 0 method(4) = (unused) 1 method(5) = source term splitting 0 method(6) = (unused) 0 method(7) = maux (should agree with parameter in driver) 2 mwaves = number of waves in each Riemann solution 4 4 mthlim(mw) = limiter for each wave (mw=1,mwaves)