ice_gather_flux_y_layer_impl Subroutine

private pure subroutine ice_gather_flux_y_layer_impl(vh, val4, layer, tr_flux_y_work, ncat, nk, nx, ny)

Meridional twin of ice_gather_flux_x_layer_impl.

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: vh(nx,ny+1,ncat)
real(kind=wp), intent(in) :: val4(nx,ny,ncat,nk)
integer, intent(in) :: layer
real(kind=wp), intent(out) :: tr_flux_y_work(nx,ny+1,ncat)
integer, intent(in) :: ncat
integer, intent(in) :: nk
integer, intent(in) :: nx
integer, intent(in) :: ny

Called by

proc~~ice_gather_flux_y_layer_impl~~CalledByGraph proc~ice_gather_flux_y_layer_impl ice_gather_flux_y_layer_impl proc~ice_pass_y ice_pass_y proc~ice_pass_y->proc~ice_gather_flux_y_layer_impl proc~ice_transport_step ice_transport_step proc~ice_transport_step->proc~ice_pass_y proc~engine_step_ice engine_step_ice proc~engine_step_ice->proc~ice_transport_step proc~driver_run_ocean driver_run_ocean proc~driver_run_ocean->proc~engine_step_ice proc~rdb_ocean_step rdb_ocean_step proc~rdb_ocean_step->proc~engine_step_ice proc~driver_run driver_run proc~driver_run->proc~driver_run_ocean

Variables

Type Visibility Attributes Name Initial
integer, private :: c
integer, private :: i
integer, private :: j

Source Code

   pure subroutine ice_gather_flux_y_layer_impl(vh, val4, layer, tr_flux_y_work, ncat, nk, nx, ny)
      !! Meridional twin of `ice_gather_flux_x_layer_impl`.
      integer, intent(in) :: layer, ncat, nk, nx, ny
      real(wp), intent(in) :: vh(nx, ny + 1, ncat)
      real(wp), intent(in) :: val4(nx, ny, ncat, nk)
      real(wp), intent(out) :: tr_flux_y_work(nx, ny + 1, ncat)
      integer :: i, j, c

      do concurrent(j=1:ny + 1, i=1:nx, c=1:ncat)
         if (j == 1) then
            tr_flux_y_work(i, j, c) = val4(i, 1, c, layer)
         else if (j == ny + 1) then
            tr_flux_y_work(i, j, c) = val4(i, ny, c, layer)
         else if (vh(i, j, c) >= 0.0_wp) then
            tr_flux_y_work(i, j, c) = val4(i, j - 1, c, layer)
         else
            tr_flux_y_work(i, j, c) = val4(i, j, c, layer)
         end if
      end do
   end subroutine ice_gather_flux_y_layer_impl