read – SPHINCS

Implementation of the methods of TYPE bns that import BNS data using $\texttt{LORENE}$

FT 25.10.2021

Uses

- utility
- iso_c_binding
- Ancestors:
- diffstar_lorene
Help

Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.

Module Procedures

is_hydro_positive read_drslorene_mass_density read_drslorene_pressure read_id_full read_id_hydro read_id_int read_id_k read_id_mass_b read_id_particles read_id_spacetime read_spatial_metric

Module Procedures

module procedure is_hydro_positive module function is_hydro_positive(this, x, y, z) result(res)

Return .FALSE. if the energy density is nonpositive or if the specific energy is nonpositive, or if the pressure is nonpositive at the specified point; .TRUE. otherwise

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(in)		::	this	diffstarlorene object which this PROCEDURE is a member of
double precision,	intent(in),	VALUE	::	x	$x$ coordinate of the desired point
double precision,	intent(in),	VALUE	::	y	$y$ coordinate of the desired point
double precision,	intent(in),	VALUE	::	z	$z$ coordinate of the desired point

Return Value logical

1 if the energy density or the specific energy or the pressure are negative, 0 otherwise

module procedure read_drslorene_mass_density module function read_drslorene_mass_density(this, x, y, z) result(res)

Returns the $\texttt{LORENE}$ mass density at the point given as argument, in units of $M_\odot/L_\odot^3$ .

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(in)		::	this	diffstarlorene object which this PROCEDURE is a member of
double precision,	intent(in),	VALUE	::	x	$x$ coordinate of the desired point
double precision,	intent(in),	VALUE	::	y	$y$ coordinate of the desired point
double precision,	intent(in),	VALUE	::	z	$z$ coordinate of the desired point

Return Value double precision

Baryon mass density at $(x,y,z)$

module procedure read_drslorene_pressure module function read_drslorene_pressure(this, x, y, z) result(res)

Returns the $\texttt{LORENE}$ pressure at the point given as argument, in units of $[\mathrm{kg}\,c^2\, \mathrm{m}^{-3}]$ .

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(in)		::	this	bnslorene object which this PROCEDURE is a member of
double precision,	intent(in),	VALUE	::	x	$x$ coordinate of the desired point
double precision,	intent(in),	VALUE	::	y	$y$ coordinate of the desired point
double precision,	intent(in),	VALUE	::	z	$z$ coordinate of the desired point

Return Value double precision

Pressure at $(x,y,z)$

module procedure read_id_full module subroutine read_id_full(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z)

Stores the ID in non-diffstarlorene-member arrays with the same shape as the diffstarlorene member arrays

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(inout)		::	this	diffstarlorene object which this PROCEDURE is a member of
integer,	intent(in)		::	n
double precision,	intent(in),	DIMENSION(:)	::	x
double precision,	intent(in),	DIMENSION(:)	::	y
double precision,	intent(in),	DIMENSION(:)	::	z
double precision,	intent(inout),	DIMENSION(:)	::	lapse
double precision,	intent(inout),	DIMENSION(:)	::	shift_x
double precision,	intent(inout),	DIMENSION(:)	::	shift_y
double precision,	intent(inout),	DIMENSION(:)	::	shift_z
double precision,	intent(inout),	DIMENSION(:)	::	g_xx
double precision,	intent(inout),	DIMENSION(:)	::	g_xy
double precision,	intent(inout),	DIMENSION(:)	::	g_xz
double precision,	intent(inout),	DIMENSION(:)	::	g_yy
double precision,	intent(inout),	DIMENSION(:)	::	g_yz
double precision,	intent(inout),	DIMENSION(:)	::	g_zz
double precision,	intent(inout),	DIMENSION(:)	::	k_xx
double precision,	intent(inout),	DIMENSION(:)	::	k_xy
double precision,	intent(inout),	DIMENSION(:)	::	k_xz
double precision,	intent(inout),	DIMENSION(:)	::	k_yy
double precision,	intent(inout),	DIMENSION(:)	::	k_yz
double precision,	intent(inout),	DIMENSION(:)	::	k_zz
double precision,	intent(inout),	DIMENSION(:)	::	baryon_density
double precision,	intent(inout),	DIMENSION(:)	::	energy_density
double precision,	intent(inout),	DIMENSION(:)	::	specific_energy
double precision,	intent(inout),	DIMENSION(:)	::	pressure
double precision,	intent(inout),	DIMENSION(:)	::	u_euler_x
double precision,	intent(inout),	DIMENSION(:)	::	u_euler_y
double precision,	intent(inout),	DIMENSION(:)	::	u_euler_z

module procedure read_id_hydro module subroutine read_id_hydro(this, nx, ny, nz, pos, baryon_density, energy_density, specific_energy, pressure, u_euler)

Stores the hydro ID in the arrays needed to compute the constraints on the refined mesh

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(inout)		::	this	diffstarlorene object which this PROCEDURE is a member of
integer,	intent(in)		::	nx
integer,	intent(in)		::	ny
integer,	intent(in)		::	nz
double precision,	intent(in),	DIMENSION(:,:,:,:)	::	pos
double precision,	intent(inout),	DIMENSION(:,:,:)	::	baryon_density
double precision,	intent(inout),	DIMENSION(:,:,:)	::	energy_density
double precision,	intent(inout),	DIMENSION(:,:,:)	::	specific_energy
double precision,	intent(inout),	DIMENSION(:,:,:)	::	pressure
double precision,	intent(inout),	DIMENSION(:,:,:,:)	::	u_euler

module procedure read_id_int module subroutine read_id_int(this, n, x, y, z)

Stores the ID in the diffstarlorene member arrays

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(inout)		::	this	diffstarlorene object which this PROCEDURE is a member of
integer,	intent(in)		::	n
double precision,	intent(in),	DIMENSION(:)	::	x
double precision,	intent(in),	DIMENSION(:)	::	y
double precision,	intent(in),	DIMENSION(:)	::	z

module procedure read_id_k module subroutine read_id_k(this, n, x, y, z, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz)

Stores the components of the extrinsic curvature in arrays

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(inout)		::	this	diffstarlorene object which this PROCEDURE is a member of
integer,	intent(in)		::	n
double precision,	intent(in),	DIMENSION(:)	::	x
double precision,	intent(in),	DIMENSION(:)	::	y
double precision,	intent(in),	DIMENSION(:)	::	z
double precision,	intent(inout),	DIMENSION(:)	::	k_xx
double precision,	intent(inout),	DIMENSION(:)	::	k_xy
double precision,	intent(inout),	DIMENSION(:)	::	k_xz
double precision,	intent(inout),	DIMENSION(:)	::	k_yy
double precision,	intent(inout),	DIMENSION(:)	::	k_yz
double precision,	intent(inout),	DIMENSION(:)	::	k_zz

module procedure read_id_mass_b module subroutine read_id_mass_b(this, x, y, z, g, baryon_density, gamma_euler)

Stores the hydro ID in the arrays needed to compute the baryon mass, storing it to variables (not arrays as the others SUBROUTINES in this SUBMODULE).

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(inout)		::	this	diffstarlorene object which this PROCEDURE is a member of
double precision,	intent(in)		::	x
double precision,	intent(in)		::	y
double precision,	intent(in)		::	z
double precision,	intent(out),	DIMENSION(6)	::	g
double precision,	intent(out)		::	baryon_density
double precision,	intent(out)		::	gamma_euler

module procedure read_id_particles module subroutine read_id_particles(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z)

Stores the hydro ID in the arrays needed to compute the SPH ID

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(inout)		::	this	diffstarlorene object which this PROCEDURE is a member of
integer,	intent(in)		::	n
real(kind=C_DOUBLE),	intent(in),	DIMENSION(:)	::	x
real(kind=C_DOUBLE),	intent(in),	DIMENSION(:)	::	y
real(kind=C_DOUBLE),	intent(in),	DIMENSION(:)	::	z
double precision,	intent(inout),	DIMENSION(:)	::	lapse
double precision,	intent(inout),	DIMENSION(:)	::	shift_x
double precision,	intent(inout),	DIMENSION(:)	::	shift_y
double precision,	intent(inout),	DIMENSION(:)	::	shift_z
double precision,	intent(inout),	DIMENSION(:)	::	g_xx
double precision,	intent(inout),	DIMENSION(:)	::	g_xy
double precision,	intent(inout),	DIMENSION(:)	::	g_xz
double precision,	intent(inout),	DIMENSION(:)	::	g_yy
double precision,	intent(inout),	DIMENSION(:)	::	g_yz
double precision,	intent(inout),	DIMENSION(:)	::	g_zz
double precision,	intent(inout),	DIMENSION(:)	::	baryon_density
double precision,	intent(inout),	DIMENSION(:)	::	energy_density
double precision,	intent(inout),	DIMENSION(:)	::	specific_energy
double precision,	intent(inout),	DIMENSION(:)	::	pressure
double precision,	intent(inout),	DIMENSION(:)	::	u_euler_x
double precision,	intent(inout),	DIMENSION(:)	::	u_euler_y
double precision,	intent(inout),	DIMENSION(:)	::	u_euler_z

module procedure read_id_spacetime module subroutine read_id_spacetime(this, nx, ny, nz, pos, lapse, shift, g, ek)

Stores the spacetime ID in multi-dimensional arrays needed to compute the BSSN variables and constraints

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(inout)		::	this	diffstarlorene object which this PROCEDURE is a member of
integer,	intent(in)		::	nx
integer,	intent(in)		::	ny
integer,	intent(in)		::	nz
double precision,	intent(in),	DIMENSION(:,:,:,:)	::	pos
double precision,	intent(inout),	DIMENSION(:,:,:)	::	lapse
double precision,	intent(inout),	DIMENSION(:,:,:,:)	::	shift
double precision,	intent(inout),	DIMENSION(:,:,:,:)	::	g
double precision,	intent(inout),	DIMENSION(:,:,:,:)	::	ek

module procedure read_spatial_metric module function read_spatial_metric(this, x, y, z) result(res)

Returns the $\texttt{LORENE}$ conformal factor to the 4th power, equal to the diagonal components of the conformally flat spatial ADM metric.

Arguments

Type	Intent	Attributes		Name
class(diffstarlorene),	intent(in)		::	this	diffstarlorene object which this PROCEDURE is a member of
real(kind=C_DOUBLE),	intent(in),	VALUE	::	x	$x$ coordinate of the desired point
real(kind=C_DOUBLE),	intent(in),	VALUE	::	y	$y$ coordinate of the desired point
real(kind=C_DOUBLE),	intent(in),	VALUE	::	z	$z$ coordinate of the desired point

Return Value real(kind=C_DOUBLE)

$g_{xx}=g_{yy}=g_{zz}$ at $(x,y,z)$

read Submodule

Contents

Module Procedures

Uses

Ancestors:

Contents

Module Procedures

Module Procedures

module procedure is_hydro_positive module function is_hydro_positive(this, x, y, z) result(res)

Arguments

Return Value logical

module procedure read_drslorene_mass_density module function read_drslorene_mass_density(this, x, y, z) result(res)

Arguments

Return Value double precision

module procedure read_drslorene_pressure module function read_drslorene_pressure(this, x, y, z) result(res)

Arguments

Return Value double precision

module procedure read_id_full module subroutine read_id_full(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z)

Arguments

module procedure read_id_hydro module subroutine read_id_hydro(this, nx, ny, nz, pos, baryon_density, energy_density, specific_energy, pressure, u_euler)

Arguments

module procedure read_id_int module subroutine read_id_int(this, n, x, y, z)

Arguments

module procedure read_id_k module subroutine read_id_k(this, n, x, y, z, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz)

Arguments

module procedure read_id_mass_b module subroutine read_id_mass_b(this, x, y, z, g, baryon_density, gamma_euler)

Arguments

module procedure read_id_particles module subroutine read_id_particles(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z)

Arguments

module procedure read_id_spacetime module subroutine read_id_spacetime(this, nx, ny, nz, pos, lapse, shift, g, ek)

Arguments

module procedure read_spatial_metric module function read_spatial_metric(this, x, y, z) result(res)

Arguments

Return Value real(kind=C_DOUBLE)

read Submodule

Contents

Module Procedures

Uses

Ancestors:

Contents

Module Procedures

Module Procedures

module procedure is_hydro_positive module function is_hydro_positive(this, x, y, z) result(res) Interface →

Arguments

Return Value logical

module procedure read_drslorene_mass_density module function read_drslorene_mass_density(this, x, y, z) result(res) Interface →

Arguments

Return Value double precision

module procedure read_drslorene_pressure module function read_drslorene_pressure(this, x, y, z) result(res) Interface →

Arguments

Return Value double precision

module procedure read_id_full module subroutine read_id_full(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z) Interface →

Arguments

module procedure read_id_hydro module subroutine read_id_hydro(this, nx, ny, nz, pos, baryon_density, energy_density, specific_energy, pressure, u_euler) Interface →

Arguments

module procedure read_id_int module subroutine read_id_int(this, n, x, y, z) Interface →

Arguments

module procedure read_id_k module subroutine read_id_k(this, n, x, y, z, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz) Interface →

Arguments

module procedure read_id_mass_b module subroutine read_id_mass_b(this, x, y, z, g, baryon_density, gamma_euler) Interface →

Arguments

module procedure read_id_particles module subroutine read_id_particles(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z) Interface →

Arguments

module procedure read_id_spacetime module subroutine read_id_spacetime(this, nx, ny, nz, pos, lapse, shift, g, ek) Interface →

Arguments

module procedure read_spatial_metric module function read_spatial_metric(this, x, y, z) result(res) Interface →

Arguments

Return Value real(kind=C_DOUBLE)

module procedure is_hydro_positive module function is_hydro_positive(this, x, y, z) result(res)

module procedure read_drslorene_mass_density module function read_drslorene_mass_density(this, x, y, z) result(res)

module procedure read_drslorene_pressure module function read_drslorene_pressure(this, x, y, z) result(res)

module procedure read_id_full module subroutine read_id_full(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z)

module procedure read_id_hydro module subroutine read_id_hydro(this, nx, ny, nz, pos, baryon_density, energy_density, specific_energy, pressure, u_euler)

module procedure read_id_int module subroutine read_id_int(this, n, x, y, z)

module procedure read_id_k module subroutine read_id_k(this, n, x, y, z, k_xx, k_xy, k_xz, k_yy, k_yz, k_zz)

module procedure read_id_mass_b module subroutine read_id_mass_b(this, x, y, z, g, baryon_density, gamma_euler)

module procedure read_id_particles module subroutine read_id_particles(this, n, x, y, z, lapse, shift_x, shift_y, shift_z, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz, baryon_density, energy_density, specific_energy, pressure, u_euler_x, u_euler_y, u_euler_z)

module procedure read_id_spacetime module subroutine read_id_spacetime(this, nx, ny, nz, pos, lapse, shift, g, ek)

module procedure read_spatial_metric module function read_spatial_metric(this, x, y, z) result(res)