int4_2D_or_int4 Function

        function int4_2D_or_int4(a, b) result(r)
          !< Bitwise 'OR' over one 2D integer array and integer of kind 4
            implicit none
            integer(kind=4), dimension(:, :), intent(in) :: a
            integer(kind=4), intent(in) :: b
            integer(kind=4), dimension(1:2) :: n
            integer(kind=4), dimension(:, :), allocatable :: r
            integer :: i, j
            n = shape(a)
            allocate(r(1:n(1), 1:n(2)))
            do j = 1, n(2)
             do i = 1, n(1)
                r(i,j) = ior(a(i, j), b)
             end do
            end do
        end function int4_2D_or_int4