subroutine sub1(n)
  real, intent(in) :: n