API tools faq
paste
Advertisement
Guest User

Untitled

a guest
Jun 12th, 2019
237
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 27.87 KB | None | 0 0
raw download clone embed print report
  1. ######################################################################
  2. # Simulation of supermassive toroidal stars #
  3. # Authors: Gabriele Bozzola, Vasileios Paschalidis #
  4. # <bozzola.gabriele@gmail.com> #
  5. # Steward Observatory, University of Arizona #
  6. ######################################################################
  7.  
  8. ######################################################################
  9. # #
  10. # VERBOSITY #
  11. # #
  12. ######################################################################
  13.  
  14. # Carpet::storage_verbose = "no"
  15. # Carpet::verbose = "no"
  16. # Carpet::veryverbose = "no"
  17.  
  18. # CarpetLib::output_bboxes = "no"
  19.  
  20. # CarpetMask::verbose = "no"
  21. # CarpetReduce::verbose = "no"
  22. # CarpetRegrid2::verbose = "no"
  23. # CarpetRegrid2::veryverbose = "no"
  24.  
  25. # NaNChecker::verbose = "all"
  26.  
  27. # # QuasiLocalMeasures::verbose = "no"
  28. # # QuasiLocalMeasures::veryverbose = "no"
  29. # SphericalSurface::verbose = "no"
  30.  
  31. NSTracker::NSTracker_verbose = "no"
  32. RhoBFinder::verbosity_level = 0
  33.  
  34. # # LapseFinder::verbosity_level = 2
  35.  
  36. # # WeylScal4::verbose = 0
  37.  
  38. # # AHFinderDirect::verbose_level = "physics details"
  39.  
  40. # # Dissipation::verbose = "no"
  41.  
  42. Trigger::Trigger_Debug = 0
  43.  
  44. Carpet::grid_structure_filename = "carpet-grid-structure"
  45. Carpet::grid_coordinates_filename = "carpet-grid-coordinates"
  46.  
  47.  
  48.  
  49. ######################################################################
  50. # #
  51. # TIME #
  52. # #
  53. ######################################################################
  54.  
  55. ActiveThorns = "MoL Time"
  56.  
  57. # Cactus::cctk_itlast = 10
  58. Cactus::terminate = "time"
  59. Cactus::cctk_final_time = 10000
  60.  
  61. Time::dtfac = 0.025
  62. # This is set so that the courant factor is never greater than 0.5,
  63. # which leads to numerical instabilities
  64. Carpet::time_refinement_factors = "[1,1,1,1,1,2,4,8,16,32,64]"
  65. MoL::ODE_Method = "rk4"
  66. MoL::MoL_Intermediate_Steps = 4
  67. MoL::MoL_Num_Scratch_Levels = 1
  68.  
  69.  
  70. ######################################################################
  71. # #
  72. # GRID #
  73. # #
  74. ######################################################################
  75.  
  76. ActiveThorns = "Coordbase CartGrid3d Boundary StaticConformal"
  77. ActiveThorns = "Carpet CarpetLib CarpetReduce CarpetRegrid2 CarpetInterp
  78. CarpetInterp2 CarpetMask"
  79.  
  80. Grid::avoid_origin = "yes"
  81. CartGrid3D::type = "coordbase"
  82. CoordBase::domainsize = "minmax"
  83. CoordBase::xmin = -4095.99
  84. CoordBase::ymin = -4096.00
  85. CoordBase::zmin = -4096.00
  86. CoordBase::xmax = 4096.01
  87. CoordBase::ymax = 4096.00
  88. CoordBase::zmax = 4096.00
  89.  
  90. # This is 3M/64 in the finest refinement level
  91. CoordBase::dx = 32.0
  92. CoordBase::dy = 32.0
  93. CoordBase::dz = 32.0
  94. # CoordBase::dx = 48.0
  95. # CoordBase::dy = 48.0
  96. # CoordBase::dz = 48.0
  97.  
  98. # This is 3M/96 in the finest refinement level
  99. # CoordBase::dx = 36.0
  100. # CoordBase::dy = 36.0
  101. # CoordBase::dz = 36.0
  102.  
  103. # This is 3M/128 in the finest refinement level
  104. # CoordBase::dx = 24.0
  105. # CoordBase::dy = 24.0
  106. # CoordBase::dz = 24.0
  107.  
  108. CoordBase::boundary_size_x_lower = 3
  109. CoordBase::boundary_size_y_lower = 3
  110. CoordBase::boundary_size_z_lower = 3
  111. CoordBase::boundary_size_x_upper = 3
  112. CoordBase::boundary_size_y_upper = 3
  113. CoordBase::boundary_size_z_upper = 3
  114.  
  115. Carpet::ghost_size = 3
  116. Carpet::domain_from_coordbase = "yes"
  117. Carpet::use_buffer_zones = "yes"
  118.  
  119. # Refinement Levels
  120. Carpet::max_refinement_levels = 11
  121.  
  122. # Fragments
  123.  
  124. CarpetRegrid2::num_levels_1 = 8
  125. # CarpetRegrid2::position_x_1 = 15
  126. # CarpetRegrid2::position_y_1 = 15
  127. CarpetRegrid2::radius_1[1] = 2048.0
  128. CarpetRegrid2::radius_1[2] = 1024.0
  129. CarpetRegrid2::radius_1[3] = 512.0
  130. CarpetRegrid2::radius_1[4] = 256.0
  131. CarpetRegrid2::radius_1[5] = 128.0
  132. CarpetRegrid2::radius_1[6] = 64.0
  133. CarpetRegrid2::radius_1[7] = 32.0
  134. CarpetRegrid2::radius_1[8] = 16.0
  135. CarpetRegrid2::radius_1[9] = 8.0
  136. CarpetRegrid2::radius_1[10] = 4.0
  137.  
  138. CarpetRegrid2::num_levels_2 = 8
  139. # CarpetRegrid2::position_x_2 = -15
  140. # CarpetRegrid2::position_y_2 = -15
  141. CarpetRegrid2::radius_2[1] = 2048.0
  142. CarpetRegrid2::radius_2[2] = 1024.0
  143. CarpetRegrid2::radius_2[3] = 512.0
  144. CarpetRegrid2::radius_2[4] = 256.0
  145. CarpetRegrid2::radius_2[5] = 128.0
  146. CarpetRegrid2::radius_2[6] = 64.0
  147. CarpetRegrid2::radius_2[7] = 32.0
  148. CarpetRegrid2::radius_2[8] = 16.0
  149. CarpetRegrid2::radius_2[9] = 8.0
  150. CarpetRegrid2::radius_2[10] = 4.0
  151.  
  152. # Star
  153.  
  154. CarpetRegrid2::num_levels_3 = 8
  155. CarpetRegrid2::radius_3[1] = 2048.0
  156. CarpetRegrid2::radius_3[2] = 1024.0
  157. CarpetRegrid2::radius_3[3] = 512.0
  158. CarpetRegrid2::radius_3[4] = 256.0
  159. CarpetRegrid2::radius_3[5] = 128.0
  160. CarpetRegrid2::radius_3[6] = 64.0
  161. CarpetRegrid2::radius_3[7] = 32.0
  162. CarpetRegrid2::radius_3[8] = 16.0
  163. CarpetRegrid2::radius_3[9] = 8.0
  164. CarpetRegrid2::radius_3[10] = 4.0
  165.  
  166. # Regridding
  167.  
  168. ActiveThorns = "CarpetTracker NSTracker SphericalSurface"
  169.  
  170. SphericalSurface::nsurfaces = 8
  171. SphericalSurface::maxntheta = 39
  172. SphericalSurface::maxnphi = 76
  173.  
  174. SphericalSurface::ntheta [0] = 39
  175. SphericalSurface::nphi [0] = 76
  176. SphericalSurface::nghoststheta[0] = 2
  177. SphericalSurface::nghostsphi [0] = 2
  178. SphericalSurface::name [0] = "BH Right"
  179.  
  180. SphericalSurface::ntheta [1] = 39
  181. SphericalSurface::nphi [1] = 76
  182. SphericalSurface::nghoststheta[1] = 2
  183. SphericalSurface::nghostsphi [1] = 2
  184. SphericalSurface::name [1] = "BH Left"
  185.  
  186. SphericalSurface::ntheta [2] = 39
  187. SphericalSurface::nphi [2] = 76
  188. SphericalSurface::nghoststheta[2] = 2
  189. SphericalSurface::nghostsphi [2] = 2
  190.  
  191. SphericalSurface::ntheta [3] = 39
  192. SphericalSurface::nphi [3] = 76
  193. SphericalSurface::nghoststheta[3] = 2
  194. SphericalSurface::nghostsphi [3] = 2
  195.  
  196. SphericalSurface::ntheta [4] = 39
  197. SphericalSurface::nphi [4] = 76
  198. SphericalSurface::nghoststheta[4] = 2
  199. SphericalSurface::nghostsphi [4] = 2
  200.  
  201.  
  202. CarpetRegrid2::regrid_every = 32
  203. CarpetRegrid2::num_centres = 3
  204. CarpetRegrid2::freeze_unaligned_levels = "yes"
  205. CarpetRegrid2::freeze_unaligned_parent_levels = "yes"
  206. CarpetRegrid2::ensure_proper_nesting = "yes"
  207. CarpetRegrid2::snap_to_coarse = "yes"
  208.  
  209. CarpetTracker::surface_name[0] = "BH Right"
  210. CarpetTracker::surface_name[1] = "BH Left"
  211.  
  212. # system specific Carpet paramters
  213. Carpet::prolongation_order_space = 5
  214. Carpet::prolongation_order_time = 2
  215.  
  216. # Fill uninitialized memory with large values
  217. CarpetLib::poison_new_memory = "yes"
  218. CarpetLib::poison_value = 113
  219.  
  220.  
  221. # Experiment with these parameters!
  222. CarpetLib::interleave_communications = "yes"
  223. CarpetLib::combine_sends = "yes"
  224.  
  225.  
  226. ######################################################################
  227. # #
  228. # ANALYSIS AND DIAGNOSTICS #
  229. # #
  230. ######################################################################
  231.  
  232.  
  233. ActiveThorns = "LapseFinder RhoBFinder"
  234.  
  235. LapseFinder::LapseFinder_comp_alp_min_every = 0
  236. # LapseFinder::LapseFinder_alp_min_loc_only_positive_x = "true"
  237. # LapseFinder::LapseFinder_average_multiple_minima_locations = "no"
  238. # LapseFinder::LapseFinder_tolerance = 0.00001
  239.  
  240.  
  241. RhoBFinder::RhoBFinder_rho_b_max_loc_only_positive_x = "true"
  242. RhoBFinder::RhoBFinder_average_multiple_maxima_locations = "no"
  243. RhoBFinder::RhoBFinder_compute_only_first = "false"
  244.  
  245. NSTracker::NSTracker_SF_Name = "BH Right"
  246. NSTracker::NSTracker_SF_Name_Opposite = "BH Left"
  247. NSTracker::NSTracker_max_distance = 3
  248. NSTracker::NSTracker_use_AHFinderDirect = "false"
  249. NSTracker::NSTracker_AHFinderDirect_horizon_number = 1
  250.  
  251. NSTracker::NSTracker_tracked_location = "RhoBFinder::RhoBFinder_rho_b_max_loc"
  252. NSTracker::NSTracker_tracked_location_first = "RhoBFinder::RhoBFinder_rho_b_max_loc"
  253. NSTracker::NSTracker_relocate_SF_at_first_iteration = "yes"
  254.  
  255. # Check all physical variables for NaNs
  256. ActiveThorns = "NaNChecker"
  257. NaNChecker::check_every = 4096
  258. NaNChecker::action_if_found = "terminate"
  259. NaNChecker::check_vars = "
  260. ADMBase::metric
  261. ADMBase::lapse
  262. ADMBase::shift
  263. IllinoisGRMHD::rho_b
  264. IllinoisGRMHD::P
  265. smallbPoynET::smallb2
  266. ML_BSSN::ML_log_confac
  267. "
  268.  
  269. ActiveThorns = "smallbPoynET"
  270.  
  271. smallbPoynET::smallbPoynET_compute_every = 1024
  272.  
  273. ActiveThorns = "ADMMass"
  274.  
  275. ADMMass::ADMMass_compute_every = 1024
  276. ADMMass::ADMMass_use_all_volume_as_volume_radius = "yes"
  277. ADMMass::ADMMass_distance_from_grid_boundary[0] = 1
  278. ADMMass::ADMMass_distance_from_grid_boundary[1] = 1
  279. ADMMass::ADMMass_distance_from_grid_boundary[2] = 1
  280.  
  281. ######################################################################
  282. # #
  283. # SPACETIME EVOLUTION #
  284. # #
  285. ######################################################################
  286.  
  287. ActiveThorns = "SymBase ADMBase InitBase ADMCoupling ADMMacros"
  288. ActiveThorns = "GenericFD NewRad"
  289. ActiveThorns = "ML_BSSN ML_BSSN_Helper"
  290. ActiveThorns = "SpaceMask CoordGauge Constants Vectors LoopControl"
  291.  
  292. SpaceMask::use_mask = "yes"
  293.  
  294.  
  295. ADMMacros::spatial_order = 4
  296.  
  297. ADMBase::evolution_method = "ML_BSSN"
  298. ADMBase::lapse_evolution_method = "ML_BSSN"
  299. ADMBase::shift_evolution_method = "ML_BSSN"
  300. ADMBase::dtlapse_evolution_method= "ML_BSSN"
  301. ADMBase::dtshift_evolution_method= "ML_BSSN"
  302.  
  303. # 1 + log and gamma driver gauge
  304. ML_BSSN::harmonicN = 1.0
  305. ML_BSSN::harmonicF = 2.0
  306.  
  307. ML_BSSN::ShiftGammaCoeff = 0.75
  308. # BetaDriver is the eta parameter. It is extremely important
  309. # to be able to resolve black holes.
  310. # It should be smaller than 1.5/max_time_step, (for CFL stability)
  311. # this can be achieved with time refinement factors, so to set
  312. # eta to ~1/M_ADM.
  313. ML_BSSN::BetaDriver = 0.10
  314.  
  315. # conformalaMethod = 1 for W, 0 for phi
  316. # ML_BSSN::conformalMethod = 1
  317.  
  318. ML_BSSN::MinimumLapse = 1.0e-8
  319. ML_BSSN::timelevels = 3
  320.  
  321. # Spacetime boundary conditions
  322. ML_BSSN::initial_boundary_condition = "extrapolate-gammas"
  323. ML_BSSN::rhs_boundary_condition = "NewRad"
  324. Boundary::radpower = 2
  325.  
  326. ActiveThorns = "Dissipation"
  327.  
  328. # Kreiss-Oliger dissipation to kill high-frequency modes
  329. Dissipation::epsdis = 0.01
  330. Dissipation::vars = "
  331. ML_BSSN::ML_log_confac
  332. ML_BSSN::ML_metric
  333. ML_BSSN::ML_curv
  334. ML_BSSN::ML_trace_curv
  335. ML_BSSN::ML_Gamma
  336. ML_BSSN::ML_lapse
  337. ML_BSSN::ML_shift
  338. "
  339.  
  340. ######################################################################
  341. # #
  342. # MATTER EVOLUTION #
  343. # #
  344. ######################################################################
  345.  
  346. ActiveThorns = "TmunuBase HydroBase LocalReduce AEILocalInterp LocalInterp"
  347. ActiveThorns = "IllinoisGRMHD"
  348.  
  349. # IllinoisGRMHD requires staggered operators
  350. CarpetLib::support_staggered_operators = "yes"
  351.  
  352. TmunuBase::stress_energy_storage = "yes"
  353. TmunuBase::stress_energy_at_RHS = "yes"
  354. TmunuBase::prolongation_type = "none"
  355.  
  356. HydroBase::timelevels = 3
  357.  
  358. # Psi6threshold activates tricks to deplete black holes of matter to
  359. # stabilize the evolution
  360. IllinoisGRMHD::Psi6threshold = 1e2
  361. # IllinoisGRMHD::gamma_th = 1.3333333333
  362. IllinoisGRMHD::gamma_th = 1.325
  363. IllinoisGRMHD::neos = 1
  364.  
  365. # Lorenz damp should be around 8/M_ADM of the system
  366. # BUT 1.5/damp_lorenz should be greater than the
  367. # maximum time step
  368. IllinoisGRMHD::damp_lorenz = 1.1
  369.  
  370. # Tau atm should be in this case 4rho
  371. IllinoisGRMHD::tau_atm = 4e-16
  372. IllinoisGRMHD::rho_b_atm = 1e-16
  373. IllinoisGRMHD::verbose = "no"
  374.  
  375. ######################################################################
  376. # #
  377. # INITIAL DATA #
  378. # #
  379. ######################################################################
  380.  
  381. ActiveThorns = "ID_converter_ILGRMHD Convert_to_HydroBase Seed_Magnetic_Fields"
  382.  
  383. Carpet::init_3_timelevels = "no"
  384. Carpet::init_fill_timelevels = "yes"
  385.  
  386.  
  387. ADMBase::initial_data = "cook"
  388. ADMBase::initial_lapse = "cook"
  389. ADMBase::initial_shift = "cook"
  390. ADMBase::initial_dtlapse = "zero"
  391. ADMBase::initial_dtshift = "zero"
  392.  
  393. ActiveThorns = "CookImporter"
  394. CookImporter::rns_dat_file = "/u/sciteam/bozzola/Cactus/rns_B.dat"
  395. CookImporter::Poly_Gamma = 1.3333333333
  396. CookImporter::rho_0_atm = 1e-16
  397.  
  398. # IL_RNSID::save_2Dmodel = "yes"
  399. # IL_RNSID::recover_2Dmodel = "yes"
  400. # # il_rnsid::model2D_file = "/u/sciteam/bozzola/et_old/Cactus/model2D.dat"
  401.  
  402. # IL_RNSID::accuracy = 1e-10
  403. # IL_RNSID::RNS_rho_min = 5e-15
  404. # IL_RNSID::RNS_Gamma = 1.3333333333
  405. # IL_RNSID::RNS_K = 1.0
  406. # IL_RNSID::rho_central = 3.38e-6
  407. # IL_RNSID::axes_ratio = 0.24
  408. # IL_RNSID::rotation_type = "diff"
  409. # IL_RNSID::A_diff = 0.333333333333
  410. # # IL_RNSID::A_2_diff = 10
  411. # # IL_RNSID::b_diff = 0.3
  412. # IL_RNSID::add_perturbation = "yes"
  413. # IL_RNSID::pert_method = "density"
  414. # IL_RNSID::pert_type = "cylindrical"
  415. # IL_RNSID::pert_m = 2
  416. # IL_RNSID::pert_amplitude = 60e-5
  417.  
  418. ID_converter_ILGRMHD::Gamma_Initial = 1.3333333333
  419. ID_converter_ILGRMHD::K_Initial = 1.0
  420. # After the first restart with Reset_Gamma = 1 the value of gamma_th
  421. # in the par file has to be changed to match Gamma_th_new, otherwise
  422. # Cactus would complain for having steered a fixed parameter
  423. ID_converter_ILGRMHD::Reset_Gamma = 1
  424. ID_converter_ILGRMHD::Gamma_th_new = 1.325
  425. # It is useful to break the symmetry by a tiny amount
  426. # This helps rhobfinder
  427. ID_converter_ILGRMHD::random_pert = 1e-12
  428. ID_converter_ILGRMHD::pert_only_rho = "yes"
  429. # ID_converter_ILGRMHD::pure_hydro = "yes"
  430.  
  431. Seed_Magnetic_Fields::enable_IllinoisGRMHD_staggered_A_fields = "yes"
  432. Seed_Magnetic_Fields::A_b = 0 # Should yield 1/beta = b^2/(2P) ~ 0.001
  433. Seed_Magnetic_Fields::P_cut = 0
  434.  
  435. ######################################################################
  436. # #
  437. # CHECKPOINT/RECOVER #
  438. # #
  439. ######################################################################
  440.  
  441. ActiveThorns = "IOUtil"
  442.  
  443. IOUtil::checkpoint_dir = $parfile
  444. IOUtil::checkpoint_every_walltime_hours = 12
  445. IOUtil::checkpoint_keep = 1
  446. IOUtil::checkpoint_on_terminate = "yes"
  447. IOUtil::out_dir = $parfile
  448. IOUtil::out_fileinfo = "none"
  449. IOUtil::parfile_write = "generate"
  450. IOUtil::recover = "autoprobe"
  451. IOUtil::recover_dir = $parfile
  452. IOHDF5::checkpoint = "yes"
  453.  
  454. ActiveThorns = "TerminationTrigger"
  455.  
  456. TerminationTrigger::max_walltime = 12.0
  457. TerminationTrigger::on_remaining_walltime = 30 # minutes
  458. TerminationTrigger::termination_from_file = "yes"
  459. TerminationTrigger::create_termination_file = "yes"
  460. TerminationTrigger::termination_file = "../TERMINATE"
  461.  
  462.  
  463. ######################################################################
  464. # #
  465. # OUTPUT #
  466. # #
  467. ######################################################################
  468.  
  469. ActiveThorns = "CarpetIOASCII CarpetIOScalar CarpetIOHDF5 CarpetIOBasic"
  470.  
  471. ActiveThorns = "TimerReport"
  472.  
  473. TimerReport::out_every = 32768
  474. TimerReport::out_filename = "TimerReport"
  475. TimerReport::output_all_timers = "yes"
  476. TimerReport::output_all_timers_together = "yes"
  477. TimerReport::output_all_timers_readable = "yes"
  478. TimerReport::n_top_timers = 40
  479.  
  480. Carpet::output_xml_timer_tree = "yes"
  481. Carpet::output_timer_tree_every = 32768
  482.  
  483. ActiveThorns = "SystemStatistics"
  484.  
  485.  
  486. # STDOUT
  487.  
  488. IOBasic::outInfo_every = 512
  489. IOBasic::outInfo_vars = "
  490. Carpet::physical_time_per_hour
  491. IllinoisGRMHD::rho_b
  492. ADMBase::alp
  493. # smallbPoynET::smallb2
  494. # ML_BSSN::H
  495. SystemStatistics::maxrss_mb
  496. "
  497.  
  498. # 2D HDF5
  499.  
  500. CarpetIOHDF5::out2D_every = 1024
  501. CarpetIOHDF5::out2D_xy = "yes"
  502. CarpetIOHDF5::out2D_xz = "yes"
  503. CarpetIOHDF5::out2D_yz = "no"
  504. CarpetIOHDF5::out2D_dir = "${parfile}/data/H5_2d"
  505. CarpetIOHDF5::out2D_vars = "
  506. IllinoisGRMHD::rho_b
  507. IllinoisGRMHD::P
  508. ADMBase::alp
  509. # smallbPoynET::smallb2
  510. # smallbPoynET::smallbx
  511. # smallbPoynET::smallby
  512. # smallbPoynET::smallbz
  513. WeylScal4::Psi4r
  514. WeylScal4::Psi4i
  515. "
  516.  
  517. CarpetIOHDF5::out3D_every = 2048
  518. CarpetIOHDF5::out3D_ghosts = no
  519. CarpetIOHDF5::out3D_outer_ghosts = no
  520. CarpetIOHDF5::out3D_dir = "${parfile}/data/H5_3d"
  521. CarpetIOHDF5::out3D_vars = "
  522. IllinoisGRMHD::rho_b
  523. IllinoisGRMHD::P
  524. # smallbPoynET::smallb2
  525. IllinoisGRMHD::vx
  526. IllinoisGRMHD::vy
  527. IllinoisGRMHD::vz
  528. # smallbPoynET::smallbx
  529. # smallbPoynET::smallby
  530. # smallbPoynET::smallbz
  531. WeylScal4::Psi4r
  532. WeylScal4::Psi4i
  533. "
  534.  
  535. # Scalar ASCII
  536.  
  537. CarpetIOScalar::outScalar_every = 1024
  538. CarpetIOScalar::all_reductions_in_one_file = "no"
  539. CarpetIOScalar::outScalar_reductions = "minimum maximum average"
  540. CarpetIOScalar::outScalar_vars = "
  541. IllinoisGRMHD::rho_b
  542. IllinoisGRMHD::P
  543. ADMBase::alp
  544. # smallbPoynET::smallb2{reductions = 'minimum maximum norm1 norm2 average'}
  545. # smallbPoynET::b2Pratio{reductions = 'maximum'}
  546. ML_BSSN::ML_Ham{reductions = 'minimum maximum norm1 norm2'}
  547. ML_BSSN::ML_Mom{reductions = 'minimum maximum norm1 norm2'}
  548. SystemStatistics::maxrss_mb{out_every = 128}
  549. SystemStatistics::swap_used_mb
  550. "
  551.  
  552. CarpetIOASCII::out0D_every = 1024
  553. CarpetIOASCII::out0D_vars = "
  554. Carpet::timing
  555. ADMMass::ADMMass_VolumeMass[0]
  556. QuasiLocalMeasures::qlm_scalars
  557. SphericalSurface::sf_origin
  558. "
  559.  
  560. # IOASCII::out1D_every = 2
  561. # IOASCII::one_file_per_group = no
  562. # IOASCII::output_symmetry_points = no
  563. # IOASCII::out3D_ghosts = no
  564. # IOASCII::out3D_outer_ghosts = no
  565. # IOASCII::out1D_vars = "
  566. # IllinoisGRMHD::rho_b
  567. # "
  568.  
  569.  
  570. ######################################################################
  571. # #
  572. # TRIGGERS #
  573. # #
  574. ######################################################################
  575.  
  576.  
  577. ActiveThorns = "Trigger"
  578. Trigger::Trigger_Number = 8
  579.  
  580. Trigger::Trigger_Once [0] = 1
  581. Trigger::Trigger_Checked_Variable [0] = "IllinoisGRMHD::rho_b"
  582. Trigger::Trigger_Reduction [0] = "maximum"
  583. Trigger::Trigger_Relation [0] = ">"
  584. Trigger::Trigger_Checked_Value [0] = 1e-5
  585. Trigger::Trigger_Reaction [0] = "steerscalar"
  586. Trigger::Trigger_Steered_Scalar [0] = "CarpetRegrid2::num_levels[1]"
  587. Trigger::Trigger_Steered_Scalar_Value[0] = "9"
  588.  
  589. Trigger::Trigger_Once [1] = 1
  590. Trigger::Trigger_Checked_Variable [1] = "IllinoisGRMHD::rho_b"
  591. Trigger::Trigger_Reduction [1] = "maximum"
  592. Trigger::Trigger_Relation [1] = ">"
  593. Trigger::Trigger_Checked_Value [1] = 1e-5
  594. Trigger::Trigger_Reaction [1] = "steerscalar"
  595. Trigger::Trigger_Steered_Scalar [1] = "CarpetRegrid2::num_levels[0]"
  596. Trigger::Trigger_Steered_Scalar_Value[1] = "9"
  597.  
  598. Trigger::Trigger_Once [2] = 1
  599. Trigger::Trigger_Checked_Variable [2] = "IllinoisGRMHD::rho_b"
  600. Trigger::Trigger_Reduction [2] = "maximum"
  601. Trigger::Trigger_Relation [2] = ">"
  602. Trigger::Trigger_Checked_Value [2] = 5e-5
  603. Trigger::Trigger_Reaction [2] = "steerscalar"
  604. Trigger::Trigger_Steered_Scalar [2] = "CarpetRegrid2::num_levels[1]"
  605. Trigger::Trigger_Steered_Scalar_Value[2] = "10"
  606.  
  607. Trigger::Trigger_Once [3] = 1
  608. Trigger::Trigger_Checked_Variable [3] = "IllinoisGRMHD::rho_b"
  609. Trigger::Trigger_Reduction [3] = "maximum"
  610. Trigger::Trigger_Relation [3] = ">"
  611. Trigger::Trigger_Checked_Value [3] = 5e-5
  612. Trigger::Trigger_Reaction [3] = "steerscalar"
  613. Trigger::Trigger_Steered_Scalar [3] = "CarpetRegrid2::num_levels[0]"
  614. Trigger::Trigger_Steered_Scalar_Value[3] = "10"
  615.  
  616. Trigger::Trigger_Once [4] = 1
  617. Trigger::Trigger_Checked_Variable [4] = "IllinoisGRMHD::rho_b"
  618. Trigger::Trigger_Reduction [4] = "maximum"
  619. Trigger::Trigger_Relation [4] = ">"
  620. Trigger::Trigger_Checked_Value [4] = 1e-4
  621. Trigger::Trigger_Reaction [4] = "steerscalar"
  622. Trigger::Trigger_Steered_Scalar [4] = "CarpetRegrid2::num_levels[1]"
  623. Trigger::Trigger_Steered_Scalar_Value[4] = "11"
  624.  
  625. Trigger::Trigger_Once [5] = 1
  626. Trigger::Trigger_Checked_Variable [5] = "IllinoisGRMHD::rho_b"
  627. Trigger::Trigger_Reduction [5] = "maximum"
  628. Trigger::Trigger_Relation [5] = ">"
  629. Trigger::Trigger_Checked_Value [5] = 1e-4
  630. Trigger::Trigger_Reaction [5] = "steerscalar"
  631. Trigger::Trigger_Steered_Scalar [5] = "CarpetRegrid2::num_levels[0]"
  632. Trigger::Trigger_Steered_Scalar_Value[5] = "11"
  633.  
  634. # if lapse < 0.2 enable AHFinderDirect
  635. Trigger::Trigger_Once [6] = 1
  636. Trigger::Trigger_Checked_Variable [6] = "ADMBase::alp"
  637. Trigger::Trigger_Reduction [6] = "minimum"
  638. Trigger::Trigger_Relation [6] = "<"
  639. Trigger::Trigger_Checked_Value [6] = 0.25
  640. Trigger::Trigger_Reaction [6] = "steerparam"
  641. Trigger::Trigger_Steered_Parameter_Thorn[6] = "AHFinderDirect"
  642. Trigger::Trigger_Steered_Parameter_Name [6] = "find_every"
  643. Trigger::Trigger_Steered_Parameter_Value[6] = "128"
  644.  
  645. # When an horizon is found use its centroid to track
  646. Trigger::Trigger_Once [7] = 1
  647. Trigger::Trigger_Checked_Variable [7] = "AHFinderDirect::ah_found_flag[1]"
  648. Trigger::Trigger_Relation [7] = "=="
  649. Trigger::Trigger_Checked_Value [7] = 1
  650. Trigger::Trigger_Reaction [7] = "steerparam"
  651. Trigger::Trigger_Steered_Parameter_Thorn[7] = "NSTracker"
  652. Trigger::Trigger_Steered_Parameter_Name [7] = "NSTracker_use_AHFinderDirect"
  653. Trigger::Trigger_Steered_Parameter_Value[7] = "true"
  654.  
  655.  
  656. ######################################################################
  657. # #
  658. # HORIZON DETECTION #
  659. # #
  660. ######################################################################
  661.  
  662. ActiveThorns = "AHFinderDirect"
  663.  
  664. # This is later steered
  665. AHFinderDirect::find_every = 0
  666.  
  667. AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial interpolation"
  668. AHFinderDirect::geometry_interpolator_pars = "order=4"
  669. AHFinderDirect::surface_interpolator_name = "Lagrange polynomial interpolation"
  670. AHFinderDirect::surface_interpolator_pars = "order=4"
  671.  
  672. AHFinderDirect::output_h_every = 512
  673.  
  674. AHFinderDirect::N_horizons = 3
  675.  
  676. AHFinderDirect::max_Newton_iterations__initial = 50
  677. AHFinderDirect::max_Newton_iterations__subsequent = 50
  678. AHFinderDirect::max_allowable_Theta_growth_iterations = 10
  679. AHFinderDirect::max_allowable_Theta_nonshrink_iterations = 10
  680. AHFinderDirect::verbose_level = "physics details"
  681. AHFinderDirect::move_origins = yes
  682. AHFinderDirect::predict_origin_movement = yes
  683.  
  684.  
  685. AHFinderDirect::initial_guess_method[1] = "coordinate sphere"
  686. AHFinderDirect::reset_horizon_after_not_finding [1] = "no"
  687. AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0
  688. AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 10
  689. AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 10
  690. AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 0
  691. AHFinderDirect::origin_x [1] = 10
  692. AHFinderDirect::origin_y [1] = 10
  693. AHFinderDirect::origin_z [1] = 0
  694. AHFinderDirect::which_surface_to_store_info_by_name [1] = "BH Right"
  695.  
  696. AHFinderDirect::initial_guess_method[2] = "coordinate sphere"
  697. AHFinderDirect::reset_horizon_after_not_finding [2] = "no"
  698. AHFinderDirect::initial_guess__coord_sphere__radius [2] = 2.0
  699. AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -10
  700. AHFinderDirect::initial_guess__coord_sphere__y_center[2] = -10
  701. AHFinderDirect::initial_guess__coord_sphere__z_center[2] = 0
  702. AHFinderDirect::origin_x [2] = -10
  703. AHFinderDirect::origin_y [2] = -10
  704. AHFinderDirect::origin_z [2] = 0
  705. AHFinderDirect::which_surface_to_store_info_by_name [2] = "BH Left"
  706.  
  707. AHFinderDirect::initial_guess_method[3] = "coordinate sphere"
  708. AHFinderDirect::origin_y [3] = 0
  709. AHFinderDirect::origin_z [3] = 0
  710. AHFinderDirect::origin_x [3] = 0
  711. AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0
  712. AHFinderDirect::initial_guess__coord_sphere__radius [3] = 1.0
  713. AHFinderDirect::which_surface_to_store_info [3] = 4
  714. AHFinderDirect::set_mask_for_individual_horizon [3] = no
  715.  
  716. # ActiveThorns = "RestMass"
  717. # RestMass::RestMass_compute_every = 1
  718.  
  719. ######################################################################
  720. # #
  721. # GRAVITATIONAL WAVE EXTRACTION #
  722. # #
  723. ######################################################################
  724.  
  725. # Wave exraction
  726.  
  727. ActiveThorns="WeylScal4 Multipole"
  728.  
  729. WeylScal4::offset = 1e-8
  730. WeylScal4::fd_order = "4th"
  731.  
  732. # Corresponing to 1/r [0.001 , 0.0013, 0.0016, 0.0019, 0.0022, 0.0025] / 0.8
  733. Multipole::nradii = 6
  734. Multipole::out_every = 1024
  735. Multipole::radius[0] = 500
  736. Multipole::radius[1] = 568.18125
  737. Multipole::radius[2] = 657.895
  738. Multipole::radius[3] = 781.25
  739. Multipole::radius[4] = 961.53875
  740. Multipole::radius[5] = 1250
  741. Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='Psi4'}"
  742. Multipole::l_max = 6
  743.  
  744. ActiveThorns = "QuasiLocalMeasures"
  745.  
  746. QuasiLocalMeasures::num_surfaces = 3
  747. QuasiLocalMeasures::spatial_order = 4
  748. QuasiLocalMeasures::interpolator = "Lagrange polynomial interpolation"
  749. QuasiLocalMeasures::interpolator_options = "order=4"
  750. QuasiLocalMeasures::surface_index [0] = 2
  751. QuasiLocalMeasures::surface_index [1] = 3
  752. QuasiLocalMeasures::surface_index [2] = 4
  753.  
  754. ######################################################################
  755. # #
  756. # PARERGA ET PARALIPOMENA #
  757. # #
  758. ######################################################################
  759.  
  760. # ActiveThorns="HTTPD HTTPDExtra Socket"
  761.  
  762. # HTTPD::user = "cactus"
  763. # HTTPD::password = "einstein"
  764.  
  765.  
  766. ### Local Variables:
  767. ### mode: sh
  768. ### End:
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!