bnsfuka Derived Type

public interface destruct_bnsfuka()

Arguments

interface

public module subroutine allocate_bnsfuka_memory(this, d) Implementation →

Allocates allocatable arrays member of a bnsfuka object

Arguments

interface

public module subroutine check_i_matter(this, i_matter) Implementation →

Checks that the given index i_matter is between 1 and n_matter, included. If not, it stops the execution of the program.

Arguments

interface

public module subroutine construct_binary(this, fukafile) Implementation →

Interface of the subroutine that constructs the \(\texttt{FUKA}\) \(\texttt{Bin_NS}\) object

Arguments

interface

public module subroutine deallocate_bnsfuka_memory(this) Implementation →

Deallocates allocatable arrays member of a bnsfuka object

Arguments

interface

public module subroutine construct_bnsfuka(derived_type, filename, eos_filenames) Implementation →

Prints a summary of the physical properties the system to the standard output and, optionally, to a formatted file whose name is given as the optional argument filename

Arguments

interface

public module subroutine destruct_binary(this) Implementation →

Destructs a \(\texttt{FUKA}\) \(\texttt{Bin_NS}\) object

Arguments

interface

public module function estimate_lengthscale_field(this, get_field, n_mat) result(scales) Implementation →

Estimate typical length scales, one per each matter object, by computing $\dfrac{f}{\partial f}$, where $f$ is a field given as input, and $\partial$ represent a derivative of it. Presently, the derivatives are computed separately along each spatial dimension, as 1D derivatives.

Arguments

Return Value double precision, DIMENSION(n_mat)

Array of the minimum $\dfrac{f}{\partial f}$ over the lattices that surround each matter object

interface

public module subroutine finalize(npart, pos, nlrf, u, pr, vel_u, theta, nstar, nu) Implementation →

Post-process the \(\mathrm{SPH}\) \(\mathrm{ID}\); for example, correct for the residual ADM linear momentum.

Arguments

interface

public module function find_center(this, separation, x_sign, get_density) result(center) Implementation →

Finds the center of a star, as the point where the density is maximal.

Arguments

Return Value double precision

Center of a star, as the point where the density is maximal.

interface

public module subroutine find_print_surfaces(this) Implementation →

Finds the surfaces of the stars, and prints them to a formatted file.

Arguments

interface

public module function find_radius(this, center, vector, get_density) result(radius) Implementation →

Finds the radius of a matter object, relative to a center and along a direction. The radius is determined as the first point where the density is zero.

Arguments

Return Value double precision

Radius of the star relative to center, along the direction specified by vector

interface

public module subroutine find_surface(this, center, n_theta, n_phi, surface) Implementation →

Finds the surface of a star, using [[bnsbase::find_radius]] along many directions.

Arguments

interface

public pure module function get_adm_mass(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_angular_momentum(this) Implementation →

Arguments

Return Value double precision, (3)

Result

interface

public pure module function get_angular_vel(this) Implementation →

Returns angular_vel

Arguments

Return Value double precision

Result

interface

public pure module function get_barycenter1_x(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_barycenter2_x(this) Implementation →

Arguments

Return Value double precision

Result

interface

public module function get_bns_identifier(this) Implementation →

Arguments

Return Value double precision

interface

public pure module function get_center1_x(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_center2_x(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_cold_system(this) Implementation →

Returns cold_system, the LOGICAL variable at specifies if the system is cold (no thermal component)

Arguments

Return Value logical

cold_system

interface

public pure module function get_distance(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_distance_com(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_energy_density_center1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_energy_density_center2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_ent_center1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_ent_center2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_eos1(this) Implementation →

Arguments

Return Value character(len=:), ALLOCATABLE

Result

interface

public module function get_eos1_fukaid(this) Implementation →

Arguments

Return Value integer

interface

public module function get_eos1_fukaid(this) Implementation →

Arguments

Return Value integer

interface

public pure module function get_eos2(this) Implementation →

Arguments

Return Value character(len=:), ALLOCATABLE

Result

interface

public module function get_eos2_fukaid(this) Implementation →

Arguments

Return Value integer

interface

public module function get_eos2_fukaid(this) Implementation →

Arguments

Return Value integer

interface

public module function get_eos_id(this, i_matter) Implementation →

Arguments

Return Value integer

Result

interface

public pure module function get_estimate_length_scale(this) Implementation →

Returns estimate_length_scale, the LOGICAL variable that specifies if a typical length scale, equal to the ratio of a field over its gradient, should be computed

Arguments

Return Value logical

estimate_length_scale

interface

public module function get_field_array(this, field) result(field_array) Implementation →

Returns the bnsfuka member arrays named field

Arguments

Return Value double precision, DIMENSION(:), ALLOCATABLE

Desired bnsfuka member array

public interface get_field_array()

Arguments

public interface get_field_value()

Arguments

interface

public module function get_field_value(this, field, n) result(field_value) Implementation →

Returns the component n of the bnsfuka member arrays named field

Arguments

Return Value double precision

Component n of the desired bnsfuka member array

interface

public pure module function get_gamma0_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma0_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma1_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma1_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma2_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma2_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma3_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma3_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_gamma_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_grav_mass1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_grav_mass2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa0_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa0_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa1_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa1_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa2_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa2_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa3_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa3_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_kappa_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_linear_momentum(this) Implementation →

Arguments

Return Value double precision, (3)

Result

interface

public pure module function get_logP1_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_logP1_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_logRho0_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_logRho0_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_logRho1_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_logRho1_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_logRho2_1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_logRho2_2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_mass1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_mass2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_n_matter(this) Implementation →

Returns n_matter, the number of matter objects in the physical system

Arguments

Return Value integer

n_matter, the number of matter objects in the

interface

public pure module function get_nbar_center1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_nbar_center2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_npeos_1(this) Implementation →

Arguments

Return Value integer

Result

interface

public pure module function get_npeos_2(this) Implementation →

Arguments

Return Value integer

Result

interface

public pure module function get_one_lapse(this) Implementation →

Returns n_matter, the number of matter objects in the physical system

Arguments

Return Value logical

n_matter, the number of matter objects in the

interface

public pure module function get_pressure_center1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_pressure_center2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_radius1_x_comp(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_radius1_x_opp(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_radius1_y(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_radius1_z(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_radius2_x_comp(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_radius2_x_opp(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_radius2_y(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_radius2_z(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_rho_center1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_rho_center2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_specific_energy_center1(this) Implementation →

Arguments

Return Value double precision

Result

interface

public pure module function get_specific_energy_center2(this) Implementation →

Arguments

Return Value double precision

Result

interface

public module function get_total_spatial_extent(this) result(box) Implementation →

INTERFACE to the SUBROUTINE that detects the spatial extent of the physical system considered, and returns a 6-dimensional array containing the coordinates $x_{\rm min},x_{\rm max},y_{\rm min},y_{\rm max},z_{\rm min},z_{\rm max}$ of a box centered at the center of the object and containing the system.

Arguments

Return Value double precision, DIMENSION(6)

6-dimensional array containing the coordinates $x_{\rm min},x_{\rm max},y_{\rm min},y_{\rm max}, z_{\rm min},z_{\rm max}$ of a box containing the physical system.

interface

public pure module function get_zero_shift(this) Implementation →

Returns n_matter, the number of matter objects in the physical system

Arguments

Return Value logical

n_matter, the number of matter objects in the

interface

public module subroutine initialize(derived_type, filename, eos_filenames) Implementation →

Initialize the DERIVED TYPE that extends idbase

Arguments

interface

public recursive module subroutine initialize_id_bnsfuka(this, flag, switch) Implementation →

Initialize the \(\texttt{FUKA}\) \(\mathrm{BNS}\) \(\mathrm{ID}\).

• If flag= [[utility:flag$sph]], set up the lattices around the stars for the \(\mathrm{BNS}\) produced with \(\texttt{FUKA}\).
• If flag= [[utility:flag$tpo]], allocate memory for the hydro grid functions.
• If flag > 0, assign its value to l_curr.
• If [[utility:flag$tpo]] < flag < 0, assign its value to tpo_curr.

Arguments

interface

public module subroutine integrate_baryon_mass_density(this, center, radius, central_density, dr, dth, dphi, mass, mass_profile, mass_profile_idx, radii, surf) Implementation →

INTERFACE to the SUBROUTINE integrating the baryon mass density to compute the radial mass profile of a single star.

Arguments

interface

public module function interpolate_fuka_spatial_metric(this, x, y, z) result(res) Implementation →

Returns the \(\texttt{FUKA}\) conformally flat spatial metric component $g_{xx}=g_{yy}=g_{zz}$ at a point $(x,y,z)$

Arguments

Return Value real(kind=C_DOUBLE)

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

interface

public module function is_hydro_positive(this, x, y, z) result(res) Implementation →

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

Arguments

Return Value logical

.TRUE. if the energy density or the specific energy or the pressure are negative, .FALSE. otherwise

interface

public module subroutine print_bns_properties(this) Implementation →

Prints the \(\mathrm{BNS}\) parameters to the standard output

Arguments

interface

public module subroutine print_summary_bnsbase(this, filename) Implementation →

Prints a summary of the physical properties the \(\mathrm{BNS}\) system to the standard output and, optionally, to a formatted file whose name is given as the optional argument filename

Arguments

interface

public module subroutine print_summary_bnsfuka(this, filename) Implementation →

Prints a summary of the physical properties of the \(\mathrm{BNS}\) produced by \(\texttt{FUKA}\) to the standard output and, optionally, to a formatted file whose name is given as the optional argument filename

Arguments

interface

public module subroutine read_bns_properties(this) Implementation →

Imports the \(\mathrm{BNS}\) parameters from \(\texttt{FUKA}\)

Arguments

interface

public module subroutine read_fuka_id_mass_b(this, x, y, z, g, baryon_density, gamma_euler) Implementation →

Stores the hydro \(\mathrm{ID}\) in the arrays needed to compute the baryon mass

Arguments

interface

public module subroutine read_fuka_id_member(this, n, x, y, z) Implementation →

Stores the \(\mathrm{ID}\) in the bnsfuka member arrays

Arguments

interface

public module subroutine read_fuka_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) Implementation →

Stores the hydro \(\mathrm{ID}\) in the arrays needed to compute the \(\mathrm{SPH}\) \(\mathrm{ID}\)

Arguments

interface

public module function read_fuka_mass_density(this, x, y, z) result(res) Implementation →

Returns the \(\texttt{FUKA}\) baryon mass density at a point $(x,y,z)$

Arguments

Return Value double precision

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

interface

public module function read_fuka_pressure(this, x, y, z) result(res) Implementation →

Returns the \(\texttt{FUKA}\) pressure at a point $(x,y,z)$

Arguments

Return Value double precision

Pressure at $(x,y,z)$

interface

public module function read_fuka_spatial_metric(this, x, y, z) result(res) Implementation →

Returns the \(\texttt{FUKA}\) conformally flat spatial metric component $g_{xx}=g_{yy}=g_{zz}$ at a point $(x,y,z)$

Arguments

Return Value real(kind=C_DOUBLE)

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

interface

public module subroutine read_fuka_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) Implementation →

Stores the \(\mathrm{ID}\) in non bnsfuka-member arrays with the same shape as the bnsfuka member arrays

Arguments