Program Main
  Implicit None
  Integer :: a(2,3), b(3), m
  Integer :: i,j

  a = Reshape( (/1,2,3,4,5,6/), (/2,3/) )
  print *, '------------------------------------------'
  print *, 'a = '
  write(*,*) (a(1,j),j=1,3)
  write(*,*) (a(2,j),j=1,3)
  print *, '------------------------------------------'
  Print *, 'Maxloc(a) = ', Maxloc(a)
  Print *, 'Maxloc(a)(1) = ', Maxloc(Maxval(a,2))
  Print *, 'Maxloc(a)(2) = ', Maxloc(Maxval(a,1))
  Print *, 'Maxloc(a, 1) = ', Maxloc(a,1)
  Print *, 'Maxloc(a, 2) = ', Maxloc(a,2)
  Print *, 'Maxloc(a, Mask = a < 3) = ', Maxloc(a, Mask = a < 3)
  print *, '------------------------------------------'
  Print *, 'Maxval(a) = ', Maxval(a)
  Print *, 'Maxval(a, 1) = ', Maxval(a,1)
  Print *, 'Maxval(a, 2) = ', Maxval(a,2)
  Print *, 'Maxval(a, a > 3) = ', Maxval(a, a>3)
  print *, '------------------------------------------'
  Print *, 'Size(a) = ', Size(a)
  Print *, 'Size(a,1) = ', Size(a,1)
  Print *, 'Size(a,2) = ', Size(a,2)
End Program