This is an old revision of the document!
An example par file
physics = adm z4 Gauge matter eos grhd2Fluid BNSdataReader Invariants AHmod hydroanalysis ############################################################################# # initial data BBID_modus = DSDS # ask Edo: ... dark star dark star, otherwise BNSS BNSdataReader_sgrid_exe = /zhome/academic/HLRS/xup/xupkuner/03_SGRID/08_reruns_2024/exe/sgrid_2F BNSdataReader_use_interpolator = yes # set to yes = for first time interpolation BNSdataReader_sgrid_datadir = "/lustre/hpe/ws10/ws10.1/ws/xupkuner- BAM_DMruns/00_NEW/01_3perc/01_BAM_R96/3perc_WT_1e10_M12_D16_SLY4C_1000MeV_EccRun_2_26_26_26" BNSdataReader_sgrid_simple_directory_structure = yes # see also: Ivans orga of sgrid folder BNSdataReader_Interpolate_verbose = no #yes for the first submission only, to see if sgrid to bam interpolation works BNSdataReader_keep_sgrid_output = no # yes for the time when it is doing the interpolation BNSdataReader_two_fluid_system = yes # for 2 Fluid data reading ############################################################################ # eos eos_ID = yes # bam would take eos table from sgrid and convert it to bam formalism eos1 = tab1dhot # eos1 refers to baryonic matter eos # hot stands for thermal corrections eos1_gamma = 1.75 # thermal coefficient, Gamma , # Gamma for thermal corrections eos2 = tab1d # eos_ID = yes, so the table is taken from sgrid ############################################################################ # checkpointing checkpoint = yes checkpoint_dt_hours = 5 checkpoint_dt_hours_quit = 23 # important to quite 1 hour before the clusters walltime checkpoint_variables = all ExitIfNAN = yes ############################################################################ # basic setup order_centered = 4 # finite differencing related order_advection = 4 advection_lopsided6 = 2 advection_lopsided = 1 order_dissipation = 6 dissipation_factor = 0.1 bampi_timer_on = yes # functionality to see in which parts of BAM it spends how many time bampi_lowlatency = no # should always be on no bampi_timer_reset_after_every_iteration = yes bampi_nghosts = 6 order_RP = 6 order_RP_shells = 6 amr_nbuffer = 6 amr_lmax = 7 # maximum number of refinement levels -1 amr_move_lcube = 2 # number of moving boxes / moving refinement levels
#spacing | always keep for numbers before the * fixed and change the Brueche nxyz = 128*2/2 #3/2 # non-moving box, normally 2 times number of amr_move_nxyz, amr_move_nxyz = 96*2/2 #3/2 # number of points across moving levels, related to resolution of BAM run # resolution of the actual Bam run | finest level dxyz = 36*2/2 #2/3 # spacing on the 0-th level | biggest level\grid amr_bo_dxmax = 6*2/2 #2/3 #
# first provide finest level information # can be calculated by using the diameter of a NS with a bit of extra space in geomtrical units and divide by number of moving boxes # then this calculated number is multiplied by a factor of 2 for each bigger level
amr = bo newfmr move amr_fmr = nestedboxes grid = box bitant #box quadrant, bitant (means the upper half of 3D space, from z [0, +inf] # quadrant should be one half of the bitant
############################################################################ # hydro
hrsc_nghosts = 4 # how many ghost points are used to go from courser to finer level hrsc_rec = WENOZ #reconstruction scheme, weighted non-oscillatory in z hrsc_TVD_limiter = MC2 # flux limiter | hr high resolution sc = shock capturing scheme hrsc_rec_metric = LAG4 # hrsc_flux = LLF # HLL, HLLC - Riemann solvers hrsc_rec_HO = WENOZ # HO = higher order
conservative_amr = yes # conservative adaptive mesh refinement, in order to working, below needs to be set to linear matter_interpolate_scheme_restriction = linear matter_interpolate_scheme_prolongation = linear matter_interpolate_order = 4
# conservative (D=density with lorentz factor, tau= , S= kinda Momentum) to primative variables ( grhd_C2P = p_root # proot = newton raphson calculation grhd_C2P_NewtonRaphsonTR = 1e-8 # threshold, accuracy grhd_C2P_NewtonRaphsonNR = 100 # maximum number of iteration
grhd_vmax = 0.99 # maximum velocity inside the code grhd_Wlor_max = 1e8 # maximum Lorentz factor
grhd_use_atmosphere = ColdStatic # how the atmosphere is modeled # belows params can be used to set the atmosphere inside the BAM grid grhd_atm_factor = 10 # take central rho and … see equations grhd_atm_level = 1e-11 # value
############################################################################# # evolution, RK = runge kutta
evolution_method = rk evolution_method_rk = rk4
dtfac = 0.25 # courant factor, relating space and time, causal finaltime = 20000 # end of simulation time in solar masses
z4_normalizedetg = yes z4_subtractA = yes z4_chi_div_floor = 1e-5
z4_initial_lapse = donothing z4_initial_shift = zero z4_lapse = 1+log withshift z4_shift = gamma0 withGadv withBadv withShiftadv z4_shiftalphapower = 0.0 z4_shiftdriver = 0.3 # should be checked, depending on mass of the system z4_shiftgammacoeff = 0.75
z4_kappa1 = 0.02 z4_kappa2 = 0.0
punctures_lapse = psiBL^(-2) # boundary
boundary = background radcentered
############################################################################ # gauge
Gauge = moving_puncture
compute_moving_puncture_distance = line #moving_puncture_fixz = none
#track_matter_var = alpha #track_matter_mode = min #track_matter_minmove = 0.01
############################################################################# # invariants for extracting waveform
ntheta = 47 # theta angle nphi = 46 # phi angle invariants_compute_modes = yes invariants_modes_r = 300 400 500 600 700 800 900 1000 1100 1200 # radii at which we want to extract the invariants (dE/dt –> GWs) invariants_modes_lmin = 0 invariants_modes_lmax = 5 invariants_energy_r = 300 400 500 600 700 800 900 1000 1100 1200 gauss_codacci_mainardi = standard invariants_order = 6
############################################################################# #ADM Mass
ADM_mass_ncircles = 101 ADM_mass_npoints = 80 ADM_mass_lmin = 1 ADM_mass_lmax = 5 ADM_mass_r = 300 400 500 600 700 800 900 1000 1100 1200
####################################################### # (better use AHmod !)
AHmod_verbose = no AHmod_ntheta = 30 AHmod_nphi = 60 AHmod_LMAX = 10 AHmod_interpolation_order = 4 AHmod_nhorizons = 1 # has to be 1, as we just have 1 BH horizon AHmod_searchMTS = 1 600 10000.0 0 AHmod_uselast = yes AHmod_initial_guess_expand = 1.0 AHmod_LevelOffset = 1 AHmod_UseOptimalLevel = no AHmod_surface0_WaitUntilClosePunctures = no AHmod_flow_iter = 500 AHmod_mass_tol = 1.0e-04 AHmod_hmean_tol = 100.0 #AHmod_time = 4.0/1.5 AHmod_output = yes AHmod_output_xyt = yes AHmod_output_lm = yes
############################################################################## # hydroanalysis
hydroanalysis_ejecta_spheres = yes hydroanalysis_ejecta_spheres_radius = 200 # constructing diff spheres at which you check if matter is inside, checking ejecta passing through shells/spheres hydroanalysis_ejecta_nradius = 5 # number of spheres, where we check for ejecta matter hydroanalysis_ejecta_dradius = 100 # stepsize of ejecta sphere radius increase
hydroanalysis_Mbar_radius = 8 hydroanalysis_Mbar_nradius = 12 hydroanalysis_Mbar_dradius = 1. hydroanalysis_rATM = 1e-25 hydroa_mode_projection = no
############################################################################ #puncture properties
moving_puncture_finboxfix = 20 moving_puncture_finboxfixv = 6
moving_puncture_properties = yes puncture_properties_r = 8 9 10 11 12 14 15 20 puncture_properties_punc = 2 puncture_properties_npoints = 40 puncture_properties_circles = 41
############################################################################# # output
#0doutiter = 1 0douttime = 10. # output time needs to be calculated and is depending on the resolution quiet carefully | need to coincide with iteration 0doutput = alpha ham momx grhd_D grhd_rho grhd_p grhd_fluid2_D grhd_fluid2_rho grhd_fluid2_p momy momz hydroa_Du hydroa_Db hydroa_fluid2_Du hydroa_fluid2_Db
#1doutiter = 1 #1douttime = 2 #1doutput = alpha betax ham momx grhd_rho flagprolong rpsi4 flagrestrict grhd_fluid2_rho #1doutinterpolate = no #1doutputall = yes
2douttime = 4 2doutput = alpha betax ham momx grhd_D grhd_rho grhd_p grhd_vx grhd_v2 grhd_fluid2_D grhd_fluid2_rho grhd_fluid2_p grhd_fluid2_vx grhd_fluid2_v2 2dformat = vtk binary float 2doutinterpolate = no 2doutputall = yes
# For now the string for the 2doutputr parameter is not used, just make # sure, it is non-empty! #2doutputr = sphere_data
3douttime = 90 3doutput = grhd_rho grhd_fluid2_rho hydroa_Du hydroa_fluid2_Du alpha rpsi4 3dformat = vtk binary float 3doutinterpolate = no 3doutputall = yes
#############################################################################
