anchor.diff:  
     implicit none
   
  -  integer i,mf1,mf2,mf3,mf4,mf5,mfile,nd2,ndim,ndpt,nis,no,nv,damap,a1,a1i,a2,a2i,f,fc,fs,rot,xy,h,hc,hs,h4,df,bb1,bb2,haux
  +  integer i,mf1,mf2,mf3,mf4,mf5,mfile,nd2,ndim,ndpt,nis,no,nv,damap&
  +       &,a1,a1i,a2,a2i,f,fc,fs,rot,xy,h,hc,hs,h4,df,bb1,bb2,haux
     real(kind=fPrec) angle,coe,radn,x2pi
     dimension damap(6),a1(6),a1i(6),a2(6),a2i(6)
     dimension rot(6),xy(6),df(6)
  
positive.code:   rewind mf3
    rewind mf4
    rewind mf5
    call daread(damap,nd2,mfile,zero)
  
    ! Normal Form Analysis
    call mapnorm(damap,f,a2,a1,xy,h,nord1)
    call dainv(a1,nv,a1i,nv)
    call dainv(a2,nv,a2i,nv)
    call ctor(f,fc,fs)
    call gettura(angle,radn)
    call taked(xy,1,rot)
  
negative.code:   write(lout,*) (angle(i),i=1,ndim)
  
    ! Clean-Up
    call dadal(damap,nv)
    call dadal(a1,nv)
    call dadal(a1i,nv)
    call dadal(a2,nv)
    call dadal(a2i,nv)
    call dadal(f,1)
    call dadal(fc,1)
    call dadal(fs,1)
    call dadal(rot,nv)
    call dadal(xy,nv)
    call dadal(h,1)
    call dadal(h4,1)
    call dadal(hc,1)
    call dadal(hs,1)
    call dadal(df,nv)
    call dadal(bb1,1)
    call dadal(bb2,1)
    call dadal(haux,1)
  