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 \(\rho _{_L}\) \(v^i_{_L}\) \(\varepsilon _{_L}\) \(\rho _{_R}\) \(v^i_{_R}\) \(\varepsilon _{_R}\)







Sod 1 0 1.5 0.12500.15
Simple100 20 1 0 \(10^{-6}\)
Blast 1 01500 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].

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 Aˆb



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ˆ3  =0, where Bˆ3  —— 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