program score
  implicit none

  integer, dimension(3,4) :: scores
  logical, dimension(3,4) :: above
  integer :: average

  print *, "score.f90"
  scores = reshape((/ 85, 71, 66,   76, 45, 45,                  &   
                      90, 50, 21,   60, 80, 55 /), shape(scores))
  print *, "scores", scores
  print *, "maximum", maxval(scores)
  print *, "stu max", maxval(scores, dim=2)
  print *, "max loc", maxloc(maxval(scores, dim=2))
  average = sum(scores) / size(scores)
  print *, "average", average
  above = scores > average
  print *, "above", above
  print *, "count", count(above)
  print *, "any(all", any(all(above, dim=2))
  print *, "any(all", any(all(above, dim=1))
end program score