This is an old revision of the document!
An example par file
physics = adm z4 Gauge matter eos grhdDY BNSdataReader Invariants ADM_mass AHmod hydroanalysis nxyz = 128*2 amr_move_nxyz = 64*2 dxyz = 16.128/2 # 0.252*2^6 amr_bo_dxmax = 0.252*8/2 0douttime = 0.252*4 2douttime = 0.252*4*30 3douttime = 0.252*4*150 Invariants_output_time = 0.252*4. AHmod_time = 0.252*32. ############################################################################ # initial data BNSdataReader_sgrid_exe = /zhome/academic/HLRS/xup/xuplgieg/sgrid_src_test/exe/sgrid BNSdataReader_sgrid_datadir = /zhome/academic/HLRS/xup/xuplgieg/ID/DD2_3D_q1_M27_b14 BNSdataReader_keep_sgrid_output = no BNSdataReader_use_interpolator = no BNSdataReader_IDfiles_dir = /lustre/hpe/ws10/ws10.1/ws/xuplgieg-BNS/ID_grid/DD2_4D BNSdataReader_reset_EOS = no #yes BNS_eos_Gamma = 1.75
############################################################################
# checkpointing checkpoint = yes checkpoint_dt_hours = 5 checkpoint_dt_hours_quit = 20 checkpoint_variables = all checkpoint_resub = yes ExitIfNAN = yes ############################################################################ # basic setup order_centered = 4 order_advection = 4 advection_lopsided6 = 2 advection_lopsided = 1 order_dissipation = 6 dissipation_factor = 0.5 dissipation_factor_level0 = 0.1 bampi_timer_on = yes bampi_lowlatency = yes bampi_nghosts = 6 order_RP = 6 order_RP_shells = 6 amr_nbuffer = 6 amr_lmax = 6 amr_move_lcube = 2 amr = bo newfmr move amr_fmr = nestedboxes grid = box bitant #for precessing cases this should be changed ############################################################################# # evolution evolution_method = rk evolution_method_rk = rk4g dtfac = 0.25 finaltime = 20000 #20000=0.098 sec.
order_timeinterpolation = 2 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 z4_shiftgammacoeff = 0.75 z4_kappa1 = 0.02 z4_kappa2 = 0.0 punctures_lapse = psiBL^(-2) ############################################################################# # eos eos = tab4D eos_tab_file = /lustre/hpe/ws10/ws10.1/ws/xuplgieg-BNS/eos/DD2_4D/DD2_4D_NEW_NEW.h5 eos_tab_file_beta = /lustre/hpe/ws10/ws10.1/ws/xuplgieg-BNS/eos/DD2_4D/DD2_4D_cold_beta_NEW.d eos_mb = 9.3046542925126903e+02 eos_T = 0.10000000000000001 eos_Y = 0.4 eos_interp = steffen ############################################################################# # hydro grhd_C2P = illinois_4d grhd_C2P_NewtonRaphsonTR = 1e-9 grhd_C2P_NewtonRaphsonNR = 200 grhd_vmax = 0.99 grhd_Wlor_max = 1e8 grhd_use_atmosphere = ColdStatic grhd_atm_factor = 10 grhd_atm_level = 10 #grhd_atm_T = 0. grhd_use_atmosphere_mask = yes grhd_recvel = bx grhd_use_nls = yes grhd_add_nls = yes nls_chi_thresh = 1e-04 nls_opt_scheme = RJS nls_num_sps = 5 nls_op_table = /lustre/hpe/ws10/ws10.1/ws/xuplgieg-BNS/eos/DD2_4D/opacities_DD2_4D_v07.h5 nls_op_table_beta = /lustre/hpe/ws10/ws10.1/ws/xuplgieg-BNS/eos/DD2_4D/opacities_DD2_4D_cold_beta_v07.d hrsc_nghosts = 4 hrsc_rec = WENOZ hrsc_TVD_limiter = MC2 hrsc_rec_metric = LAG6 hrsc_flux = HLL hrsc_rec_safety = no grhd_use_excision = no grhd_excision_rfct = 0.9 grhd_excision_modus = atm matter_interpolate_scheme_restriction = linear matter_interpolate_scheme_prolongation = linear matter_interpolate_order = 4 conservative_amr = yes camr_treshold = 1e20 rhs_matter_null_ONlevel0 = yes ################################################################################ # M1
grrhd_m1_implicit_sources = linear grrhd_m1_fluxes = Radice
grrhd_m1_C2P_NewtonRaphsonTR = 1e-8 grrhd_m1_C2P_NewtonRaphsonNR = 100
rad_hrsc_rec = CONST rad_hrsc_flux = LLF
grrhd_m1_decouple_fac = 1e5
radiation_use_mask = no optical_depth_mask = 100
grrhd_m1_NuHeat = yes grrhd_m1_nls_rates = no
radiation_interpolate_order = 4 radiation_interpolate_scheme_restriction = LINTVD radiation_interpolate_scheme_prolongation = LINTVD
rhs_radiation_null_ONlevel0 = yes
radiation_at_merger = no radiation_activate_dist = 18
rad_evolution_method = rk rad_evolution_method_rk = rk3
############################################################################# # output
stdout_flush = yes
0doutput = alpha adm_rho ham momx grhd_rho grhd_D grhd_DY grhd_Y grhd_Tau grhd_T grhd_p hydroa_Du hydroa_ut grhd_Ym grhd_DYm
nls_Qeff_nue nls_Reff_nue nls_lum_nue nls_lum_nua nls_lum_nux nls_Qeff_nua nls_Reff_nua nls_Qeff_nux nls_Reff_nux nls_lum0_nue nls_lum0_nua nls_lum0_nux nls_Qeff_num nls_Reff_num nls_Qeff_numa nls_Reff_numa nls_lum_num nls_lum_numa nls_lum0_num nls_lum0_numa
#1doutput = grhd_rho grhd_Y grhd_T grhd_p grhd_epsl
2doutput = grhd_rho grhd_p grhd_vx grhd_v2 grhd_Y grhd_Ym grhd_DYm grhd_D grhd_T grhd_DY hydroa_ut hydroa_Du hydroa_DYu
nls_Qeff_nue nls_Qeff_nua nls_Qeff_nux nls_Reff_nua nls_Reff_nue nls_Reff_nux nls_Qeff_num nls_Reff_num nls_Qeff_numa nls_Reff_numa nls_tau1_nue nls_tau1_nua nls_tau1_nux nls_tau1_num nls_tau1_numa nls_tau0_nue nls_tau0_nua nls_tau0_nux nls_tau0_num nls_tau0_numa
2doutputmode = xy xz 2dformat = xdmf binary float 2doutinterpolate = no 2doutputall = yes 2doutputr = rpsi4 ipsi4 2doutputr_type = SphericalDF
3doutput = grhd_rho grhd_Y grhd_T grhd_D hydroa_ut hydroa_Du grhd_vx grhd_Ym grhd_vy grhd_vz 3doutputall = yes
############################################################################# # 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
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 # (extract same radii as waves)
ADM_mass_ncircles = 101 ADM_mass_npoints = 80 ADM_mass_lmin = 0 ADM_mass_lmax = 5 ADM_mass_r = 300 400 500 600 700 800 900 1000 1100 1200
####################################################### # AH finder
AHmod_verbose = no AHmod_ntheta = 30 AHmod_nphi = 60 AHmod_LMAX = 10 AHmod_interpolation_order = 4 AHmod_nhorizons = 1 AHmod_searchMTS = 1 2000 20000.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_output = yes AHmod_output_xyt = yes AHmod_output_lm = yes
####################################################### #hydroanalysis
hydroanalysis_ejecta_spheres = yes hydroanalysis_ejecta_spheres_radius = 200 hydroanalysis_ejecta_nradius = 5 hydroanalysis_ejecta_dradius = 100
hydroanalysis_radiation_spheres = yes
hydroanalysis_ejecta_angular = yes hydroanalysis_ejecta_angular_radius = 300
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
