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:48]
theoastro 		bam:example_parfile [2024/07/23 14:16] (current)
theoastro
Line 1: Line 1:
 === An example par file ===  === An example par file === 
  
-  physics = adm z4 Gauge matter eos grhdDY BNSdataReader Invariants ADM_mass AHmod hydroanalysis+  physics = adm z4 Gauge matter eos grhd2Fluid BNSdataReader Invariants AHmod hydroanalysis 
      
-  nxyz                        = 128*2 +  nxyz = 128*2/2   
-  amr_move_nxyz               = 64*2 +  amr_move_nxyz               = 96*2/2    
-  dxyz                        = 16.128/2                           # 0.252*2^6 +  dxyz                        = 36*2/2   
-  amr_bo_dxmax                = 0.252*8/2 +  amr_bo_dxmax                = 6*2/2      
-  0douttime                   = 0.252*4 +   
-  2douttime                   = 0.252*4*30 +  0douttime               = 10. 
-  3douttime                   = 0.252*4*150+  2douttime               = 4 
 +  3douttime               = 90
   Invariants_output_time      = 0.252*4.   Invariants_output_time      = 0.252*4.
-  AHmod_time                  = 0.252*32.+     
 +  ############################################################################# 
 +  # initial data 
 +   
 +  BBID_modus = DSDS  
 +   
 +  BNSdataReader_sgrid_exe = /path/to/BAM/executable  
 +  BNSdataReader_use_interpolator = yes  
 +  BNSdataReader_sgrid_datadir = "/path/to/sgrid/folder/.../"  
 +  BNSdataReader_sgrid_simple_directory_structure = yes  
 +  BNSdataReader_Interpolate_verbose = no 
 +  BNSdataReader_keep_sgrid_output = no  
 +  BNSdataReader_two_fluid_system = yes 
      
   ############################################################################   ############################################################################
-  # initial data+  # eos 
 +   
 +  eos_ID              = yes  
 +  eos1                = tab1dhot  
 +  eos1_gamma          = 1.75  
 +  eos2                = tab1d 
      
-  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   # checkpointing
      
   checkpoint                  = yes   checkpoint                  = yes
-  checkpoint_dt_hours         = 5 +  checkpoint_dt_hours         = 5  
-  checkpoint_dt_hours_quit    = 20 +  checkpoint_dt_hours_quit    = 23  
-  checkpoint_variables        = all  +  checkpoint_variables        = all 
-  checkpoint_resub            = yes +  ExitIfNAN                   = yes
-  ExitIfNAN     = yes+
      
   ############################################################################   ############################################################################
-  # basic setup+  # basic setup 
      
-  order_centered              = 4+  order_centered              = 4 
   order_advection             = 4   order_advection             = 4
   advection_lopsided6         = 2   advection_lopsided6         = 2
   advection_lopsided          = 1    advection_lopsided          = 1 
   order_dissipation           = 6   order_dissipation           = 6
-  dissipation_factor          = 0.5 +  dissipation_factor          = 0.1
-  dissipation_factor_level0   = 0.1+
      
-  bampi_timer_on              = yes +  bampi_timer_on              = yes  
-  bampi_lowlatency            = yes +  bampi_lowlatency            = no 
-  bampi_nghosts               = 6+  bampi_timer_reset_after_every_iteration = yes 
 +  bampi_nghosts               = 6 
   order_RP                    = 6   order_RP                    = 6
   order_RP_shells             = 6   order_RP_shells             = 6
-  amr_nbuffer                 = 6  +  amr_nbuffer                 = 6 
-  amr_lmax                    = 6 +  amr_lmax                    = 7 
-  amr_move_lcube              = 2 +  amr_move_lcube     = 2 
-  +
   amr                         = bo newfmr move   amr                         = bo newfmr move
-  amr_fmr                     = nestedboxes +  amr_fmr               = nestedboxes 
-  grid                        = box bitant #for precessing cases this should be changed +  grid                         = box bitant  
 +   
 +  ############################################################################ 
 +  # hydro 
 +   
 +  hrsc_nghosts = 4 # how many ghost points are used to go from courser to finer level 
 +  hrsc_rec = WENOZ  
 +  hrsc_TVD_limiter = MC2  
 +  hrsc_rec_metric = LAG4  
 +  hrsc_flux = LLF # HLL, HLLC - Riemann solvers 
 +  hrsc_rec_HO = WENOZ  
 +   
 +  conservative_amr = yes  
 +  matter_interpolate_scheme_restriction  = linear 
 +  matter_interpolate_scheme_prolongation = linear 
 +  matter_interpolate_order = 4 
 +   
 +  grhd_C2P = p_root  
 +  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  +  # evolution 
-  evolution_method            = rk +   
-  evolution_method_rk         = rk4g +  evolution_method        = rk 
-  dtfac                       = 0.25 +  evolution_method_rk     = rk4 
-  finaltime                   = 20000 #20000=0.098 sec. +  dtfac                   = 0.25   
- +  finaltime  = 20000  
-  order_timeinterpolation     = 2+
      
   z4_normalizedetg            = yes   z4_normalizedetg            = yes
Line 79: Line 110:
   z4_kappa1                   = 0.02   z4_kappa1                   = 0.02
   z4_kappa2                   = 0.0   z4_kappa2                   = 0.0
 +  
   punctures_lapse             = psiBL^(-2)   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
      
   #############################################################################   #############################################################################
-  # eos +  # invariants  
-  eos                         = tab4D  +   
-  eos_tab_file                = /lustre/hpe/ws10/ws10.1/ws/xuplgieg-BNS/eos/DD2_4D/DD2_4D_NEW_NEW.h5 +  ntheta                      = 47 
-  eos_tab_file_beta           = /lustre/hpe/ws10/ws10.1/ws/xuplgieg-BNS/eos/DD2_4D/DD2_4D_cold_beta_NEW.d +  nphi                        = 46  
-  eos_mb                      = 9.3046542925126903e+02 +  invariants_compute_modes    = yes 
-  eos_T     = 0.10000000000000001 +  invariants_modes_r          = 300 400 500 600 700 800 900 1000 1100 1200  
-  eos_Y     = 0.4 +  invariants_modes_lmin       = 0 
-  eos_interp     = steffen+  invariants_modes_lmax       = 5 
 +  invariants_energy_r         = 300 400 500 600 700 800 900 1000 1100 1200 
 +  gauss_codacci_mainardi      = standard 
 +  invariants_order            = 6 
 +  
   #############################################################################   #############################################################################
-  # hydro+  #ADM Mass
      
-  grhd_C2P                    = illinois_4d +  ADM_mass_ncircles     = 101 
-  grhd_C2P_NewtonRaphsonTR    = 1e-9 +  ADM_mass_npoints      = 80 
-  grhd_C2P_NewtonRaphsonNR    = 200  +  ADM_mass_lmin         = 1 
-  grhd_vmax                   = 0.99 +  ADM_mass_lmax         = 5 
-  grhd_Wlor_max               = 1e8 +  ADM_mass_r            = 300 400 500 600 700 800 900 1000 1100 1200 
-  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 +  # AHmod
-  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+  AHmod_verbose     = no 
 +  AHmod_ntheta      = 30 
 +  AHmod_nphi        = 60 
 +  AHmod_LMAX        = 10 
 +  AHmod_interpolation_order = 4 
 +  AHmod_nhorizons   = 1  
 +  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 
      
-  ################################################################################ +  ############################################################################## 
-  # M1+  # hydroanalysis
      
-grrhd_m1_implicit_sources = linear +  hydroanalysis_ejecta_spheres  = yes 
-grrhd_m1_fluxes = Radice +  hydroanalysis_ejecta_spheres_radius = 200  
- +  hydroanalysis_ejecta_nradius = 5  
-grrhd_m1_C2P_NewtonRaphsonTR = 1e-8 +  hydroanalysis_ejecta_dradius = 100 
-grrhd_m1_C2P_NewtonRaphsonNR = 100 +   
- +  hydroanalysis_Mbar_radius  = 8  
-rad_hrsc_rec              = CONST +  hydroanalysis_Mbar_nradius = 12 
-rad_hrsc_flux             = LLF +  hydroanalysis_Mbar_dradius = 1. 
- +  hydroanalysis_rATM         = 1e-25 
-grrhd_m1_decouple_fac = 1e5 +  hydroa_mode_projection     = no 
- +   
-radiation_use_mask = no +  ############################################################################ 
-optical_depth_mask = 100 +  #puncture properties 
- +   
-grrhd_m1_NuHeat = yes +  moving_puncture_finboxfix = 20  
-grrhd_m1_nls_rates = no +  moving_puncture_finboxfixv = 6  
- +  moving_puncture_properties  = yes 
-radiation_interpolate_order = 4 +  puncture_properties_r       = 8 9 10 11 12 14 15 20 
-radiation_interpolate_scheme_restriction = LINTVD +  puncture_properties_punc    = 2 
-radiation_interpolate_scheme_prolongation = LINTVD +  puncture_properties_npoints = 40 
- +  puncture_properties_circles = 41 
-rhs_radiation_null_ONlevel0 = yes +   
- +  ############################################################################# 
-radiation_at_merger     = no +  # output 
-radiation_activate_dist = 18 +   
- +  0doutput                = alpha ham momx grhd_D grhd_rho grhd_p grhd_fluid2_D grhd_fluid2_rho grhd_fluid2_p momy momz hydroa_Du  
-rad_evolution_method            = rk +  hydroa_Db hydroa_fluid2_Du hydroa_fluid2_Db 
-rad_evolution_method_rk         = rk3 +   
- +  #1doutiter               = 1 
- +  #1douttime               = 2 
-############################################################################# +  #1doutput                = alpha betax ham momx grhd_rho flagprolong rpsi4 flagrestrict grhd_fluid2_rho 
-# output +  #1doutinterpolate        = no 
- +  #1doutputall             = yes 
-stdout_flush = yes +   
- +  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 
-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 +  2dformat                = vtk binary float 
-                          nls_Qeff_nue nls_Reff_nue nls_lum_nue nls_lum_nua nls_lum_nux nls_Qeff_nua nls_Reff_nua +  2doutinterpolate        = no 
-   nls_Qeff_nux nls_Reff_nux nls_lum0_nue nls_lum0_nua nls_lum0_nux nls_Qeff_num nls_Reff_num +  2doutputall             = yes 
-   nls_Qeff_numa nls_Reff_numa nls_lum_num nls_lum_numa nls_lum0_num nls_lum0_numa +   
- +  3doutput                = grhd_rho grhd_fluid2_rho hydroa_Du hydroa_fluid2_Du alpha rpsi4 
-#1doutput = grhd_rho grhd_Y grhd_T grhd_p grhd_epsl +  3dformat                = vtk binary float 
- +  3doutinterpolate        = no 
-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 +  3doutputall = yes 
-                          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 +  #2doutputr               = sphere_data 
-   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+
  
  
Last modified: le 2024/07/23 13:48

Page Tools