c****************************************************************
c
      subroutine init_phi0(n)
c Initial constant for Euler angles.
c
c          phi0 = (idlk-dlk)*pi2/n
c          phii = idlk*pi2/n

      implicit double precision(a-h,o-z)
      parameter(pi2=2.d0*3.1415926535d0)
       common/ischeme/icircular,iele,ihydro,iclose,iknot
     +,neuler
       common/tor_param/dr,cg,phi0,dlk,phii,subseg

         write(*,*) 'dlk=',dlk
         idlk = dlk
         if (abs(dlk-idlk) .gt. 0.5d0) then
            if (dlk .gt. 0) then
                idlk = idlk + 1
            else
                idlk = idlk-1
            endif
         endif
         phi0 = (idlk-dlk)*pi2/n
         phii = idlk*pi2/n

14    format('phi0 = ',f9.5,4x,'dlk = ',f9.5,4x,'phii = ',f9.5)

      return
      end