[Commits] [svn:einsteintoolkit] GRHydro/branches/hot_and_MHD_temp_dev/ (Rev. 171)
cott at tapir.caltech.edu
cott at tapir.caltech.edu
Wed Nov 3 11:50:43 CDT 2010
User: cott
Date: 2010/11/03 11:50 AM
Removed:
/branches/hot_and_MHD_temp_dev/
configuration.ccl.omni, param.ccl.omni
Log:
* remove old files
File Changes:
Directory: /branches/hot_and_MHD_temp_dev/
==========================================
File [removed]: configuration.ccl.omni
Delta lines: +0 -6
===================================================================
--- branches/hot_and_MHD_temp_dev/configuration.ccl.omni 2010-11-03 16:49:49 UTC (rev 170)
+++ branches/hot_and_MHD_temp_dev/configuration.ccl.omni 2010-11-03 16:50:43 UTC (rev 171)
@@ -1,6 +0,0 @@
-# Configuration definition for thorn GRHydro
-# $Header$
-
-REQUIRES Boundary CartGrid3D SpaceMask ADMMacros
-REQUIRES THORNS: EOS_Omni
-
File [removed]: param.ccl.omni
Delta lines: +0 -528
===================================================================
--- branches/hot_and_MHD_temp_dev/param.ccl.omni 2010-11-03 16:49:49 UTC (rev 170)
+++ branches/hot_and_MHD_temp_dev/param.ccl.omni 2010-11-03 16:50:43 UTC (rev 171)
@@ -1,528 +0,0 @@
-# Parameter definitions for thorn GRHydro
-# $Header$
-
-shares:ADMBase
-
-USES int initial_shift ""
-{
-}
-
-USES KEYWORD metric_type ""
-{
-}
-
-EXTENDS KEYWORD shift_evolution_method ""
-{
- "Comoving" :: "Set the shift so the matter is approximately comoving"
-}
-
-
-#########################################
-### PARAMETRES SHARED FROM HydroBase ###
-#########################################
-
-shares: HydroBase
-
-USES CCTK_INT timelevels
-USES KEYWORD prolongation_type
-USES INT hydro_excision
-EXTENDS KEYWORD evolution_method ""
-{
- "GRHydro" :: "Use GRHydro to evolve the hydro variables"
-}
-EXTENDS KEYWORD Bvec_evolution_method ""
-{
- "GRHydro" :: "Use GRHydro to evolve the MHD variables"
-}
-
-#########################################
-### PARAMETRES SHARED FROM MoL ###
-#########################################
-
-shares: MethodOfLines
-
-USES CCTK_INT MoL_Num_Evolved_Vars
-USES CCTK_INT MoL_Num_Constrained_Vars
-USES CCTK_INT MoL_Num_SaveAndRestore_Vars
-USES CCTK_INT MoL_Max_Evolved_Array_Size
-USES CCTK_INT MoL_Num_ArrayEvolved_Vars
-
-shares:SpaceMask
-
-USES boolean use_mask
-
-shares: ADMMacros
-
-USES INT spatial_order
-
-restricted:
-
-BOOLEAN GRHydro_enable_internal_excision "Set this to 'false' to disable the thorn-internal excision." STEERABLE=ALWAYS
-{
-} "true"
-
-CCTK_INT GRHydro_hydro_excision "Turns excision automatically on in HydroBase" ACCUMULATOR-BASE=HydroBase::hydro_excision
-{
- 1:1 :: "Only '1' allowed"
-} 1
-
-CCTK_INT GRHydro_MaxNumEvolvedVars "The maximum number of evolved variables used by GRHydro" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars
-{
- 5:8 :: "dens scon[3] tau"
-} 5
-
-CCTK_INT GRHydro_MaxNumConstrainedVars "The maximum number of constrained variables used by GRHydro" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Constrained_Vars
-{
- 7:13 :: "A small range, depending on testing or not"
-} 8
-
-CCTK_INT GRHydro_MaxNumSandRVars "The maximum number of save and restore variables used by GRHydro" ACCUMULATOR-BASE=MethodofLines::MoL_Num_SaveAndRestore_Vars
-{
- 0:16 :: "A small range, depending on testing or not"
-} 16
-
-keyword recon_method "Which reconstruction method to use"
-{
- "tvd" :: "Slope limited TVD"
- "ppm" :: "PPM reconstruction"
- "eno" :: "ENO reconstruction"
-} "tvd"
-
-keyword method_type "Which type of method to use"
-{
- "RSA FV" :: "Reconstruct-Solve-Average finite volume method"
- "Flux Split FD" :: "Finite difference using Lax-Friedrichs flux splitting"
-} "RSA FV"
-
-
-boolean check_for_trivial_rp "Do check for trivial Riemann Problem"
-{
-} "yes"
-
-keyword recon_vars "Which type of variables to reconstruct"
-{
- "primitive" :: "Reconstruct the primitive variables"
- "conservative" :: "Reconstruct the conserved variables"
-} "primitive"
-
-keyword tvd_limiter "Which slope limiter to use"
-{
- "minmod" :: "Minmod"
- "vanleerMC" :: "Van Leer MC - Ian"
- "vanleerMC2" :: "Van Leer MC - Luca"
- "minmod2" :: "Minmod - Ian"
- "minmod3" :: "Minmod - Luca (or vice versa). Expect I've made a typo"
- "Superbee" :: "Superbee"
-} "minmod"
-
-real myfloor "A minimum number for the TVD reconstruction routine"
-{
- 0.0: :: "Must be positive"
-} 1.e-10
-
-boolean ppm_detect "Should the PPM solver use shock detection?"
-{
-} "no"
-
-keyword ppm_flatten "Which flattening procedure should the PPM solver use?"
-{
- "stencil_3" :: "our flattening procedure, which requires only stencil 3 and which may work"
- "stencil_4" :: "original C&W flattening procedure, which requires stencil 4"
-} "stencil_3"
-
-real ppm_epsilon "Epsilon for PPM zone flattening"
-{
- 0.0: :: "Must be positive. Default is from Colella & Woodward"
-} 0.33
-
-real ppm_omega1 "Omega1 for PPM zone flattening"
-{
- : :: "Anything goes. Default is from Colella & Woodward"
-} 0.75
-
-real ppm_omega2 "Omega2 for PPM zone flattening"
-{
- : :: "Anything goes. Default is from Colella & Woodward"
-} 10.0
-
-real ppm_omega_tracer "Omega for tracer PPM zone flattening"
-{
- : :: "Anything goes. Default is from Plewa & Mueller"
-} 0.50
-
-real ppm_epsilon_shock "Epsilon for PPM shock detection"
-{
- : :: "Anything goes. Default is from Colella & Woodward"
-} 0.01
-
-real ppm_eta1 "Eta1 for PPM shock detection"
-{
- : :: "Anything goes. Default is from Colella & Woodward"
-} 20.0
-
-real ppm_eta2 "Eta2 for PPM shock detection"
-{
- : :: "Anything goes. Default is from Colella & Woodward"
-} 0.05
-
-real ppm_k0 "K0 for PPM shock detection"
-{
- : :: "Anything goes. Default suggested by Colella & Woodward is: (polytropic constant)/10.0"
-} 0.2
-
-real ppm_small "A floor used by PPM shock detection"
-{
- 0.0:1.0 :: "In [0,1]"
-} 1.e-7
-
-int ppm_mppm "Use modified (enhanced) ppm scheme"
-{
- 0:1 :: "0 (off, default) or 1 (on)"
-} 0
-
-int ppm_mppm_debug_eigenvalues "write eigenvalues into debug grid variables"
-{
- 0:1 :: "0 (off, default) or 1 (on)"
-} 0
-
-int eno_order "The order of accuracy of the ENO reconstruction"
-{
- 1:* :: "Anything above 1, but above 5 is pointless"
-} 2
-
-keyword riemann_solver "Which Riemann solver to use"
-{
- "Roe" :: "Standard Roe solver"
- "Marquina" :: "Marquina flux"
- "HLLE" :: "HLLE"
-} "HLLE"
-
-
-keyword HLLE_type "Which HLLE type to use"
-{
- "Standard" :: "Standard HLLE solver"
- "Tadmor" :: "Tadmor's simplification of HLLE"
-} "Standard"
-
-
-
-keyword left_eigenvectors "How to compute the left eigenvectors"
-{
- "analytical" :: "Analytical left eigenvectors"
- "numerical" :: "Numerical left eigenvectors"
-} "analytical"
-
-keyword numerical_viscosity "How to compute the numerical viscosity"
-{
- "fast" :: "Fast numerical viscosity"
- "normal" :: "Normal numerical viscosity"
-} "fast"
-
-keyword bound "Which boundary condition to use - FIXME"
-{
- "flat" :: "Zero order extrapolation"
- "none" :: "None"
- "static" :: "Static, no longer supported"
- "scalar" :: "Constant"
-} "none"
-
-int GRHydro_stencil "Width of the stencil"
-{
- 0: :: "Must be positive"
-} 2
-
-# not used for EOS_Omni
-keyword GRHydro_eos_type "Type of Equation of State"
-{
- "Polytype" :: "P = P(rho) or P = P(eps)"
- "General" :: "P = P(rho, eps)"
-} "General"
-
-# Used with EOS_Omni
-string GRHydro_eos_table "Name for the Equation of State"
-{
- .* :: "Can be anything"
-} "Ideal_Fluid"
-
-real GRHydro_eos_rf_prec "Precision to which root finding should be carried out"
-{
- (0.0:* :: "anything larger than 0 goes"
-} 1.0e-8
-
-
-####################### Other Parameters ##############################
-
-#Parameters controlling conservative <-> primitive change.
-
-int GRHydro_c2p_warnlevel "Warnlevel for Con2Prim warnings" STEERABLE=ALWAYS
-{
- 0:1 :: "Either 0 or 1"
-} 0
-
-# This parameter is not used anymore in maintained code: the warnings are produced only for really evolved points.
-int GRHydro_c2p_warn_from_reflevel "Start warning on given refinement level and on higher levels" STEERABLE=ALWAYS
-{
- 0: :: "Greater or equal to 0"
-} 0
-
-boolean c2p_reset_pressure "If the pressure guess is unphysical should we recompute it?"
-{
-} "no"
-
-real c2p_reset_pressure_to_value "The value to which the pressure is reset to when a failure occurrs in C2P"
-{
- 0: :: "greater than zero"
-} 1.e-20
-
-real GRHydro_eps_min "Minimum value of internal energy - this is now only used in GRHydro_InitData's GRHydro_Only_Atmo routine"
-{
- 0: :: "Positive"
-} 1.e-10
-
-real GRHydro_perc_ptol "Tolerance for primitive variable solve (percent)"
-{
- 0: :: "Do we really want both tolerances?"
-} 1.e-8
-
-real GRHydro_del_ptol "Tolerance for primitive variable solve (absolute)"
-{
- 0: :: "Do we really want both tolerances?"
-} 1.e-18
-
-int GRHydro_countmax "Maximum number of iterations for Con2Prim solve"
-{
- 1:* :: "Must be positive"
-} 100
-
-int GRHydro_countmin "Minimum number of iterations for Con2Prim solve"
-{
- 0:* :: "Must be non negative"
-} 1
-
-int GRHydro_polish "Number of extra iterations after root found"
-{
- 0:* :: "Must be non negative"
-} 0
-
-###Parameters for compact fluid objects
-
-REAL GRHydro_rho_central "Central Density for Star"
-{
- : :: ""
-} 1.e-5
-
-##Parameters for the atmosphere.
-
-REAL rho_abs_min "A minimum rho below which evolution is turned off (atmosphere). If negative, the relative minimum will be used instead."
-{
- -1.0: :: ""
-} -1.0
-
-REAL rho_rel_min "A minimum relative rho (rhomin = centden * rho_rel_min) below which evolution is turned off (atmosphere). Only used if rho_abs_min < 0.0"
-{
- 0: :: ""
-} 1.e-9
-
-REAL initial_rho_abs_min "An absolute value for rho in the atmosphere. To be used by initial data routines only. Unused if negative."
-{
- -1.0: :: ""
-} -1.0
-
-REAL initial_rho_rel_min "A relative (to the central density) value for rho in the atmosphere. To be used by initial data routines only. Unused if negative."
-{
- -1.0: :: ""
-} -1.0
-
-REAL initial_atmosphere_factor "A relative (to the initial atmosphere) value for rho in the atmosphere. To be used by initial data routines only. Unused if negative."
-{
- -1.0: :: ""
-} -1.0
-
-REAL rho_abs_min_after_recovery "A new absolute value for rho in the atmosphere. To be used after recovering. Unused if negative." STEERABLE=ALWAYS
-{
- -1.0: :: ""
-} -1.0
-
-REAL GRHydro_atmo_tolerance "A point is set to atmosphere in the Con2Prim's if its rho < GRHydro_rho_min *(1+GRHydro_atmo_tolerance). This avoids occasional spurious oscillations in carpet buffer zones lying in the atmosphere (because prolongation happens on conserved variables)" STEERABLE=ALWAYS
-{
- 0.0: :: "Zero or larger. A useful value could be 0.0001"
-} 0.0
-
-BOOLEAN wk_atmosphere "Use some of Wolfgang Kastauns atmosphere tricks"
-{
-} "no"
-
-BOOLEAN Check_Rho_Minimum "Should a check on rho < GRHydro_rho_min be performed and written as WARNING level 2?"
-{
-} "no"
-
-boolean EoS_Change "Recalculate fluid quantities if changing the EoS"
-{
-} "no"
-
-keyword EoS_Change_type "Change polytropic K or Gamma?"
-{
- "K" :: "Change the K"
- "Gamma" :: "Change the Gamma"
- "GammaKS" :: "Change K and Gamma, Shibata et al. 2004 3-D GR Core Collapse style"
-} "Gamma"
-
-real initial_Gamma "If changing Gamma, what was the value used in the initial data routine?"
-{
- (0.0: :: "Positive"
-} 1.3333
-
-real initial_k "If changing K, what was the value used in the initial data routine?"
-{
- (0.0: :: "Positive"
-} 100.0
-
-#If you want to use FishEye, you need the weighted fluxes
-
-BOOLEAN use_weighted_fluxes "Weight the flux terms by the cell surface areas"
-{
-} "no"
-
-#For the comoving shift
-
-real comoving_factor "Factor multiplying the velocity for the comoving shift"
-{
- 0.0:2.0 :: "[0,2] is allowed, but [0,1] is probably reasonable"
-} 0.0
-
-keyword comoving_attenuate "Which attenuation method for the comoving shift"
-{
- "sqrt" :: "Multiply by sqrt(rho/rho_max)"
- "tanh" :: "Multiply by 1/2*(1+tanh(factor*(rho/rho_max - offset)))"
-} "tanh"
-
-keyword comoving_v_method "Which method for getting the radial velocity"
-{
- "projected" :: "vr = x_i . v^i / r"
- "components" :: "vr = sqrt(v_i v^i)"
-} "projected"
-
-real comoving_tanh_factor "The factor in the above tanh attenuation"
-{
- (0.0:* :: "Any positive number"
-} 10.0
-
-real comoving_tanh_offset "The offset in the above tanh attenuation"
-{
- 0.0:1.0 :: "Only makes sense in [0,1]"
-} 0.05
-
-int set_trivial_rp_grid_function "set gf for triv. rp (only for debugging)"
-{
- 0:1 :: "0 for no (default), 1 for yes"
-} 0
-
-#Advect a tracer?
-
-boolean evolve_tracer "Should we advect tracers?"
-{
-} "no"
-
-int number_of_tracers "Number of tracer variables to be advected"
-{
- 0:* :: "positive or zero"
-} 0
-
-boolean use_min_tracer "Should there be a floor on the tracer?"
-{
-} "no"
-
-real min_tracer "The floor placed on the tracer"
-{
- *:* :: "Anything"
-} 0.0
-
-#Track particles?
-int number_of_particles "Number of particles to track" ACCUMULATOR-BASE=MethodOfLines::MoL_Max_Evolved_Array_Size
-{
- 0:* :: "0 switches off particle tracking"
-} 0
-
-int number_of_arrays "Number of arrays to evolve" ACCUMULATOR-BASE=MethodOfLines::MoL_Num_ArrayEvolved_Vars
-{
- 0:3 :: "Either zero or three, depending on the particles"
-} 0
-
-STRING particle_interpolator "What interpolator should be used for the particles"
-{
- ".+" :: "A valid interpolator name"
-} "Lagrange polynomial interpolation"
-
-CCTK_INT particle_interpolation_order "What order should be used for the particle interpolation"
-{
- 1:* :: "A valid positive interpolation order"
-} 2
-
-#Finding the local gradient
-
-keyword gradient_method "Which method is used to set GRHydro::DiffRho?"
-{
- "First diff" :: "Standard first differences"
- "Curvature" :: "Curvature based method of Paramesh / FLASH"
- "Rho weighted" :: "Based on the size of rho"
-} "First diff"
-
-#For the new EOS interface
-
-BOOLEAN use_eosgeneral "Should we use the new EOS interface?"
-{
-} "no"
-
-STRING eosgeneral_name "The name of the EOS to be used"
-{
- ".*" :: "Any registered EOS"
-} "Polytrope"
-
-STRING eosgeneral_indep_names "The names of additional variables to be used in the EOS call"
-{
- ".*" :: "Any variables consistent with the EOS"
-} ""
-
-INT eosgeneral_n_indeps "The number of additional variables to be used in the EOS call"
-{
- 0:* :: "Anything non-negative"
-} 0
-
-STRING eosgeneral_indep_gfs "The names of the grid functions holding the additional variables"
-{
- ".*" :: "Should be in the full Thorn::GF format"
-} ""
-
-STRING eosgeneral_indep_plus_gfs "The names of the grid functions holding the additional variables that will be correct at the plus cell boundary"
-{
- ".*" :: "Should be in the full Thorn::GF format"
-} ""
-
-STRING eosgeneral_indep_minus_gfs "The names of the grid functions holding the additional variables that will be correct at the minus cell boundary"
-{
- ".*" :: "Should be in the full Thorn::GF format"
-} ""
-
-#Setting outflow boundaries over only some of the domain for multipatch runs
-
-BOOLEAN outflow_boundaries "Use the outflow boundary conditions?"
-{
-} "no"
-
-REAL outflowboundary_normal[3] "The normal to the plane; to the right of this outflow conditions will be set"
-{
- *:* :: "Anything"
-} 0.0
-
-# NaN detection level
-
-INT GRHydro_NaN_verbose "The warning level for NaNs occuring within GRHydro"
-{
- 0:* :: "The warning level"
-} 2
-
-CCTK_REAL GRHydro_lorentz_overshoot_cutoff "Set the Lorentz factor to this value in case it overshoots (1/0)"
-{
- 0:* :: "Some big value"
-} 1.e100
-
More information about the Commits
mailing list