GRHydro_Init_Data

Luca Baiotti, Ian Hawke, Scott Hawley

24/07/2008

Abstract

GRHydro_Init_Data - some initial data for GRHydro

1 Introduction

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 .

1.1 Tests

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.

1.2 Shocktube tests

There are three possible shock-tube problems, referred to as Sod, Simple and Blast, with initial data

Type ρLvLi𝜀 LρRvRi𝜀 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 ⋅ 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].

1.3 Only atmosphere

For testing purposes, this routine sets all the points to the values of the atmosphere.

1.4 Simple wave

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).

References

[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.

2 Parameters

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 Ab^
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 (B3^=0, where B3^ —— 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

3 Interfaces

General

Implements:

grhydro_init_data

Inherits:

grhydro

grid

eos_omni

Grid Variables

3.0.1 PRIVATE GROUPS
  Group Names     Variable Names     Details   
grhydro_init_data_reflevelGRHydro_init_data_reflevel  compact0
  descriptionRefinement level GRHydro is working on right now
  dimensions0
  distribution CONSTANT
  group typeSCALAR
  tagscheckpoint=”no”
  timelevels1
  variable type INT

3.0.2 PROTECTED GROUPS
  Group Names     Variable Names    Details   
simple_wave_grid_functions   compact0
simple_tmp   description1D arrays for the simple-wave routine
c_s   dimensions3
  distribution DEFAULT
  group typeGF
  tagscheckpoint=”no”
  timelevels1
  variable type REAL
simple_wave_scalars   compact0
simple_rho_0   descriptionvalues at v=0
simple_eps_0   dimensions0
  distributionCONSTANT
  group type SCALAR
  timelevels1
  variable type REAL
simple_wave_output   compact0
simple_rho   descriptionoutput variables for the simple-wave routine
simple_eps   dimensions3
  distributionDEFAULT
  group type GF
  tagscheckpoint=”no”
  timelevels1
  variable type REAL

Uses header:

SpaceMask.h

4 Schedule

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.

Storage

 

 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
   

Scheduled Functions

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

Aliased Functions

 

Alias Name:        Function Name:
c2p2c c2p2c_call
c2p2cM c2p2c_call
p2c2p p2c2p_call
p2c2pM p2c2p_call
p2c2pM_polytype p2c2p_call