| 1 | !============================================================================================
|
|---|
| 2 | ! AstroBEAR global.data File
|
|---|
| 3 | ! Global parameters, Contains GlobalData and LevelData and PoissonData namelists.
|
|---|
| 4 | ! Configured for: Template
|
|---|
| 5 | ! Default values for variables in [brackets]
|
|---|
| 6 | !============================================================================================
|
|---|
| 7 |
|
|---|
| 8 | &GlobalData
|
|---|
| 9 |
|
|---|
| 10 | !=============================================================================================
|
|---|
| 11 | ! Parameters Related to Space
|
|---|
| 12 | !=============================================================================================
|
|---|
| 13 | nDim = 2 ! number of dimensions for this problem (1-3)
|
|---|
| 14 | GmX = 640,640,1 ! Base grid resolution [x,y,z]
|
|---|
| 15 | MaxLevel = 0 ! Maximum level for this simulation (0 is fixed grid)
|
|---|
| 16 | LastStaticLevel = -1 ! Use static AMR for levels through LastStaticLevel [-1]
|
|---|
| 17 | GxBounds = 0.0,0.0,0.0,10.0,10.0,0.0 ! Problem boundaries in computational units,format:
|
|---|
| 18 | ! (xlower, ylower, zlower, xupper, yupper, zupper)
|
|---|
| 19 | ! For 2D problems, set zlower and zupper to 0.d0.
|
|---|
| 20 |
|
|---|
| 21 | Gmthbc = 1,1,1,1,1,1 ! Sets the physical boundary conditions at the edge of the problem domain
|
|---|
| 22 | ! format: (x1, y1, z1, x2, y2, z2)
|
|---|
| 23 | ! 1-Extrapolated,
|
|---|
| 24 | ! 2-Periodic,
|
|---|
| 25 | ! 3-ReflectingWall (Field lines do not penetrate)
|
|---|
| 26 | ! 4-Reflect_BParallel (Field lines held normal
|
|---|
| 27 | ! 5-Reflect_Cylindrical (Like reflect wall, but also changes sign of phi
|
|---|
| 28 | ! components of velocity and magnetic fields)
|
|---|
| 29 | ! [1,1,1,1,1,1]
|
|---|
| 30 |
|
|---|
| 31 | !=============================================================================================
|
|---|
| 32 | ! Parameters Related to Time and output
|
|---|
| 33 | !=============================================================================================
|
|---|
| 34 | start_time = 0.0d0 ! The initial time in computational units [0d0]. Do NOT change this on restarts.
|
|---|
| 35 | final_time = 0.00123958795 ! The final time in computational units.
|
|---|
| 36 | final_frame = 155 ! The final frame [10]
|
|---|
| 37 | lRestart = F ! restart flag [F]
|
|---|
| 38 | lPostProcess = F ! Post Process (preempts lrestart) [F]
|
|---|
| 39 | restart_frame = -1 !The frame from which restarts should begin (-1 sets restart from last chombo dump). Read only if lRestart or lPostProcess = .true.
|
|---|
| 40 | initial_maxspeed = 10000 ! 1000.d0 ! A lower bound on the initial maximum wave speed in computational units.
|
|---|
| 41 | cfl_vars = 1, .3, .5 ! Maximum allowed CFL before restarting step, desired CFL to run at, relaxation parameter. Smaller numbers give more gradual relaxation [1, .3, .5]
|
|---|
| 42 | iDataFileType = 0 ! Type of output: 0 is chombo, 1 is BOV, & 2 is curve files (1D fixed grid only) [0]
|
|---|
| 43 | lParallelHDF5 = F ! Use parallel IO for chombo dumps (requires phdf5 library and compile flags)
|
|---|
| 44 | lUniformTimeStepping = F ! Advance all levels by same time step. (Turns off child levels taking multiple time steps for each parent step).
|
|---|
| 45 |
|
|---|
| 46 | !=============================================================================================
|
|---|
| 47 | ! Parameters Related to Performance/Algorithms
|
|---|
| 48 | !=============================================================================================
|
|---|
| 49 | lSkipProfile = T ! Flag for skipping profiling [T] - should set this to false for good performance on production runs.
|
|---|
| 50 | lUseOriginalNewSubGrids = F ! T for Old Subgrid Algo, F for New Subgrid Algo. [F]
|
|---|
| 51 | MinimumGridPoints = 2 ! The minimum number of cells a grid can have along one dimension [2] (ignored unless lUseOriginalNewSubGrids=T).
|
|---|
| 52 | iThreaded = -1 ! -1 is non-threaded, 0 is pseudo-threaded, 1 is threaded [0]. If you are using threaded you must compile with pThreadFlag=2
|
|---|
| 53 | LevelBalance = 0d0,0d0 ! Try to balance load on every level? [0d0,0d0]
|
|---|
| 54 |
|
|---|
| 55 |
|
|---|
| 56 | !=============================================================================================
|
|---|
| 57 | ! Parameters Related to Testing/Debugging
|
|---|
| 58 | !=============================================================================================
|
|---|
| 59 |
|
|---|
| 60 | lTest = F ! Use testing objects to compare results with reference data [F]
|
|---|
| 61 | lTimingLog = F ! Output log showing timing information [F]
|
|---|
| 62 | lStressTest = F ! Used to manually configure AMR patches and produce random fluxes and divergence [F]
|
|---|
| 63 | lSkipHydro = F ! Used if you want to turn of hydrodynamics and just see results of multi-physics (diffusion, radtransfer, etc...) [F]
|
|---|
| 64 | /
|
|---|
| 65 |
|
|---|
| 66 | !============================================================================================
|
|---|
| 67 | !Level-specific variables (All LevelData lists should be initialized from 0 to MaxLevel).
|
|---|
| 68 | !============================================================================================
|
|---|
| 69 |
|
|---|
| 70 | &LevelData
|
|---|
| 71 | qTolerance = 1.d-1,1.d-1,1.d-1,1.d-1,1.d-1,1.d-1,1.d-1,1.d-1 ! field variable tolerances on each level
|
|---|
| 72 | DesiredFillRatios = .70,.70,.70,.70,.70,.70,.70,.70 ! Determines how much of the new subgrids must be filled with refined cells (also ignored unless lUseOriginalNewSubGrids=T)
|
|---|
| 73 | buffers = 2 2 2 2 2 2 2 2 ! Number of zones to buffer ghost cells
|
|---|
| 74 | /
|
|---|
| 75 |
|
|---|
| 76 | !============================================================================================
|
|---|
| 77 | !Variables used for self gravity module
|
|---|
| 78 | !============================================================================================
|
|---|
| 79 |
|
|---|
| 80 | &PoissonData
|
|---|
| 81 | solver = 1 ! Which solver to use for linear system solve? StructPCG=1, StructGMRes=2 [1]
|
|---|
| 82 | tolerance = 1e-6 ! Solver convergence tolerance [1e-6]
|
|---|
| 83 | printlevel = 0 ! Sets print level for hypre solver. Increasing this will make solvers more verbose. [0]
|
|---|
| 84 | MaxIters = 1000 ! Maximum number of iterations to attempt convergence. [1000]
|
|---|
| 85 | hverbosity = 0 ! Anything > 0 will cause linear systems to be outputted as matrix, rhs, lhs files. [0]
|
|---|
| 86 | mthbc = 2 2 2 2 2 2 ! What type of boundary conditions to use for gravitational potential
|
|---|
| 87 | ! format: (x1, y1, z1, x2, y2, z2)
|
|---|
| 88 | ! 0-Zero Slope,
|
|---|
| 89 | ! 1-Reflecting,
|
|---|
| 90 | ! 2-Periodic,
|
|---|
| 91 | ! 3-Multipole Expansion,
|
|---|
| 92 | ! 4-User Specified (User is reponsible for setting Phi at physical boundaries)
|
|---|
| 93 | /
|
|---|
| 94 |
|
|---|
| 95 | !============================================================================================
|
|---|
| 96 | !Variables used for radiative transfer module
|
|---|
| 97 | !============================================================================================
|
|---|
| 98 |
|
|---|
| 99 | &RadTransferData
|
|---|
| 100 | solver = 1 ! Which solver to use for linear system solve? StructPCG=1, StructGMRes=2, StructSMG=3 [1]
|
|---|
| 101 | tolerance = 1e-6 ! Solver convergence tolerance [1e-6]
|
|---|
| 102 | printlevel = 0 ! Sets print level for hypre solver. Increasing this will make solvers more verbose. [0]
|
|---|
| 103 | MaxIters = 1000 ! Maximum number of iterations to attempt convergence. [1000]
|
|---|
| 104 | hverbosity = 0 ! Anything > 0 will cause linear systems to be outputted as matrix, rhs, lhs files. [0]
|
|---|
| 105 | mthbc = 2 2 2 2 2 2 ! What type of boundary conditions to use for gravitational potential
|
|---|
| 106 | ! format: (x1, y1, z1, x2, y2, z2)
|
|---|
| 107 | ! 0-Free Streaming,
|
|---|
| 108 | ! 1-Extrapolated,
|
|---|
| 109 | ! 2-Periodic
|
|---|
| 110 | ! 3-Reflecting
|
|---|
| 111 | ! 4-User defined boundary
|
|---|
| 112 | ! 5-User defined radiative flux
|
|---|
| 113 | RadFlux = 0d0 ! User defined radiative flux used for mthbc 5
|
|---|
| 114 | Crankiness = 0 ! parameter that adjusts from backward euler (0) to crank-nicholson(1) time stepping. (2 is forward euler but nothing over 1 should be used - and 0 should be more stable but less accurate)
|
|---|
| 115 | cfl_rad = .01 ! Fractional change in temperature allowed before recalculating source function
|
|---|
| 116 | MaxSubCycles = 1000 ! Maximum number of sub-cycles to allow. May result in increase in cfl_rad.
|
|---|
| 117 | /
|
|---|
