Guest User

flash.par

a guest
Jun 24th, 2014
281
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 5.08 KB | None | 0 0
  1. # Parameters file for core-collapse SN
  2.  
  3. # Simulation basics
  4. basenm = "a140_"
  5. restart = .false.
  6. checkpointFileNumber = 0
  7. plotFileNumber = 0
  8. run_comment = "a140 1D with radiation"
  9. #output_directory = "output"
  10.  
  11.  
  12. # IO
  13. checkpointFileIntervalStep = 0
  14. checkpointFileIntervalTime = 0.05
  15. plotFileIntervalStep = 0
  16. plotFileIntervalTime = 0.
  17. # wall_clock_checkpoint = 43200.
  18. wall_clock_time_limit = 21000.
  19.  
  20. # Time
  21. tinitial = 0.0
  22. tmax = 3.9
  23. nend = 100000000
  24. tstep_change_factor = 1.5
  25. dtinit = 1.E-8
  26. dtmax = 1.E5
  27. dtmin = 1.E-20
  28.  
  29. # Domain
  30. geometry = "spherical"
  31. xmax = 1.9e11
  32. xmin = 0.0
  33. xl_boundary_type = "reflect"
  34. xr_boundary_type = "outflow"
  35.  
  36. # Grid/Refinement
  37. nblockx = 100
  38. nblocky = 1
  39. nblockz = 1
  40.  
  41. gr_lrefineMaxRedDoByLogR = .true.
  42. gr_lrefineMaxRedDoByTime = .FALSE.
  43. gr_lrefineMaxRedLogBase = 10.0
  44. gr_lrefineMaxRedRadiusFact = 0.15
  45. gr_lrefineMaxRedTRef = 2.0
  46. gr_lrefineMaxRedTimeScale = 0.5
  47.  
  48. lrefine_max = 3
  49. lrefine_min = 1
  50. refine_var_1 = "dens"
  51. refine_var_2 = "pres"
  52. refine_var_3 = "none"
  53. refine_var_4 = "none"
  54. refine_cutoff_1 = 0.8
  55. refine_cutoff_2 = 0.8
  56. refine_cutoff_3 = 0.8
  57. refine_cutoff_4 = 0.8
  58.  
  59. # Simulation
  60. model_file = "a140.flash"
  61. nsub = 4
  62. vel_mult = 1.0
  63.  
  64. # Hydro
  65. useHydro = .TRUE.
  66. cfl = 0.5
  67. interpol_order = 1
  68. updateHydroFluxes = .TRUE.
  69. eintSwitch = 0.0 # Always use Etot
  70. convertToConsvdForMeshCalls = .false.
  71. converttoconsvdinmeshinterp = .false.
  72.  
  73. ## SWITCHES SPECIFIC TO THE UNSPLIT HYDRO SOLVER ##
  74. # I. INTERPOLATION SCHEME:
  75. order = 3 # Interpolation order (first/second/third/fifth order)
  76. slopeLimiter = "hybrid" # Slope limiters (minmod, mc, vanLeer, hybrid, limited)
  77. LimitedSlopeBeta= 1. # Slope parameter for the "limited" slope by Toro
  78. charLimiting = .true. # Characteristic limiting vs.Primitive limiting
  79.  
  80. use_avisc = .true. # use artificial viscosity (originally for PPM)
  81. cvisc = 0.1 # coefficient for artificial viscosity
  82. use_flattening = .true. # use flattening (dissipative) (originally for PPM)
  83. use_steepening = .false. # use contact steepening (originally for PPM)
  84. use_upwindTVD = .false. # use upwind biased TVD slope for PPM (need nguard=6)
  85. flux_correct = .true.
  86. EOSforRiemann = .false.
  87. transOrder = 1
  88. use_auxEintEqn = .TRUE.
  89. hydroComputeDtOption = -1
  90.  
  91. # II. RIEMANN SOLVERS:
  92. RiemannSolver = "hybrid" # Roe, HLL, HLLC, LLF, Marquina
  93. entropy = .false. # Entropy fix for the Roe solver
  94.  
  95. # III. STRONG SHOCK HANDELING SCHEME:
  96. shockDetect = .true. # Shock Detect for numerical stability
  97. ## -------------------------------------------------------------##
  98.  
  99. # Gravity
  100. useGravity = .true.
  101. updateGravity = .TRUE.
  102. grav_boundary_type = "isolated"
  103. mpole_3daxisymmetry = .false.
  104. mpole_dumpMoments = .FALSE.
  105. mpole_PrintRadialInfo = .false.
  106. mpole_IgnoreInnerZone = .false.
  107. mpole_lmax = 16
  108. mpole_ZoneRadiusFraction_1 = 1.0
  109. mpole_ZoneExponent_1 = 0.005
  110. mpole_ZoneScalar_1 = 0.5
  111. mpole_ZoneType_1 = "logarithmic"
  112. point_mass = 0.0
  113. point_mass_rsoft = 0.e0
  114. use_gravHalfUpdate = .TRUE.
  115. use_gravConsv = .FALSE.
  116. use_gravPotUpdate = .FALSE.
  117.  
  118.  
  119. # EOS
  120. eosMode = "dens_ie"
  121. eosModeInit = "dens_temp_equi"
  122. gamma = 1.66666666666667
  123.  
  124. # Small numbers
  125. smallt = 1.0e-20
  126. smlrho = 1.0e-27
  127. smltemp = 1.0e-2
  128. smallp = 1.E-20
  129. smalle = 1.E-20
  130. smallu = 1.E-10
  131. smallx = 1.E-100
  132. small = 1.E-100
  133.  
  134.  
  135.  
  136. ### Multigroup diffusion options ###
  137. rt_useMGD = .true.
  138. rt_mgdNumGroups = 1
  139. rt_mgdBounds_1 = 0.0
  140. rt_mgdBounds_2 = 1.0e14
  141.  
  142.  
  143. op_emitConst = 10.0
  144. op_transConst = 1.0e-06
  145. op_absorbConst = 10.0
  146.  
  147.  
  148. # Set boundary conditions:
  149. rt_mgdXlBoundaryType = "vacuum"
  150. rt_mgdXrBoundaryType = "vacuum"
  151. rt_mgdYlBoundaryType = "vacuum"
  152. rt_mgdYrBoundaryType = "vacuum"
  153. rt_mgdZlBoundaryType = "vacuum"
  154. rt_mgdZrBoundaryType = "vacuum"
Advertisement
Add Comment
Please, Sign In to add comment