cavity_status_impl Subroutine

private pure subroutine cavity_status_impl(status, active, buf)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: status(:,:)
real(kind=wp), intent(in) :: active(:,:)
real(kind=wp), intent(inout) :: buf(:,:,:)

Called by

proc~~cavity_status_impl~~CalledByGraph proc~cavity_status_impl cavity_status_impl proc~fill_cavity_melt_status fill_cavity_melt_status proc~fill_cavity_melt_status->proc~cavity_status_impl

Variables

Type Visibility Attributes Name Initial
integer, private :: i
integer, private :: j
integer, private :: nx
integer, private :: ny
real(kind=wp), private :: qnan

Source Code

   pure subroutine cavity_status_impl(status, active, buf)
      ! assumed-shape-ok: diag fill — fires once per output frame (cadence-bounded).
      integer, intent(in)     :: status(:, :)  ! assumed-shape-ok: diag fill — cadence-bounded
      real(wp), intent(in)    :: active(:, :)  ! assumed-shape-ok: diag fill — cadence-bounded
      real(wp), intent(inout) :: buf(:, :, :)  ! assumed-shape-ok: diag fill — cadence-bounded
      integer :: i, j, nx, ny
      real(wp) :: qnan
      nx = min(size(buf, 1), size(status, 1), size(active, 1))
      ny = min(size(buf, 2), size(status, 2), size(active, 2))
      qnan = ieee_value(0.0_wp, ieee_quiet_nan)
      do concurrent(j=1:ny, i=1:nx)
         if (active(i, j) > 0.5_wp) then
            buf(i, j, 1) = real(status(i, j), wp)
         else
            buf(i, j, 1) = qnan
         end if
      end do
   end subroutine cavity_status_impl