compute_and_print_quality_indicators Module Procedure

module procedure compute_and_print_quality_indicators module subroutine compute_and_print_quality_indicators(npart, pos, h, nu, nstar, path)

Uses

    • options
    • set_h
    • utility
    • RCB_tree_3D
    • sphincs_sph
    • kernel_table
  • proc~~compute_and_print_quality_indicators~~UsesGraph proc~compute_and_print_quality_indicators compute_and_print_quality_indicators RCB_tree_3D RCB_tree_3D proc~compute_and_print_quality_indicators->RCB_tree_3D kernel_table kernel_table proc~compute_and_print_quality_indicators->kernel_table module~utility utility proc~compute_and_print_quality_indicators->module~utility options options proc~compute_and_print_quality_indicators->options set_h set_h proc~compute_and_print_quality_indicators->set_h sphincs_sph sphincs_sph proc~compute_and_print_quality_indicators->sphincs_sph constants constants module~utility->constants matrix matrix module~utility->matrix

Compute the quality indicators, referring to

Daniel J. Price, Smoothed Particle Hydrodynamics and Magnetohydrodynamics. Journal of Computational Physics, 231, 3, 759-794 (2012). DOI: 10.1016/j.jcp.2010.12.011, eqs.(64), (67) and (74-75)

Rosswog, S. SPH Methods in the Modelling of Compact Objects. Living Rev Comput Astrophys 1, 1 (2015)., eqs.(6) and (9)

FT 05.10.2022


Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: npart
double precision, intent(in), DIMENSION(3,npart) :: pos
double precision, intent(in), DIMENSION(npart) :: h
double precision, intent(in), DIMENSION(npart) :: nu
double precision, intent(in), DIMENSION(npart) :: nstar
character(len=*), intent(in), optional :: path

Path to which saving the output file containing the quality indicators.


Calls

proc~~compute_and_print_quality_indicators~~CallsGraph proc~compute_and_print_quality_indicators compute_and_print_quality_indicators all_clists all_clists proc~compute_and_print_quality_indicators->all_clists exact_nei_tree_update exact_nei_tree_update proc~compute_and_print_quality_indicators->exact_nei_tree_update interp_gradw_table interp_gradw_table proc~compute_and_print_quality_indicators->interp_gradw_table interp_w_gradw_table interp_w_gradw_table proc~compute_and_print_quality_indicators->interp_w_gradw_table iorig iorig proc~compute_and_print_quality_indicators->iorig lpart lpart proc~compute_and_print_quality_indicators->lpart ncand ncand proc~compute_and_print_quality_indicators->ncand nic nic proc~compute_and_print_quality_indicators->nic rpart rpart proc~compute_and_print_quality_indicators->rpart w_no_norm w_no_norm proc~compute_and_print_quality_indicators->w_no_norm

Called by

proc~~compute_and_print_quality_indicators~~CalledByGraph proc~compute_and_print_quality_indicators compute_and_print_quality_indicators interface~compute_and_print_quality_indicators compute_and_print_quality_indicators interface~compute_and_print_quality_indicators->proc~compute_and_print_quality_indicators proc~compute_and_print_sph_variables compute_and_print_sph_variables proc~compute_and_print_sph_variables->interface~compute_and_print_quality_indicators interface~compute_and_print_sph_variables compute_and_print_sph_variables interface~compute_and_print_sph_variables->proc~compute_and_print_sph_variables

Contents


Variables

Type Visibility Attributes Name Initial
double precision, private :: Wa
double precision, private :: Wab_ha
double precision, private :: Wi
double precision, private :: Wi1
integer, private :: a
integer, private :: b
logical, private, parameter :: debug = .TRUE.
double precision, private :: dvv
double precision, private :: dx
double precision, private :: dy
double precision, private :: dz
double precision, private :: eab(3)
character(len=:), private, ALLOCATABLE :: err_msg
logical, private :: exist
double precision, private :: grW
double precision, private :: grW_ha_x
double precision, private :: grW_ha_y
double precision, private :: grW_ha_z
double precision, private :: grW_hb_x
double precision, private :: grW_hb_y
double precision, private :: grW_hb_z
double precision, private :: grWa
double precision, private :: grWb
double precision, private :: ha
double precision, private :: ha2_4
double precision, private :: ha_1
double precision, private :: ha_3
double precision, private :: ha_4
double precision, private :: hb
double precision, private :: hb2_4
double precision, private :: hb_1
double precision, private :: hb_3
double precision, private :: hb_4
integer, private :: ill
integer, private :: inde
integer, private :: index1
integer, private :: ios
integer, private :: itot
integer, private :: k
integer, private :: l
character(len=:), private, ALLOCATABLE :: namefile
double precision, private, DIMENSION(npart) :: qi_1
double precision, private, DIMENSION(3,npart) :: qi_2
double precision, private, DIMENSION(3,npart) :: qi_3
double precision, private, DIMENSION(3,3,npart) :: qi_4
double precision, private, DIMENSION(3,npart) :: qi_ham
double precision, private :: rab
double precision, private :: rab2
double precision, private :: rab_1
integer, private :: unit_qi
double precision, private :: va
double precision, private :: vb
double precision, private :: vol_b
double precision, private :: xa
double precision, private :: xb
double precision, private :: ya
double precision, private :: yb
double precision, private :: za
double precision, private :: zb