Differences

This shows you the differences between two versions of the page.

Link to this comparison view

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

Page Tools