subroutine surfacl (v, rij)
      real*8 v, rij
      call surface(v, rij)
      return
      end