anchor.diff:      
     
       ! internal 
  +    integer, dimension (9) :: idrif
  +    
  +    
       ! constants
       real*8 :: EPSLON=0.000001,PMX = 999.0
       integer :: MAXROT
  
positive.code:             do il=1,nlag
                  i = (id-1)*nvarg*nlag+(iv-1)*nlag+il
                  if(np(i) == 0.) go to 6
                  rnum   = np(i)
                  gam(i) = gam(i) / dble(rnum)
                  hm(i)  = hm(i)  / dble(rnum)
                  tm(i)  = tm(i)  / dble(rnum)
  
negative.code:                 tmin, tmax, nlag,   &                            ! lag definition
                  ndir, ixd, iyd, izd,&                              ! directions and parameters
                  isill, sills, nvarg, ivtail, ivhead,ivtype,  &   ! variograms types and parameters
                  xsiz,ysiz,zsiz, &                                ! WARNING TODO dummy variable not used here, remove? 
                  np, gam, hm, tm, hv, tv)                           ! output
  
  
  