place_particles_lattice

module procedure place_particles_lattice module subroutine place_particles_lattice(this, central_density, xmin, xmax, ymin, ymax, zmin, zmax, npart_des, npart_out, stretch, thres, pos, pvol, nu, h, get_density, get_id, validate_position)

Uses

- constants
- utility
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.

Places paricles on a 3D lattice containing both stars

FT 5.10.2020

Arguments

Type Intent Optional Attributes Name

class(particles),

intent(inout)

this

particles object which this PROCEDURE is a member of

double precision,

intent(in)

central_density

Maximum baryon mass density of the system

double precision,

intent(in)

xmin

Left $x$ boundary of the lattice

double precision,

intent(in)

xmax

Right $x$ boundary of the lattice

double precision,

intent(in)

ymin

Left $y$ boundary of the lattice

double precision,

intent(in)

ymax

Right $y$ boundary of the lattice

double precision,

intent(in)

zmin

Left $z$ boundary of the lattice

double precision,

intent(in)

zmax

Right $z$ boundary of the lattice

integer,

intent(in)

npart_des

Desired particle number

integer,

intent(out)

npart_out

Real, output particle number

double precision,

intent(in)

stretch

Stretching factor fo the lattice. xmin to zmax are multiplied by it

double precision,

intent(in)

thres

(~rho_max)/thres is the minimum mass density considered when placing particles. Used only when redistribute_nu is .FALSE. . When redistribute_nu is .TRUE. thres= 100*nu_ratio

double precision, intent(inout), DIMENSION(:,:), ALLOCATABLE :: pos

double precision,

intent(inout),

DIMENSION(:), ALLOCATABLE

pvol

Array storing the particle positions Array storing the particle volumes

double precision,

intent(inout),

DIMENSION(:), ALLOCATABLE

Array storing the particle baryon masses

double precision,

intent(inout),

DIMENSION(:), ALLOCATABLE

Array storing the initial guess for the particle smoothing lengths

function get_density(x, y, z) result(density)

Returns the baryon mass density at the desired point

Arguments

Type	Intent		Name
double precision,	intent(in)	::	x	$x$ coordinate of the desired point
double precision,	intent(in)	::	y	$y$ coordinate of the desired point
double precision,	intent(in)	::	z	$z$ coordinate of the desired point

Return Value double precision

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

subroutine get_id(x, y, z, sqdetg, baryon_density, gamma_euler)

Returns the baryon mass density at the desired point

Arguments

Type	Intent		Name
double precision,	intent(in)	::	x	$x$ coordinate of the desired point
double precision,	intent(in)	::	y	$y$ coordinate of the desired point
double precision,	intent(in)	::	z	$z$ coordinate of the desired point
double precision,	intent(out)	::	sqdetg
double precision,	intent(out)	::	baryon_density
double precision,	intent(out)	::	gamma_euler

procedure(validate_position_int),

optional

validate_position

Returns 1 if the position is not valid, 0 otherwise

Calls

Help

Called by

Help

Variables

Type	Visibility	Attributes		Name
double precision,	private,	DIMENSION(:,:,:), ALLOCATABLE	::	bar_density_tmp
double precision,	private		::	dx
double precision,	private		::	dy
double precision,	private		::	dz
double precision,	private,	DIMENSION(:,:,:), ALLOCATABLE	::	gam_euler_tmp
integer,	private		::	i
integer,	private		::	j
integer,	private		::	k
integer,	private		::	npart_half
integer,	private		::	npart_tmp
double precision,	private,	DIMENSION(:), ALLOCATABLE	::	nstar
double precision,	private,	DIMENSION(:,:,:), ALLOCATABLE	::	nstar_tmp
integer,	private		::	nx
integer,	private		::	ny
integer,	private		::	nz
double precision,	private,	DIMENSION(:,:,:,:), ALLOCATABLE	::	pos_tmp
integer,	private		::	sgn
double precision,	private,	DIMENSION(:,:,:), ALLOCATABLE	::	sqdetg_tmp
double precision,	private		::	thres_baryon_density
double precision,	private		::	vol
double precision,	private		::	vol_a
double precision,	private		::	xtemp
double precision,	private		::	ytemp
double precision,	private		::	zlim
double precision,	private		::	ztemp

Functions

function validate_position_final(x, y, z) result(answer)

Returns validate_position( x, y, z ) if the latter is present, .TRUE. otherwise

FT 21.04.2021

Arguments

Type	Intent		Name
double precision,	intent(in)	::	x	$x$ coordinate of the desired point
double precision,	intent(in)	::	y	$y$ coordinate of the desired point
double precision,	intent(in)	::	z	$z$ coordinate of the desired point

Return Value logical

validate_position( x, y, z ) if the latter is present, .TRUE. otherwise

place_particles_lattice Module Procedure

Contents

Variables

Functions

module procedure place_particles_lattice module subroutine place_particles_lattice(this, central_density, xmin, xmax, ymin, ymax, zmin, zmax, npart_des, npart_out, stretch, thres, pos, pvol, nu, h, get_density, get_id, validate_position)

Uses

Arguments

function get_density(x, y, z) result(density)

Arguments

Return Value double precision

subroutine get_id(x, y, z, sqdetg, baryon_density, gamma_euler)

Arguments

Calls

Called by

Contents

Variables

Functions

Variables

Functions

function validate_position_final(x, y, z) result(answer)

Arguments

Return Value logical

place_particles_lattice Module Procedure

Contents

Variables

Functions

module procedure place_particles_lattice module subroutine place_particles_lattice(this, central_density, xmin, xmax, ymin, ymax, zmin, zmax, npart_des, npart_out, stretch, thres, pos, pvol, nu, h, get_density, get_id, validate_position) Interface →

Uses

Arguments

function get_density(x, y, z) result(density)

Arguments

Return Value double precision

subroutine get_id(x, y, z, sqdetg, baryon_density, gamma_euler)

Arguments

Calls

Called by

Contents

Variables

Functions

Variables

Functions

function validate_position_final(x, y, z) result(answer)

Arguments

Return Value logical

module procedure place_particles_lattice module subroutine place_particles_lattice(this, central_density, xmin, xmax, ymin, ymax, zmin, zmax, npart_des, npart_out, stretch, thres, pos, pvol, nu, h, get_density, get_id, validate_position)