anchor.diff:  !FOX  CORRNEW(8)=SMIDA(2)-DPDAV ;
       endif
       call dacct(corrau1,nvar,corrnew,nvar,corrau2,nvar)
  -#ifdef DEBUG
  -!     call wda('b436   ',0d0,2,0,0,0)
  -!     if (umcalls.eq.8) then
  -!       call dumpda('b436   ',0,0)
  -!     call abend('b436                                              ')
  -!     endif
  -#endif
       do 436 kkk=1,6
   !FOX  CORRAU2(KKK)=CORRAU2(KKK)+DPDAV2(KKK) ;
   436   continue
  
positive.code:     corr(3,2)=coefh1/det1
    endif
  470  continue
  #ifdef DEBUG
  !     call wda('a470   ',0d0,2,0,0,0)
  !     if (umcalls.eq.8) then
  !       call dumpda('a470   ',0,0)
  !     call abend('a470                                              ')
  !     endif
  #endif
  
    call dadal(damap,6)
    call dadal(damapi,6)
  
negative.code:       aml6(i,i)=aml6(i,i)-one
      enddo
    endif
    call dacop(yp(1),damap(2))
    call dacop(yp(2),damap(4))
    if(ndimf.eq.3) then
  