Update params

Examples/Fluid_Kerr/params.txt

@@ -8,21 +8,21 @@ verbosity = 0
 
 # location / naming of output files
 # output_path = "" # Main path for all files. Must exist!
-chk_prefix = PerfectFluid_
-plot_prefix = PerfectFluidp_
-# restart_file = PerfectFluid_000000.3d.hdf5
+chk_prefix = FluidKerrBH_
+plot_prefix = FluidKerrBHp_
+#restart_file = 
 
 # HDF5files are written every dt = L/N*dt_multiplier*checkpoint_interval
-checkpoint_interval = 500
+checkpoint_interval = 1
 # set to 0 to turn off plot files (except at t=0 and t=stop_time)
 # set to -1 to never ever print plotfiles
-plot_interval = 1
-num_plot_vars = 2
-plot_vars = rho eps
+plot_interval = 0
+num_plot_vars = 0
+# plot_vars = 
 
 # subpaths - specific directories for hdf5, pout, extraction data
 # (these are created at runtime)
-hdf5_subpath = "hdf5"
+hdf5_subpath = "/sakura/ptmp/dtraykova/fluid_test_Kerr_4lvl/hdf5"
 pout_subpath = "pout"
 data_subpath = "data"
 
@@ -41,16 +41,17 @@ print_progress_only_to_rank_0 = 1
 # Here we decouple the evolution so the scalar evolved on the
 # metric background without backreaction (this avoids the need
 # to solve the constaints)
-G_Newton = 0.0
-
-# Scalar field initial data
-#scalar_amplitude = 0.1
-#scalar_width = 5.0
-#scalar_mass = 0.2
+G_Newton = 1.0
 
+# Fluid params
+fluid_rho0 = 1.0
+fluid_width = 5.
+fluid_delta = 0.2
+
 # Kerr BH data
-kerr_mass = 0.
-kerr_spin = 0.0
+kerr_mass = 1.0
+kerr_spin = 0.5
+# kerr_center = 64 64 64 # defauts to center
 
 #################################################
 # Grid parameters
@@ -62,28 +63,25 @@ kerr_spin = 0.0
 # NB - if you have a non-cubic grid, you can specify 'N1' or 'N1_full',
 # 'N2' or 'N2_full' and 'N3' or 'N3_full' ( then dx_coarsest = L/N(max) )
 # NB - the N values need to be multiples of the block_factor
-N_full = 64 #128
-L_full = 128 #256
+N_full = 32
+L_full = 128
 
 # Maximum number of times you can regrid above coarsest level
-max_level = 5 # There are (max_level+1) grids, so min is zero
+max_level = 4 # There are (max_level+1) grids, so min is zero
 
 # Frequency of regridding at each level and thresholds on the tagging
 # Need one for each level except the top one, ie max_level items
 # Generally you do not need to regrid frequently on every level
-# in this example turn off regridding on all levels
-# Level Regridding: 0   1   2   3   4   5
-regrid_interval   = 0   0   0   0   0   0
-# regrid_threshold  = 0.5
-
-regrid_interval = 0 0 0 0 0 0
+# Level Regridding:   0   1   2   3   4
+#regrid_interval   =  50  25   5   5   5
+#regrid_thresholds = 0.3 0.3 0.3 0.3 0.3
+regrid_interval   = 0 0 0 0
 
 # Max and min box sizes
 max_box_size = 16
 min_box_size = 16
 
-tag_buffer_size = 0 # this example uses a fixed grid
-
+# tag_buffer_size = 3
 # grid_buffer_size = 8
 # fill_ratio = 0.7
 # num_ghosts = 3
@@ -92,14 +90,13 @@ tag_buffer_size = 0 # this example uses a fixed grid
 #################################################
 # Boundary Conditions parameters
 
-#Periodic directions - 0 = false, 1 = true
-isPeriodic = 1 1 1
+# Periodic directions - 0 = false, 1 = true
+isPeriodic = 0 0 0
 # if not periodic, then specify the boundary type
 # 0 = static, 1 = sommerfeld, 2 = reflective
-# 3 = extrapolating, 4 = mixed
 # (see BoundaryConditions.hpp for details)
-hi_boundary = 4 4 4
-lo_boundary = 2 2 2
+hi_boundary = 3 3 3
+lo_boundary = 3 3 2
 
 # if reflective boundaries selected, must set
 # parity of all vars (in order given by UserVariables.hpp)
@@ -111,8 +108,8 @@ vars_parity            = 0 0 4 6 0 5 0    #chi and hij
                          0 0 4 6 0 5 0    #K and Aij
                          0 1 2 3          #Theta and Gamma
                          0 1 2 3 1 2 3    #lapse shift and B
-                         0 1 2 3 0	  #D Sj and tau
-			 1 2 3 0 0        #vi rho and epsilon
+                         0 1 2 3 0        #D Sj and tau
+                         1 2 3 0 0        #vi rho and epsilon
 vars_parity_diagnostic = 0 1 2 3          #Ham and Mom
 
 # if sommerfeld boundaries selected, must select
@@ -122,17 +119,17 @@ nonzero_asymptotic_vars = chi h11 h22 h33 lapse
 nonzero_asymptotic_values = 1.0 1.0 1.0 1.0 1.0
 
 # if you are using extrapolating BC:
-extrapolation_order = 0 #1
-num_extrapolating_vars = 2
-extrapolating_vars = phi Pi
+extrapolation_order = 0
+num_extrapolating_vars = 10
+extrapolating_vars = D Sj1 Sj2 Sj3 tau vi1 vi2 vi3 rho eps
 
 #################################################
 # Evolution parameters
 
 # dt will be dx*dt_multiplier on each grid level
-dt_multiplier = 0.078125 #0.25
-stop_time = 5.0
-# max_steps = 4
+dt_multiplier = 0.1
+stop_time = 100.0
+#max_steps = 4
 
 # Spatial derivative order (only affects CCZ4 RHS)
 max_spatial_derivative_order = 4 # can be 4 or 6

Examples/PerfectFluid_Kerr/test_Gauss_small/params.txt

@@ -0,0 +1,162 @@
+# See the wiki page for an explanation of the params!
+# https://github.com/GRChombo/GRChombo/wiki/Guide-to-parameters
+
+#################################################
+# Filesystem parameters
+
+verbosity = 0
+
+# location / naming of output files
+# output_path = "" # Main path for all files. Must exist!
+chk_prefix = FluidKerrBH_
+plot_prefix = FluidKerrBHp_
+#restart_file = /sakura/ptmp/dtraykova/Kerr/hdf5/KerrBH_000060.3d.hdf5
+
+# HDF5files are written every dt = L/N*dt_multiplier*checkpoint_interval
+checkpoint_interval = 1
+# set to 0 to turn off plot files (except at t=0 and t=stop_time)
+# set to -1 to never ever print plotfiles
+plot_interval = 0
+num_plot_vars = 0
+# plot_vars = 
+
+# subpaths - specific directories for hdf5, pout, extraction data
+# (these are created at runtime)
+hdf5_subpath = "/sakura/ptmp/dtraykova/fluid_test_Kerr_4lvl/hdf5"
+pout_subpath = "pout"
+data_subpath = "data"
+
+# change the name of output files
+# pout_prefix = "pout"
+print_progress_only_to_rank_0 = 1
+
+# ignore_checkpoint_name_mismatch = 0
+# write_plot_ghosts = 0
+
+#################################################
+# Initial Data parameters
+
+# Change the gravitational constant of the Universe!
+# Default is 1.0, for standard geometric units
+# Here we decouple the evolution so the scalar evolved on the
+# metric background without backreaction (this avoids the need
+# to solve the constaints)
+G_Newton = 1.0
+
+# Fluid params
+fluid_rho0 = 1.0
+fluid_width = 5.
+fluid_delta = 0.2
+
+# Kerr BH data
+kerr_mass = 1.0
+kerr_spin = 0.5
+# kerr_center = 64 64 64 # defauts to center
+
+#################################################
+# Grid parameters
+
+# 'N' is the number of subdivisions in each direction of a cubic box
+# 'L' is the length of the longest side of the box, dx_coarsest = L/N
+# NB - If you use reflective BC and want to specify the subdivisions and side
+# of the box were there are no symmetries, specify 'N_full' and 'L_full' instead
+# NB - if you have a non-cubic grid, you can specify 'N1' or 'N1_full',
+# 'N2' or 'N2_full' and 'N3' or 'N3_full' ( then dx_coarsest = L/N(max) )
+# NB - the N values need to be multiples of the block_factor
+N_full = 32
+L_full = 128
+
+# Maximum number of times you can regrid above coarsest level
+max_level = 4 # There are (max_level+1) grids, so min is zero
+
+# Frequency of regridding at each level and thresholds on the tagging
+# Need one for each level except the top one, ie max_level items
+# Generally you do not need to regrid frequently on every level
+# Level Regridding:   0   1   2   3   4
+#regrid_interval   =  50  25   5   5   5
+#regrid_thresholds = 0.3 0.3 0.3 0.3 0.3
+regrid_interval   = 0 0 0 0
+
+# Max and min box sizes
+max_box_size = 16
+min_box_size = 16
+
+# tag_buffer_size = 3
+# grid_buffer_size = 8
+# fill_ratio = 0.7
+# num_ghosts = 3
+# center = 256.0 256.0 256.0 # defaults to center of the grid
+
+#################################################
+# Boundary Conditions parameters
+
+# Periodic directions - 0 = false, 1 = true
+isPeriodic = 0 0 0
+# if not periodic, then specify the boundary type
+# 0 = static, 1 = sommerfeld, 2 = reflective
+# (see BoundaryConditions.hpp for details)
+hi_boundary = 3 3 3
+lo_boundary = 3 3 2
+
+# if reflective boundaries selected, must set
+# parity of all vars (in order given by UserVariables.hpp)
+# 0 = even
+# 1,2,3 = odd x, y, z
+# 4,5,6 = odd xy, yz, xz
+# 7     = odd xyz
+vars_parity            = 0 0 4 6 0 5 0    #chi and hij
+                         0 0 4 6 0 5 0    #K and Aij
+                         0 1 2 3          #Theta and Gamma
+                         0 1 2 3 1 2 3    #lapse shift and B
+                         0 1 2 3 0        #D Sj and tau
+                         1 2 3 0 0        #vi rho and epsilon
+vars_parity_diagnostic = 0 1 2 3          #Ham and Mom
+
+# if sommerfeld boundaries selected, must select
+# non zero asymptotic values
+num_nonzero_asymptotic_vars = 5
+nonzero_asymptotic_vars = chi h11 h22 h33 lapse
+nonzero_asymptotic_values = 1.0 1.0 1.0 1.0 1.0
+
+# if you are using extrapolating BC:
+extrapolation_order = 0
+num_extrapolating_vars = 10
+extrapolating_vars = D Sj1 Sj2 Sj3 tau vi1 vi2 vi3 rho eps
+
+#################################################
+# Evolution parameters
+
+# dt will be dx*dt_multiplier on each grid level
+dt_multiplier = 0.1
+stop_time = 100.0
+#max_steps = 4
+
+# Spatial derivative order (only affects CCZ4 RHS)
+max_spatial_derivative_order = 4 # can be 4 or 6
+
+nan_check = 1
+
+# Lapse evolution
+lapse_advec_coeff = 1.0
+lapse_coeff = 2.0
+lapse_power = 1.0
+
+# Shift evolution
+shift_advec_coeff = 0.0 # Usually no advection for beta
+shift_Gamma_coeff = 0.75
+eta = 1.0 # eta of gamma driver, should be of order ~1/M_ADM of spacetime
+
+# CCZ4 parameters
+formulation = 1 # 1 for BSSN, 0 for CCZ4
+kappa1 = 0.
+kappa2 = 0.
+kappa3 = 0.
+covariantZ4 = 1 # 0: keep kappa1; 1 [default]: replace kappa1 -> kappa1/lapse
+
+# coefficient for KO numerical dissipation
+sigma = 0.3
+
+# min_chi = 1.e-4
+# min_lapse = 1.e-4
+
+#################################################