evaluate_zhao_rho_hat Subroutine

public subroutine evaluate_zhao_rho_hat(p, branch, side, phi_hat, phi0_hat, phi_m_hat, n_swe_inf_hat, rho_hat)

Arguments

Type IntentOptional Attributes Name
type(zhao_params_type), intent(in) :: p
character(len=1), intent(in) :: branch
character(len=*), intent(in) :: side
real(kind=dp), intent(in) :: phi_hat
real(kind=dp), intent(in) :: phi0_hat
real(kind=dp), intent(in) :: phi_m_hat
real(kind=dp), intent(in) :: n_swe_inf_hat
real(kind=dp), intent(out) :: rho_hat

Calls

proc~~evaluate_zhao_rho_hat~~CallsGraph proc~evaluate_zhao_rho_hat evaluate_zhao_rho_hat proc~evaluate_zhao_density_hat evaluate_zhao_density_hat proc~evaluate_zhao_rho_hat->proc~evaluate_zhao_density_hat

Called by

proc~~evaluate_zhao_rho_hat~~CalledByGraph proc~evaluate_zhao_rho_hat evaluate_zhao_rho_hat proc~build_type_a_branch_from_minimum build_type_a_branch_from_minimum proc~build_type_a_branch_from_minimum->proc~evaluate_zhao_rho_hat proc~sample_monotonic_phi_hat_at_z sample_monotonic_phi_hat_at_z proc~sample_monotonic_phi_hat_at_z->proc~evaluate_zhao_rho_hat proc~sample_type_a_phi_hat_at_z sample_type_a_phi_hat_at_z proc~sample_type_a_phi_hat_at_z->proc~build_type_a_branch_from_minimum proc~sample_zhao_state_at_z sample_zhao_state_at_z proc~sample_zhao_state_at_z->proc~sample_monotonic_phi_hat_at_z proc~sample_zhao_state_at_z->proc~sample_type_a_phi_hat_at_z proc~sample_zhao_reservoir_state sample_zhao_reservoir_state proc~sample_zhao_reservoir_state->proc~sample_zhao_state_at_z proc~resolve_sheath_injection_context resolve_sheath_injection_context proc~resolve_sheath_injection_context->proc~sample_zhao_reservoir_state proc~init_particle_batch_from_config init_particle_batch_from_config proc~init_particle_batch_from_config->proc~resolve_sheath_injection_context

Source Code

  subroutine evaluate_zhao_rho_hat(p, branch, side, phi_hat, phi0_hat, phi_m_hat, n_swe_inf_hat, rho_hat)
    type(zhao_params_type), intent(in) :: p
    character(len=1), intent(in) :: branch
    character(len=*), intent(in) :: side
    real(dp), intent(in) :: phi_hat, phi0_hat, phi_m_hat, n_swe_inf_hat
    real(dp), intent(out) :: rho_hat

    real(dp) :: n_swi_hat, n_swe_f_hat, n_swe_r_hat, n_phe_f_hat, n_phe_c_hat

    call evaluate_zhao_density_hat( &
      p, branch, side, phi_hat, phi0_hat, phi_m_hat, n_swe_inf_hat, &
      n_swi_hat, n_swe_f_hat, n_swe_r_hat, n_phe_f_hat, n_phe_c_hat &
      )
    rho_hat = n_swi_hat - n_swe_f_hat - n_swe_r_hat - n_phe_f_hat - n_phe_c_hat
  end subroutine evaluate_zhao_rho_hat