program mlp3
      use, intrinsic :: iso_c_binding
      interface
        subroutine p(a,i,j) bind(c,name='p')
        use, intrinsic :: iso_c_binding
          real (c_float)    :: a(2,3)
          integer (c_int)   :: i, j
        end subroutine p
      end interface

      real (c_float) ::  a(2,3) ! Order 2 3
      call p(a,1,3)
      write (6,10) a(1,3)
 10   format (1x,dc,f9.1)
      stop
      end program mlp3