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
grhd_C2P_NewtonRaphsonNR = 100
grhd_vmax = 0.99
grhd_Wlor_max = 1e8
grhd_use_atmosphere = ColdStatic
grhd_atm_factor = 10
grhd_atm_level = 1e-11
#############################################################################
# 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 = 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
nphi = 46
invariants_compute_modes = yes
invariants_modes_r = 300 400 500 600 700 800 900 1000 1100 1200
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
0douttime = 10.
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
3douttime = 90
3doutput = grhd_rho grhd_fluid2_rho hydroa_Du hydroa_fluid2_Du alpha rpsi4
3dformat = vtk binary float
3doutinterpolate = no
3doutputall = yes
# For now the string for the 2doutputr parameter is not used, just make sure, it is non-empty!
#2doutputr = sphere_data
#############################################################################
Last modified: le 2024/07/23 14:02
