zero_array_real_array_4 Subroutine

  subroutine zero_array_real_array_4(array_in)
    implicit none
    real, dimension(:, :, :, :), intent(in out) :: array_in
    integer :: counter, start, end

    start = lbound(array_in, dim = 4)
    end = ubound(array_in, dim = 4)

    !$OMP PARALLEL DO DEFAULT(none) &
    !$OMP SHARED(start, end, array_in) &
    !$OMP SCHEDULE(static)
    do counter = start, end
       array_in(:, :, :, counter) = 0.0
    end do
    !$OMP END PARALLEL DO
  end subroutine zero_array_real_array_4