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 14:02]
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   +
-   +
-  #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                         = bo newfmr move
-  amr_fmr                     = nestedboxes +  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  +  grid                         = box bitant 
-  of the bitant+
      
   ############################################################################   ############################################################################
Line 75: Line 70:
      
   hrsc_nghosts = 4 # how many ghost points are used to go from courser to finer level   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_rec = WENOZ  
-  hrsc_TVD_limiter = MC2 # flux limiter | hr high resolution sc = shock capturing scheme +  hrsc_TVD_limiter = MC2  
-  hrsc_rec_metric = LAG4 # +  hrsc_rec_metric = LAG4 
   hrsc_flux = LLF # HLL, HLLC - Riemann solvers   hrsc_flux = LLF # HLL, HLLC - Riemann solvers
-  hrsc_rec_HO = WENOZ # HO = higher order +  hrsc_rec_HO = WENOZ 
      
-  conservative_amr = yes # conservative adaptive mesh refinement, in order to working, below needs to be set to linear  +  conservative_amr = yes
   matter_interpolate_scheme_restriction  = linear   matter_interpolate_scheme_restriction  = linear
   matter_interpolate_scheme_prolongation = linear   matter_interpolate_scheme_prolongation = linear
   matter_interpolate_order = 4   matter_interpolate_order = 4
      
-  # conservative (D=density with lorentz factor, tau= , S= kinda Momentum) to primative variables (  +  grhd_C2P = p_root 
-  grhd_C2P = p_root # proot = newton raphson calculation+
   grhd_C2P_NewtonRaphsonTR = 1e-8    grhd_C2P_NewtonRaphsonTR = 1e-8 
   grhd_C2P_NewtonRaphsonNR = 100   grhd_C2P_NewtonRaphsonNR = 100
Line 97: Line 91:
      
   #############################################################################   #############################################################################
-  # evolution, RK = runge kutta+  # evolution
      
   evolution_method        = rk   evolution_method        = rk
   evolution_method_rk     = rk4   evolution_method_rk     = rk4
-  dtfac                   = 0.25    # courant factor, relating space and time, causal +  dtfac                   = 0.25   
-  finaltime = 20000   # end of simulation time in solar masses+  finaltime = 20000  
      
   z4_normalizedetg            = yes   z4_normalizedetg            = yes
Line 112: Line 106:
   z4_shift                    = gamma0 withGadv withBadv withShiftadv   z4_shift                    = gamma0 withGadv withBadv withShiftadv
   z4_shiftalphapower          = 0.0   z4_shiftalphapower          = 0.0
-  z4_shiftdriver              = 0.3 # should be checked, depending on mass of the system+  z4_shiftdriver              = 0.3
   z4_shiftgammacoeff          = 0.75   z4_shiftgammacoeff          = 0.75
   z4_kappa1                   = 0.02   z4_kappa1                   = 0.02
Line 133: Line 127:
      
   #############################################################################   #############################################################################
-  # invariants for extracting waveform+  # invariants 
      
   ntheta                      = 47   ntheta                      = 47
Line 155: Line 149:
      
   #######################################################   #######################################################
-  # (better use AHmod !)+  # AHmod
      
   AHmod_verbose     = no   AHmod_verbose     = no
Line 162: Line 156:
   AHmod_LMAX        = 10   AHmod_LMAX        = 10
   AHmod_interpolation_order = 4   AHmod_interpolation_order = 4
-  AHmod_nhorizons   = 1 # has to be 1, as we just have 1 BH horizon+  AHmod_nhorizons   = 1
   AHmod_searchMTS   = 1 600  10000.0  0   AHmod_searchMTS   = 1 600  10000.0  0
   AHmod_uselast     = yes   AHmod_uselast     = yes
Line 172: Line 166:
   AHmod_mass_tol    = 1.0e-04   AHmod_mass_tol    = 1.0e-04
   AHmod_hmean_tol   = 100.0   AHmod_hmean_tol   = 100.0
-  #AHmod_time        = 4.0/1.5+  AHmod_time        = 4.0/1.5
   AHmod_output      = yes   AHmod_output      = yes
   AHmod_output_xyt  = yes   AHmod_output_xyt  = yes
Line 181: Line 175:
      
   hydroanalysis_ejecta_spheres  = yes   hydroanalysis_ejecta_spheres  = yes
-  hydroanalysis_ejecta_spheres_radius = 200 # constructing diff spheres at which you check if matter is inside, checking ejecta  +  hydroanalysis_ejecta_spheres_radius = 200  
-  passing through shells/spheres +  hydroanalysis_ejecta_nradius = 5  
-  hydroanalysis_ejecta_nradius = 5 # number of spheres, where we check for ejecta matter +  hydroanalysis_ejecta_dradius = 100
-  hydroanalysis_ejecta_dradius = 100 # stepsize of ejecta sphere radius increase +
      
   hydroanalysis_Mbar_radius  = 8    hydroanalysis_Mbar_radius  = 8 
Line 206: Line 199:
   # output   # 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    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   hydroa_Db hydroa_fluid2_Du hydroa_fluid2_Db
Line 216: Line 208:
   #1doutputall             = yes   #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    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   grhd_fluid2_vx grhd_fluid2_v2
Line 223: Line 214:
   2doutputall             = yes   2doutputall             = yes
      
-  3douttime               = 90 
   3doutput                = grhd_rho grhd_fluid2_rho hydroa_Du hydroa_fluid2_Du alpha rpsi4   3doutput                = grhd_rho grhd_fluid2_rho hydroa_Du hydroa_fluid2_Du alpha rpsi4
   3dformat                = vtk binary float   3dformat                = vtk binary float
Line 229: Line 219:
   3doutputall = yes   3doutputall = yes
              
-   # For now the string for the 2doutputr parameter is not used, just make sure, it is non-empty! 
   #2doutputr               = sphere_data   #2doutputr               = sphere_data
   #############################################################################   #############################################################################
  
  
Last modified: le 2024/07/23 14:02

Page Tools