subroutine test_sub(a,b)
      integer :: a
      double precision :: b(3,2)
      a = 3
      b(1,1) = 1.0
      b(1,2) = 2.0
      b(2,1) = 3.0
      b(2,2) = 4.0
      b(3,1) = 5.0
      b(3,2) = 6.0
      end subroutine