Wrapper SUBROUTINE to get the value of the Hamiltonian constraint
computed using using the hydro ID mapped from the particles to the refined
mesh, at a given mesh point (i,j,k)
, on the refinement level l
, for the
tpo object tpof
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
class(tpo), | intent(inout) | :: | tpof |
tpo object to use |
||
integer, | intent(in) | :: | i |
index of the desired point |
||
integer, | intent(in) | :: | j |
index of the desired point |
||
integer, | intent(in) | :: | k |
index of the desired point |
||
integer, | intent(in) | :: | l |
Index of the refinement level |
||
double precision, | intent(out) | :: | hc |
Value of the Hamiltonian constraint at |
SUBROUTINE get_ham_parts( tpof, i, j, k, l, hc )
!# Wrapper SUBROUTINE to get the value of the Hamiltonian constraint
! computed using using the hydro ID mapped from the particles to the refined
! mesh, at a given mesh point `(i,j,k)`, on the refinement level `l`, for the
! [[tpo]] object `tpof`
IMPLICIT NONE
CLASS(tpo), INTENT(INOUT):: tpof
!! [[tpo]] object to use
INTEGER, INTENT(IN):: i
!! \(x\) index of the desired point
INTEGER, INTENT(IN):: j
!! \(y\) index of the desired point
INTEGER, INTENT(IN):: k
!! \(z\) index of the desired point
INTEGER, INTENT(IN):: l
!! Index of the refinement level
DOUBLE PRECISION, INTENT(OUT):: hc
!! Value of the Hamiltonian constraint at \((i,j,k)\)
hc= tpof% get_HC_parts( i, j, k, l )
END SUBROUTINE get_ham_parts