anchor.diff:    rewind mf3
     rewind mf4
     rewind mf5
  -  call daread(damap,nd2,mfile,zero)
  +  call daread(damap,nv,mfile,zero)
   
     ! Normal Form Analysis
  +  !call mapnorm(damap,f,a2,a1,xy,h,nord1)
     call mapnorm(damap,f,a2,a1,xy,h,nord1)
  -  call dainv(a1,nv,a1i,nv)
  -  call dainv(a2,nv,a2i,nv)
  +  call daprid(a1,1,nd2,66)
  +  call daprid(a2,1,nd2,66)
  +  call etinv(a1,a1i)
  +  call etinv(a2,a2i)
     call ctor(f,fc,fs)
     call gettura(angle,radn)
     call taked(xy,1,rot)
  
positive.code:   write(lout,*) (angle(i),i=1,ndim)
  
    ! Clean-Up
    call dadal(damap,nd2)
    call dadal(a1,nd2)
    call dadal(a1i,nd2)
    call dadal(a2,nd2)
    call dadal(a2i,nd2)
    call dadal(f,1)
    call dadal(fc,1)
    call dadal(fs,1)
    call dadal(rot,nd2)
    call dadal(xy,nd2)
    call dadal(h,1)
    call dadal(h4,1)
    call dadal(hc,1)
    call dadal(hs,1)
    call dadal(df,nd2)
    call dadal(bb1,1)
    call dadal(bb2,1)
    call dadal(haux,1)
  
negative.code:   call idprset(-102)
    call mld_allocArrays(.true.)
    call lieinit(no,nv,ndim,ndpt,0,nis)
    call etall(damap,nv)
    call etall(xy,nv)
    call etall(bb1,1)
    call etall(bb2,1)
    call etall(haux,1)
    call etallnom(a1,nv,'A1        ')
    call etallnom(a1i,nv,'A1I       ')
    call etallnom(a2,nv,'A2        ')
    call etallnom(a2i,nv,'A2I       ')
    call etallnom(f,1,'F         ')
    call etallnom(fc,1,'FC        ')
    call etallnom(fs,1,'FS        ')
    call etallnom(rot,nv,'ROT       ')
    call etallnom(h,1,'H         ')
    call etallnom(h4,1,'H4        ')
    call etallnom(hc,1,'HC        ')
    call etallnom(hs,1,'HS        ')
    call etallnom(df,nv,'DF        ')
    rewind mfile
    rewind 26
    rewind mf1
  