Abstract
GRHydro_Init_Data - some initial data for GRHydro
This thorn generates some initial data for the GRHydro code. There are more (and more physically interesting) initial-data codes in other thorns. As with the GRHydro code itself, please feel free to add, alter or extend any part of this code. However please keep the documentation up to date (even, or especially, if it’s just to say what doesn’t work).
Currently this thorn contains a few tests that should really be test suites, some shock-tube initial data, a routine to set atmopshere everywhere on the grid (for tests), a routine to read initial data from files (not very generic, tough) and a routine to set up the simple-wave initial data .
There are tests of the TVD reconstruction routine and of the routines that convert between conservative and primitive variables. These all just produce output to the screen or to fort.* files. The reconstruction test outputs the function to be reconstructed and the boundary-extended values. The conservative-to-primitive test just outputs the two sets of variables. If you haven’t altered the code an if you set
eos_polytrope::eos_gamma = 2.0 eos_polytrope::eos_k = 100.0
(which are the defaults), the output should be
primitive variables: rho : 1.29047172182043 velx : 9.902578465178671E-004 vely : 9.902578465178671E-004 velz : 9.902578465178671E-004 eps : 0.374770481293314 press : 166.531726481819 w_lor : 1.00000147091915
The conservative to primitive to conservative test outputs the initial and final data which should agree.
There are three possible shock-tube problems, referred to as Sod, Simple and Blast, with initial data
Type | \(\rho _{_L}\) | \(v^i_{_L}\) | \(\varepsilon _{_L}\) | \(\rho _{_R}\) | \(v^i_{_R}\) | \(\varepsilon _{_R}\) |
Sod | 1 | 0 | 1.5 | 0.125 | 0 | 0.15 |
Simple | 10 | 0 | 20 | 1 | 0 | \(10^{-6}\) |
Blast | 1 | 0 | 1500 | 1 | 0 | \(1.5\cdot 10^{-2}\) |
The shock shape can be planar (along each axis or along the main diagonal) or spherical and the position of the plane or of the center of the sphere can be chosen though parameters. If a diagonal shock is selected, the initial data is set to either the left or right state depending on where the centre of the cell falls. Cleverer routines that weight the initial data to avoid “staircasing” may be added if there is demand. For more discussion on shock tubes see [1].
For testing purposes, this routine sets all the points to the values of the atmosphere.
This routine testes initial data for a simple wave with sinusoidal initial function for the velocity, as described in Anile, Miller, Motta, Formation and damping of relativistic strong shocks,Phys. Fluids 26, 1450 (1983).
[1] J. M. Martí and E. Müller. Numerical hydrodynamics in Special Relativity. Living Rev. Relativity, 3, 1999. [Article in online journal], cited on 31/7/01, http://www.livingreviews.org/Articles/Volume2/1999-3marti/index.html.
[2] J. A. Font, M. Miller, W. Suen and M. Tobias. Three Dimensional Numerical General Relativistic Hydrodynamics I: Formulations, Methods, and Code Tests Phys. Rev., D61, 044011, 2000.
advectedloop_case | Scope: private | KEYWORD |
Description: Vẑ=0 or not?
| ||
Range | Default: Vẑ=0 | |
Vẑ=0 | Useful to evaluate divB deviations
| |
Vẑ/=0 | Useful to evaluate con2prim robustness in keeping Vẑ const.
| |
advectedloop_dela | Scope: private | KEYWORD |
Description: How to calculate Bî field from the potential A
| ||
Range | Default: Exact | |
Exact | Analytic, exact closed formula applied
| |
Numeric | Finite difference approximation of the derivatives applied
| |
advectedloop_type | Scope: private | KEYWORD |
Description: 2-dimensional or 3-dimensional?
| ||
Range | Default: 2D | |
2D | 2-dimensional (Bẑ=0)
| |
3D | 3-dimensional (B=0, where B —— oblique cylinder axis.
| |
alfvenwave_pressure | Scope: private | REAL |
Description: P_gas for the Alfven wave
| ||
Range | Default: 1.0 | |
(0:* | positive
| |
alfvenwave_type | Scope: private | KEYWORD |
Description: 1-dimensional or 2-dimensional?
| ||
Range | Default: 1D | |
1D | 1-dimensional
| |
2D | 2-dimensional (in x-y plane)
| |
atmosphere_vel | Scope: private | REAL |
Description: Velocity of the atmosphere if non-trivial
| ||
Range | Default: 0.0 | |
*:* | Anything
| |
attenuate_atmosphere | Scope: private | BOOLEAN |
Description: Attenuate the velocity in the atmosphere
| ||
Default: no | ||
bh_bondi_pos_x | Scope: private | REAL |
Description: X-coordinate of black hole in Bondi solution
| ||
Range | Default: 0.0 | |
*:* | anything
| |
bh_bondi_pos_y | Scope: private | REAL |
Description: Y-coordinate of black hole in Bondi solution
| ||
Range | Default: 0.0 | |
*:* | anything
| |
bh_bondi_pos_z | Scope: private | REAL |
Description: Z-coordinate of black hole in Bondi solution
| ||
Range | Default: 0.0 | |
*:* | anything
| |
bondi_beta_sonicpt | Scope: private | REAL |
Description: Plasma beta parameter at the sonic point. Calculate bondi_bmag afterwards.
| ||
Range | Default: 1.0 | |
(0:* | positive
| |
bondi_bmag | Scope: private | REAL |
Description: B_0 parameter for magnetized Bondi
| ||
Range | Default: 0.01 | |
0:* | Anything positive
| |
bondi_bvec_method | Scope: private | KEYWORD |
Description: how to compute the magnetic field vector
| ||
Range | Default: direct | |
direct | directly from Cartesian metric
| |
transform | transform Schwarzschild solution to Kerr Schild | |
bondi_central_mass | Scope: private | REAL |
Description: Mass of central object to find Bondi solution about
| ||
Range | Default: 1.0 | |
(0:* | positive
| |
bondi_central_spin | Scope: private | REAL |
Description: Dimensionless spin of central object within Bondi solution
| ||
Range | Default: 0.0 | |
(-1.:1.) | dimensionless spin so any real number between -1 and 1
| |
bondi_coordinates | Scope: private | KEYWORD |
Description: Which coordinate system to use
| ||
Range | Default: Isotropic | |
Boyer-Lindquist | Schwarzschild or Boyer-Lindquist (Cartesian) Coordinates
| |
Kerr-Schild | Kerr-Schild (Cartesian) Coordinates
| |
Isotropic | Isotropic (Cartesian) Coordinates
| |
bondi_evolve_only_annulus | Scope: private | BOOLEAN |
Description: reset to initial data outside of bondi_freeze_inner_radius and bondi_freeze_outer_radius
| ||
Default: no | ||
bondi_freeze_inner_radius | Scope: private | REAL |
Description: reset to initial at radii below this
| ||
Range | Default: -1. | |
*:* | any value
| |
bondi_freeze_outer_radius | Scope: private | REAL |
Description: reset to initial at radii above this
| ||
Range | Default: 1e300 | |
*:* | any value
| |
bondi_overwrite_boundary | Scope: private | BOOLEAN |
Description: reset data to initial data in outer boundary in boundary condition
| ||
Default: no | ||
bondi_radial_offset | Scope: private | REAL |
Description: redefine r_grid=r_KS-r0 to avoid singularity on grid
| ||
Range | Default: 0.0 | |
0:* | Any positive number
| |
bondi_rmax | Scope: private | REAL |
Description: Largest radius in units of central mass at which the solution is found
| ||
Range | Default: 400. | |
(0:* | dimensionless outer radius for Bondi solution
| |
bondi_rmin | Scope: private | REAL |
Description: Smallest radius in units of central mass at which the solution is found
| ||
Range | Default: 1.e-15 | |
(0:* | dimensionless inner radius for Bondi solution
| |
bx_init | Scope: private | REAL |
Description: Initial B-field in the x-dir
| ||
Range | Default: 0.0 | |
*:* | Anything
| |
by_init | Scope: private | REAL |
Description: Initial B-field in the y-dir
| ||
Range | Default: 0.0 | |
*:* | Anything
| |
bz_init | Scope: private | REAL |
Description: Initial B-field in the z-dir
| ||
Range | Default: 0.0 | |
*:* | Anything
| |
change_shock_direction | Scope: private | BOOLEAN |
Description: Change the shock direction
| ||
Default: no | ||
cyl_press_inner | Scope: private | REAL |
Description: pressure in inner core
| ||
Range | Default: 1.d0 | |
(0:* | any positive number
| |
cyl_press_outer | Scope: private | REAL |
Description: pressure in outer region
| ||
Range | Default: 3.d-5 | |
(0:* | any positive number
| |
cyl_r_inner | Scope: private | REAL |
Description: Inner Radius
| ||
Range | Default: 0.8 | |
(0:* | Any positive number
| |
cyl_r_outer | Scope: private | REAL |
Description: Outer Radius
| ||
Range | Default: 1.0 | |
(0:* | Any positive number
| |
cyl_rho_inner | Scope: private | REAL |
Description: density in inner core
| ||
Range | Default: 1.d-2 | |
(0:* | any positive number
| |
cyl_rho_outer | Scope: private | REAL |
Description: density in outer region
| ||
Range | Default: 1.d-4 | |
(0:* | any positive number
| |
dens_init | Scope: private | REAL |
Description: Initial conserved mass density
| ||
Range | Default: 1.29047362 | |
(0:* | Anything positive.
| |
eps_init | Scope: private | REAL |
Description: Initial specific internal energy
| ||
Range | Default: 1.0d-6 | |
(0:* | Anything positive.
| |
gxx_init | Scope: private | REAL |
Description: Initial xx metric componenent
| ||
Range | Default: 1.0 | |
*:* | Anything, but be carefull to set a positive definite 3-metric!
| |
gxy_init | Scope: private | REAL |
Description: Initial xy metric componenent
| ||
Range | Default: 0.0 | |
*:* | Anything, but be carefull to set a positive definite 3-metric!
| |
gxz_init | Scope: private | REAL |
Description: Initial xz metric componenent
| ||
Range | Default: 0.0 | |
*:* | Anything, but be carefull to set a positive definite 3-metric!
| |
gyy_init | Scope: private | REAL |
Description: Initial yy metric componenent
| ||
Range | Default: 1.0 | |
*:* | Anything, but be carefull to set a positive definite 3-metric!
| |
gyz_init | Scope: private | REAL |
Description: Initial yz metric componenent
| ||
Range | Default: 0.0 | |
*:* | Anything, but be carefull to set a positive definite 3-metric!
| |
gzz_init | Scope: private | REAL |
Description: Initial zz metric componenent
| ||
Range | Default: 1.0 | |
*:* | Anything, but be carefull to set a positive definite 3-metric!
| |
mdot_sonicpt_bondi | Scope: private | REAL |
Description: Accretion rate at sonic point in hydro units
| ||
Range | Default: 12.566370614359172954 | |
(0:* | positive
| |
monopole_point_bx | Scope: private | REAL |
Description: Pointlike Monopole Bx value
| ||
Range | Default: 1.0 | |
*:* | Any number
| |
monopole_type | Scope: private | KEYWORD |
Description: Which kind of monopole?
| ||
Range | Default: Point | |
Point | Single point with Bx /= 0
| |
Gauss | Gaussian w/radius R_Gauss
| |
1dalt | 1-d alternating
| |
2dalt | 2-d alternating
| |
3dalt | 3-d alternating
| |
n_bondi_pts | Scope: private | REAL |
Description: Number of points to use in determining global Bondi solution
| ||
Range | Default: 2000 | |
(1:* | number of points in global Bondi solution
| |
num_bondi_sols | Scope: private | INT |
Description: Number of central masses about which to calculate Bondi solutions
| ||
Range | Default: 1 | |
1:100 | positive
| |
poloidal_a_b | Scope: private | REAL |
Description: Vector potential strength
| ||
Range | Default: 0.1 | |
*:* | Anything.
| |
poloidal_n_p | Scope: private | INT |
Description: Vector potential strength
| ||
Range | Default: 3 | |
0:* | Any positive integer.
| |
poloidal_p_cut | Scope: private | REAL |
Description: Pressure used to confine the B field inside a star
| ||
Range | Default: 1.0e-8 | |
(0:* | Anything positive.
| |
poloidal_p_p | Scope: private | INT |
Description: Index of pressure factor
| ||
Range | Default: 1 | |
(0:* | Any non-negative integer
| |
poloidal_rho_max | Scope: private | REAL |
Description: Maximum initial density
| ||
Range | Default: 1.0e-3 | |
(0:* | Anything positive.
| |
press_init | Scope: private | REAL |
Description: Initial pressure
| ||
Range | Default: 6.666666666666667d-7 | |
(0:* | Anything positive.
| |
r_gauss | Scope: private | REAL |
Description: Radius for a Gaussian monopole
| ||
Range | Default: 1.0 | |
0:* | Any positive number
| |
r_sonicpt_bondi | Scope: private | REAL |
Description: Radial distance of the sonic point from the black hole in units of mass_bh_bondi
| ||
Range | Default: 8.0 | |
(0:* | positive
| |
rho_init | Scope: private | REAL |
Description: Initial rest mass density
| ||
Range | Default: 1.0d-6 | |
(0:* | Anything positive.
| |
rotor_bvcxl | Scope: private | REAL |
Description: intial component of Bvec[0]
| ||
Range | Default: 1.0 | |
*:* | any real number
| |
rotor_bvcyl | Scope: private | REAL |
Description: intial component of Bvec[1]
| ||
Range | Default: 0.0 | |
*:* | any real number
| |
rotor_bvczl | Scope: private | REAL |
Description: intial component of Bvec[2]
| ||
Range | Default: 0.0 | |
*:* | any real number
| |
rotor_pressin | Scope: private | REAL |
Description: initial pressure inside rotor
| ||
Range | Default: 1.d0 | |
(0:* | any positive number
| |
rotor_pressout | Scope: private | REAL |
Description: initial pressure outside rotor
| ||
Range | Default: 1.d0 | |
(0:* | any positive number
| |
rotor_r_rot | Scope: private | REAL |
Description: radius of rotor
| ||
Range | Default: 0.1 | |
(0:* | any positive number
| |
rotor_rhoin | Scope: private | REAL |
Description: initial density inside rotor
| ||
Range | Default: 10.d0 | |
(0:* | any positive number
| |
rotor_rhoout | Scope: private | REAL |
Description: initial density outside rotor
| ||
Range | Default: 1.d0 | |
(0:* | any positive number
| |
rotor_rsmooth_rel | Scope: private | REAL |
Description: Define the radius in relative terms if so
| ||
Range | Default: 0.05 | |
(0:* | any positive number
| |
rotor_use_smoothing | Scope: private | BOOLEAN |
Description: Smooth the edge?
| ||
Default: yes | ||
rotor_v_max | Scope: private | REAL |
Description: Maximum velocity
| ||
Range | Default: 0.995 | |
(-1:1) | any subluminal speed (negative is clockwise)
| |
rotor_xc | Scope: private | REAL |
Description: center of rotation
| ||
Range | Default: 0.5 | |
*:* | Any location
| |
rotor_yc | Scope: private | REAL |
Description: center of rotation
| ||
Range | Default: 0.5 | |
*:* | Any location
| |
set_bondi_beta_sonicpt | Scope: private | BOOLEAN |
Description: Set plasma beta parameter instead of bondi_bmag
| ||
Default: no | ||
shock_case | Scope: private | KEYWORD |
Description: Simple, Sod’s problem or other?
| ||
Range | Default: Sod | |
Simple | GRAstro_Hydro test case
| |
Sod | Sod’s problem
| |
Blast | Strong blast wave
| |
Balsaralike1 | Hydro version of Balsara Test #1
| |
Balsara0 | Balsara Test #1, but unmagnetized
| |
Balsara1 | Balsara Test #1
| |
Balsara2 | Balsara Test #2
| |
Balsara3 | Balsara Test #3
| |
Balsara4 | Balsara Test #4
| |
Balsara5 | Balsara Test #5
| |
Alfven | Generical Alfven Test | |
Komissarov1 | Komissarov Test #1
| |
Komissarov2 | Komissarov Test #2
| |
Komissarov3 | Komissarov Test #3
| |
Komissarov4 | Komissarov Test #4
| |
Komissarov5 | Komissarov Test #5
| |
Komissarov6 | Komissarov Test #6
| |
Komissarov7 | Komissarov Test #7
| |
Komissarov8 | Komissarov Test #8
| |
Komissarov9 | Komissarov Test #9
| |
shock_radius | Scope: private | REAL |
Description: Radius of sperical shock
| ||
Range | Default: 1.0 | |
0.0:* | Anything positive
| |
shock_xpos | Scope: private | REAL |
Description: Position of shock plane: x
| ||
Range | Default: 0.0 | |
*:* | Anything
| |
shock_ypos | Scope: private | REAL |
Description: Position of shock plane: y
| ||
Range | Default: 0.0 | |
*:* | Anything
| |
shock_zpos | Scope: private | REAL |
Description: Position of shock plane: z
| ||
Range | Default: 0.0 | |
*:* | Anything
| |
shocktube_type | Scope: private | KEYWORD |
Description: Diagonal or parallel shock?
| ||
Range | Default: xshock | |
diagshock | Diagonal across all axes
| |
diagshock2d | Diagonal across x-y axes
| |
xshock | Parallel to x axis | |
yshock | Parallel to y axis
| |
zshock | Parallel to z axis | |
sphere | spherically symmetric shock
| |
simple_wave_constant_c_0 | Scope: private | REAL |
Description: The c_0 constant in Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)
| ||
Range | Default: 0.3 | |
0:1 | It is the sound speed where the fluid velocity is zero
| |
simple_wave_v_max | Scope: private | REAL |
Description: The v_max constant in Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)
| ||
Range | Default: 0.7 | |
0:1 | It is the maximum velocity in the initial configuration (see p. 1457, bottom
of first column)
| |
sx_init | Scope: private | REAL |
Description: Initial x component of conserved momentum density
| ||
Range | Default: 0.166666658 | |
*:* | Anything.
| |
sy_init | Scope: private | REAL |
Description: Initial y component of conserved momentum density
| ||
Range | Default: 0.166666658 | |
*:* | Anything.
| |
sz_init | Scope: private | REAL |
Description: Initial z component of conserved momentum density
| ||
Range | Default: 0.166666658 | |
*:* | Anything.
| |
tau_init | Scope: private | REAL |
Description: Initial conserved total energy density
| ||
Range | Default: 0.484123939 | |
(0:* | Anything positive.
| |
use_c2p_with_entropy_eqn | Scope: private | BOOLEAN |
Description: Use the con2prim routine that uses the entropy equation instead of the energy equation
| ||
Default: no | ||
velx_init | Scope: private | REAL |
Description: Initial x velocity
| ||
Range | Default: 1.0d-1 | |
*:* | Anything.
| |
vely_init | Scope: private | REAL |
Description: Initial y velocity
| ||
Range | Default: 1.0d-1 | |
*:* | Anything.
| |
velz_init | Scope: private | REAL |
Description: Initial z velocity
| ||
Range | Default: 1.0d-1 | |
*:* | Anything.
| |
bvec_evolution_method | Scope: shared from HYDROBASE | KEYWORD |
entropy_evolution_method | Scope: shared from HYDROBASE | KEYWORD |
initial_avec | Scope: shared from HYDROBASE | KEYWORD |
Extends ranges:
| ||
poloidalmagfield | Poloidal Magnetic Field
| |
initial_bvec | Scope: shared from HYDROBASE | KEYWORD |
Extends ranges:
| ||
shocktube | Shocktube type
| |
cylexp | Poloidal Magnetic Field
| |
poloidalmagfield | Poloidal Magnetic Field
| |
magnetized Bondi | radial magnetic field appropriate for Bondi test
| |
initial_entropy | Scope: shared from HYDROBASE | KEYWORD |
Extends ranges:
| ||
magnetized Bondi | Initial entropy for a radial magnetic field appropriate for Bondi test
| |
initial_hydro | Scope: shared from HYDROBASE | KEYWORD |
Extends ranges:
| ||
shocktube | Shocktube type
| |
shocktube_hot | Shocktube with hot nuclear EOS | |
only_atmo | Set only a low atmosphere
| |
read_conformal | After reading in initial alp, rho and gxx from h5 files, sets the other quantities
| |
simple_wave | Set initial data from Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)
| |
monopole | Monopole at the center
| |
cylexp | Cylindrical Explosion
| |
rotor | ”Magnetic Rotor test from DelZanna,Bucciantini , and Londrillo A&A 400,
397-413 (2003)”
| |
advectedloop | Magnetic advected loop test
| |
alfvenwave | Circularly polarized Alfven wave
| |
see [1] below | Spherical single black hole Bondi solution
| |
see [1] below | Spherical single black hole Bondi solution - TEST ISO CASE!!!!!!
| |
see [1] below | Magnetized Spherical single black hole Bondi solution
| |
see [1] below | Magnetized Spherical single black hole Bondi solution - TEST ISO CASE!!!!!!
| |
[1]
hydro\_bondi\_solution
[1]
hydro\_bondi\_solution\_iso
[1]
magnetized\_bondi\_solution
[1]
magnetized\_bondi\_solution\_iso
temperature_evolution_method | Scope: shared from HYDROBASE | KEYWORD |
timelevels | Scope: shared from HYDROBASE | INT |
y_e_evolution_method | Scope: shared from HYDROBASE | KEYWORD |
Implements:
grhydro_init_data
Inherits:
grhydro
grid
eos_omni
Group Names | Variable Names | Details | |
grhydro_init_data_reflevel | GRHydro_init_data_reflevel | compact | 0 |
description | Refinement level GRHydro is working on right now | ||
dimensions | 0 | ||
distribution | CONSTANT | ||
group type | SCALAR | ||
tags | checkpoint=”no” | ||
timelevels | 1 | ||
variable type | INT | ||
Group Names | Variable Names | Details | |
simple_wave_grid_functions | compact | 0 | |
simple_tmp | description | 1D arrays for the simple-wave routine | |
c_s | dimensions | 3 | |
distribution | DEFAULT | ||
group type | GF | ||
tags | checkpoint=”no” | ||
timelevels | 1 | ||
variable type | REAL | ||
simple_wave_scalars | compact | 0 | |
simple_rho_0 | description | values at v=0 | |
simple_eps_0 | dimensions | 0 | |
distribution | CONSTANT | ||
group type | SCALAR | ||
timelevels | 1 | ||
variable type | REAL | ||
simple_wave_output | compact | 0 | |
simple_rho | description | output variables for the simple-wave routine | |
simple_eps | dimensions | 3 | |
distribution | DEFAULT | ||
group type | GF | ||
tags | checkpoint=”no” | ||
timelevels | 1 | ||
variable type | REAL | ||
Uses header:
SpaceMask.h
This section lists all the variables which are assigned storage by thorn EinsteinEvolve/GRHydro_InitData. Storage can either last for the duration of the run (Always means that if this thorn is activated storage will be assigned, Conditional means that if this thorn is activated storage will be assigned for the duration of the run if some condition is met), or can be turned on for the duration of a schedule function.
Conditional: | |
GRHydro_init_data_reflevel | |
GRHydro_init_data_reflevel | |
GRHydro_init_data_reflevel | |
GRHydro_init_data_reflevel | |
simple_wave_grid_functions | |
simple_wave_scalars | |
simple_wave_output | |
CCTK_PARAMCHECK
grhydro_initdata_checkparameters
check parameters
Language: | c | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_monopolem
monopole initial data
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_shocktube
shocktube initial data
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_cylindricalexplosionm
cylindrical explosion initial data - mhd-only
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_init_data_refinementlevel
calculate current refinement level
Before: | grhydro_con2primtest | |
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_con2primtest
testing the conservative to primitive solver
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_init_data_refinementlevel
calculate current refinement level
Before: | c2p2c_call | |
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
c2p2cm
testing conservative to primitive to conservative - mhd version
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
c2p2c
testing conservative to primitive to conservative
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_init_data_refinementlevel
calculate current refinement level
Before: | p2c2p_call | |
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
p2c2pm
testing primitive to conservative to primitive - mhd version
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
p2c2p
testing primitive to conservative to primitive
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_rotorm
mhd rotor initial data
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_init_data_refinementlevel
calculate current refinement level
Before: | p2c2p_call | |
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
p2c2pm_polytype
testing primitive to conservative to primitive - mhd polytype version
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_reconstruction_test
testing the reconstruction
Language: | fortran | |
Options: | global | |
loop-local | ||
Storage: | grhydro_prim_bext | |
grhydro_scalars | ||
Type: | function | |
HydroBase_Initial (conditional)
grhydro_only_atmo
only atmosphere as initial data
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_readconformaldata
set the missing quantities, after reading in from file initial data from conformally-flat codes (garching)
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_simplewave
set initial data from anile miller motta, phys.fluids. 26, 1450 (1983)
Language: | fortran | |
Type: | function | |
CCTK_ANALYSIS (conditional)
grhydro_simplewave_analysis
compute some output variables for the simple wave
After: | grhydro_entropy | |
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_bondi_iso
setup grhydro vars for the hydrodynamic bondi solution
After: | hydrobase_excisionmasksetup | |
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_bondim_iso
setup grhydro vars for the magnetized bondi solution
After: | hydrobase_excisionmasksetup | |
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_bondi
setup grhydro vars for the hydrodynamic bondi solution
After: | hydrobase_excisionmasksetup | |
Language: | c | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_advectedloopm
mhd advected loop initial data
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_bondim
setup grhydro vars for the magnetized bondi solution
After: | hydrobase_excisionmasksetup | |
Language: | c | |
Type: | function | |
HydroBase_Con2Prim (conditional)
grhydro_bondim_range
force analytic solution outside anulus
Before: | con2prim | |
Language: | c | |
Type: | function | |
HydroBase_Boundaries (conditional)
grhydro_bondim_boundary
force analytic solution in boundaries
Before: | hydrobase_select_boundaries | |
Language: | c | |
Type: | function | |
CCTK_INITIAL (conditional)
grhydro_poloidalmagfieldm
set up a poloidal magnetic field. it expects the other fluid variables already to be set, as for example in the tov solution
After: | hydrobase_initial | |
Before: | grhydrotransformprimtolocalbasis | |
Language: | fortran | |
Type: | function | |
CCTK_INITIAL (conditional)
hydrobase_boundaries
call boundary conditions after magnetic field initial data setup
After: | grhydro_poloidalmagfieldm | |
Before: | grhydrotransformprimtolocalbasis | |
Type: | group | |
HydroBase_Initial (conditional)
grhydro_alfvenwavem
circularly polarized alfven wave initial data
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_shocktube_hot
hot shocktube initial data
After: | hydrobase_y_e_one | |
hydrobase_zero | ||
Language: | fortran | |
Type: | function | |
HydroBase_Initial (conditional)
grhydro_shocktubem
shocktube initial data - mhd version
Language: | fortran | |
Type: | function | |
ApplyBCs (conditional)
grhydro_diagshock_boundarym
diagonal shock boundary conditions
After: | boundaryconditions | |
boundary::boundary_clearselection | ||
Language: | fortran | |
Sync: | grhydro::dens | |
grhydro::tau | ||
grhydro::scon | ||
grhydro::bcons | ||
grhydro::psidc | ||
Type: | function | |
ApplyBCs (conditional)
grhydro_diagshock_boundarym
diagonal shock boundary conditions
After: | boundaryconditions | |
boundary::boundary_clearselection | ||
Language: | fortran | |
Sync: | grhydro::dens | |
grhydro::tau | ||
grhydro::scon | ||
grhydro::bcons | ||
Type: | function | |
ApplyBCs (conditional)
grhydro_diagshock2d_boundarym
2-d diagonal shock boundary conditions
After: | boundaryconditions | |
boundary::boundary_clearselection | ||
Language: | fortran | |
Type: | function | |
Alias Name: | Function Name: |
c2p2c | c2p2c_call |
c2p2cM | c2p2c_call |
p2c2p | p2c2p_call |
p2c2pM | p2c2p_call |
p2c2pM_polytype | p2c2p_call |