This is a _temporary_ repository to be able to start to work on the

code right now. I have put in the public version of Whisky to start from.

Everybody with commit rights should get commit messages (and the other
way around). It should not be a problem to add people to that list, just
ask. I don't want to get into political problems because someone feels
excluded, but I also don't want to give everyone access per se.

Frank



git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@3 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

3 files changed, 254 insertions, 0 deletions

diff --git a/par/BenchRNS_Carpet_3level.par b/par/BenchRNS_Carpet_3level.par
new file mode 100644
index 0000000..ad54a30
--- /dev/null
+++ b/par/BenchRNS_Carpet_3level.par
@@ -0,0 +1,196 @@
+ 
+#--- Base and Driver Thorns
+
+ActiveThorns	=  "	time 
+			carpet 
+			carpetlib 
+			carpetregrid 
+			carpetreduce 
+			cartgrid3d 
+			coordbase 
+			mol 
+			boundary 
+			spacemask 
+			carpetiobasic 
+			ioutil 
+			symbase 
+			aeilocalinterp 
+			carpetinterp 
+		   "
+
+#--- Spacetime :-)
+
+ActiveThorns	=  "	admbase 
+			staticconformal 
+			admconstraints 
+			admcoupling 
+			coordgauge 
+			admmacros 
+			bssn_mol 
+		   "
+
+#--- Hydrodynamics
+
+ActiveThorns	=  "	whisky 
+			whisky_init_data 
+			whisky_rnsid 
+			hydrobase
+			eos_base 
+			eos_polytrope 
+		   "
+
+
+################## Output Related Stuff ##################################
+
+IOBasic::outInfo_vars			=	"HydroBase::rho admbase::alp"
+iobasic::outInfo_every			= 	1
+
+
+#--- general
+
+
+#--- output directories
+IO::out_dir					=	$parfile
+
+#--- checkpoint recovery
+io::recover 					= 	"no"
+
+
+#--- output frequency
+carpetiobasic::real_min                         =       1.0e-2
+io::out_single_precision 			= 	"yes"
+
+
+
+############################# DRIVER  #####################################
+
+#--- Benchmarking
+
+carpet::constant_load_per_processor		= "yes"
+
+#--- performance & infrastructure
+
+carpet::enable_all_storage 			= "no"
+CarpetLib::save_memory_during_regridding 	= "yes"
+
+#--- Carpet prolongation order settings
+
+Carpet::prolongation_order_space		= 3
+Carpet::prolongation_order_time			= 2
+
+
+#--- Carpet initial data setup parameters
+
+Carpet::init_3_timelevels 		= "no"
+Carpet::init_each_timelevel 		= "no"
+
+MoL::initial_data_is_crap		= "yes"
+
+
+#--- Grid Setup and carpet parameters
+
+
+carpet::buffer_width 				= 6
+
+driver::ghost_size                      	= 3
+
+cartgrid3d::type                        	= "box"
+cartgrid3d::domain                      	= "full"
+cartgrid3d::avoid_origin			= "no"
+
+cartgrid3d::dxyz				= 1.0
+driver::global_nx                  		= 30
+driver::global_ny                  		= 30
+driver::global_nz                  		= 30
+
+cactus::cctk_itlast                     	= 128
+
+Carpet::domain_from_coordbase 			= "no"
+
+carpet::max_refinement_levels           	= 3
+carpetregrid::refinement_levels         	= 3
+
+carpet::veryverbose 				= "no"
+carpet::verbose					= "no"
+
+
+################## Spacetime Related Stuff #############################
+
+
+admbase::evolution_method			=	"adm_bssn"
+admbase::lapse_evolution_method			=	"1+log"
+admbase::shift_evolution_method 		= 	"static"
+
+boundary::radpower				=	2
+
+admconstraints::bound				=	"static"
+
+bssn_mol::timelevels				=	3
+
+
+######################## Hydro Related Stuff ###########################
+
+
+HydroBase::timelevels				=	3
+
+whisky::whisky_rho_central			=	1.61930347e-08
+whisky::riemann_solver				=	"Marquina"
+whisky::whisky_eos_type				=	"Polytype"
+whisky::whisky_eos_table			=	"2D_Polytrope"
+whisky::recon_method	            		=       "ppm"
+whisky::ppm_detect				=	"yes"
+whisky::whisky_stencil          		=       3
+whisky::bound		                	=       "none"
+
+whisky::EoS_Change				=	"no"
+EOS_Polytrope::eos_gamma                	=       2.0
+EOS_Polytrope::eos_k                    	=       1.0
+EOS_Polytrope::use_cgs				=	"yes"
+EOS_Polytrope::gamma_ini			=	2.0
+
+#------ Atmosphere treatment
+
+SpaceMask::use_mask 				= 	yes
+
+whisky::rho_rel_min				=	1.e-5
+
+#------ Atmosphere only for
+whisky::initial_atmosphere_factor		=	0.99
+
+
+
+####################### Initial Data #####################################
+
+
+admbase::initial_data				=	"rnsid"
+admbase::initial_lapse				=	"rnsid"
+admbase::initial_shift				=	"rnsid"
+admbase::metric_type				=	"physical"
+
+rnsid::axes_ratio				=	1.0
+rnsid::accuracy					=	1e-7
+
+rnsid::RNS_Gamma				=	2.0
+rnsid::RNS_K					=	100.0	
+
+rnsid::recover_2Dmodel				=	"no"
+rnsid::save_2Dmodel				=	"no"
+
+rnsid::mdiv					=	151
+rnsid::sdiv					=	301
+
+
+
+
+##################### Infrastructure Related Stuff ########################
+
+#--- Time integration
+
+MoL::ODE_Method					=	"icn"
+MoL::MoL_Intermediate_Steps 			= 	3
+time::dtfac                             	= 	0.375
+
+
+
+#---- Timer output
+cactus::cctk_timer_output               	=       "FULL"
diff --git a/par/Whisky.th b/par/Whisky.th
new file mode 100644
index 0000000..30d60e6
--- /dev/null
+++ b/par/Whisky.th
@@ -0,0 +1,19 @@
+CactusMoL/MoL
+CactusBase/Boundary        
+CactusBase/CartGrid3D      
+CactusBase/IOUtil          
+CactusBase/IOBasic         
+CactusBase/Time            
+CactusBase/IOASCII       
+CactusPUGH/PUGH            
+CactusPUGH/PUGHSlab
+CactusPUGH/PUGHReduce
+EUHydro/Einstein
+HydroBase/HydroBase
+Whisky/Whisky
+Whisky/Whisky_Init_Data
+CactusUtils/NaNChecker
+CactusEOS/EOS_Base
+CactusEOS/EOS_Ideal_Fluid
+CactusEOS/EOS_1D_Base
+CactusEOS/EOS_1D_Ideal_Polytrope
diff --git a/par/whisky.par b/par/whisky.par
new file mode 100644
index 0000000..33f8f12
--- /dev/null
+++ b/par/whisky.par
@@ -0,0 +1,39 @@
+ActiveThorns = "time MoL pugh pughreduce pughslab cartgrid3d ioutil ioascii iobasic boundary einstein hydrobase whisky whisky::whisky_init_data nanchecker eos_base eos_ideal_fluid"
+
+#This is a hack. ghost_size should be the same as whisky::whisky_stencil.
+
+pugh::ghost_size=2
+
+time::dtfac = 0.25
+
+methodoflines::ode_method = "rk2"
+
+einstein::evolution_system="none"
+einstein::shift="none"
+einstein::use_conformal = "no"
+einstein::use_conformal_derivs = "no"
+einstein::initial_data="shocktube"
+whisky_init_data::shocktube_type="diagshock"
+whisky_init_data::shock_case="Simple"
+
+EOS_Ideal_Fluid::eos_ideal_fluid_gamma=1.66666666666666
+whisky::riemann_solver="Roe"
+whisky::tvd_limiter="minmod"
+whisky::bound="flat"
+whisky::whisky_eos_table="Ideal_Fluid"
+
+grid::type = "BySpacing"
+grid::domain = "full"
+grid::dxyz = 0.5
+
+driver::global_nx = 15
+driver::global_ny = 15
+driver::global_nz = 15
+
+cactus::cctk_itlast = 5
+
+IO::outdir = "whisky_shocktube"
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_vars = "HydroBase::rho"
+IOASCII::out1D_every = 1
+IOASCII::out1D_vars = "HydroBase::rhoHydroBase::press HydroBase::eps  HydroBase::vel whisky::dens whisky::tau whisky::scon"