anchor.diff:        tl=zero
   #include "include/umlalid.f90"
       endif
  -#ifdef DEBUG
  -!     call wda('biflag',0d0,2,0,0,0)
  -!     if (umcalls.eq.8) then
  -!       call wda('biflag',0d0,iflag,0,0,0)
  -!       call dumpda('biflag',0,0)
  -!     call abend('biflag                                            ')
  -!     endif
  -#endif
       if(iflag.eq.1) then
   !FOX  EJF1=E0F*(ONE+DPDA)/(NUCM0/NUCMDA) ;
   !FOX  EJ1=SQRT(EJF1*EJF1+NUCMDA*NUCMDA) ;
  
positive.code:         clon(6)=dare(dpda)
        endif
      endif
  #ifdef DEBUG
  !     call wda('aclon',0d0,2,0,0,0)
  !     if (umcalls.eq.8) then
  !       call dumpda('aclon',0,0)
  !     call abend('aclon                                             ')
  !     endif
  #endif
      if(ilinc.eq.1.and.i.eq.1) then
        write(lout,10000) nd2
        if(iprint.eq.1) write(lout,10130)
  
negative.code:     ivar=ivar+2
    endif
  !--Normal Form Analysis for calculation of linear lattice functions
    if(ilinc.eq.1.or.ilinc.eq.2) then
      mfile=18
  !Eric
  