problem = Keplerian disk in spherical polar coordinates reference = configure = python configure.py --prob=disk_athena --coord=spherical_polar -hdf5 --hdf5_path=/c13scratch/apps/hdf5-parallel/1.8.20 -mpi python configure.py --prob=disk_mhd_new_Dec --coord=spherical_polar -b -hdf5 --hdf5_path=/usr/lib/x86_64-linux-gnu/hdf5/openmpi -mpi -h5double problem_id = disk_mhd_Jan24 # problem ID: basename of output filenames file_type = hdf5 # HDF5 data dump variable = cons,b,bcc,uov # variables to be output dt = 0.1 # time increment between outputs cfl_number = 0.3 # The Courant, Friedrichs, & Lewy (CFL) Number nlim = -1 # cycle limit tlim = 10000.0 # time limit integrator = vl2 # time integration algorithm xorder = 2 # order of spatial reconstruction ncycle_out = 1 # interval for stdout summary info nx1 = 1024 # Number of zones in X1-direction x1min = 0.1 # minimum value of X1 x1max = 70.0 # maximum value of X1 ix1_bc = outflow # Inner-X1 boundary condition flag ox1_bc = outflow # Outer-X1 boundary condition flag x1rat = 1.0219 # non-uniform grid ratio nx2 = 512 # Number of zones in X2-direction x2min = 0.1 # minimum value of X2 x2max = 3.0 # maximum value of X2 ix2_bc = outflow # Inner-X2 boundary condition flag ox2_bc = outflow # Outer-X2 boundary condition flag nx3 = 32 # Number of zones in X3-direction x3min = 0.0 # minimum value of X3 x3max = 6.2831853071795862 # maximum value of X3 ix3_bc = periodic # Inner-X3 boundary condition flag ox3_bc = periodic # Outer-X3 boundary condition flag refinement = none nx1 = 32 nx2 = 32 nx3 = 32 x1min = 0.1 x1max = 40.0 x2min = 0.8 x2max = 1.8 x3min = 0.0 x3max = 6.283185307179586 level=1 x1min = 0.1 x1max = 30.0 x2min = 1.2 x2max = 1.77 x3min = 0.0 x3max = 6.283185307179586 level=2 x1min = 0.1 x1max = 20.0 x2min = 1.4 x2max = 1.6 x3min = 0.0 x3max = 6.283185307179586 level=3 gamma = 1.66666667 # gamma = C_p/C_v iso_sound_speed = 0.1 # isothermal sound speed dfloor = 1.e-7 pfloor = 1.e-10 OAorder = 0 # 0: w/o OA, 1: w/ OA (1st), 2: w/ OA (2nd) Omega0 = 0.0 # angular velocity of the system GM = 1.0 r0 = 1.0 rho0 = 1.0 p0_over_r0 = 0.03 pslope = -1.0 dslope = -1.5 rmin = 0.1 beta0 = 5000