Advertisement
Guest User

Untitled

a guest
Jun 12th, 2019
211
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  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
RAW Paste Data Copied
Advertisement