compute_z_high_box_potential_statistics Subroutine

public subroutine compute_z_high_box_potential_statistics(mesh, sim, snapshot, phi_mean, phi_std, phi_min, phi_max)

box の z-high 全断面における電位統計をセル中心 N x N 格子で評価する。

N は reservoir face 平均と同じ sim.injection_face_phi_grid_n を使う。 注入speciesの開口には依存せず、周期セル全断面を基準面として評価する。

Arguments

Type IntentOptional Attributes Name
type(mesh_type), intent(in) :: mesh
type(sim_config), intent(in) :: sim
type(electrostatic_snapshot_type), intent(inout) :: snapshot
real(kind=dp), intent(out) :: phi_mean
real(kind=dp), intent(out) :: phi_std
real(kind=dp), intent(out) :: phi_min
real(kind=dp), intent(out) :: phi_max

Calls

proc~~compute_z_high_box_potential_statistics~~CallsGraph proc~compute_z_high_box_potential_statistics compute_z_high_box_potential_statistics none~eval_local_phi electrostatic_snapshot_type%eval_local_phi proc~compute_z_high_box_potential_statistics->none~eval_local_phi proc~resolve_face_sampling_geometry resolve_face_sampling_geometry proc~compute_z_high_box_potential_statistics->proc~resolve_face_sampling_geometry proc~lower_ascii lower_ascii proc~resolve_face_sampling_geometry->proc~lower_ascii

Source Code

  subroutine compute_z_high_box_potential_statistics( &
    mesh, sim, snapshot, phi_mean, phi_std, phi_min, phi_max &
    )
    type(mesh_type), intent(in) :: mesh
    type(sim_config), intent(in) :: sim
    type(electrostatic_snapshot_type), intent(inout) :: snapshot
    real(dp), intent(out) :: phi_mean, phi_std, phi_min, phi_max

    if (.not. sim%use_box) then
      error stop 'z-high box potential statistics require sim.use_box=true.'
    end if
    call compute_rectangular_face_potential_statistics( &
      mesh, sim, snapshot, 'z_high', sim%box_min, sim%box_max, &
      sim%injection_face_phi_grid_n, phi_mean, phi_std, phi_min, phi_max &
      )
  end subroutine compute_z_high_box_potential_statistics