anchor.diff:  ! ================================
               dd_here=Gt(i,j)
               
  +            if (dd_here.GE. 0.01) then
  +                far = dabs((dd_here+2.0D0)**2-0.25D0)*10.0D0
  +            else
  +                far = infinity
  +            end if
   
   
               if (C(i,j) .EQ. 0.0D0) then
  
positive.code:               fl = INT(dd_here)
                rem = dd_here - fl
                N = fl
                              
                C(i,j) = C(i,j) * snu
                CALL RKBESL(C(i,j),rem,fl+1,1,BK,N)              
                C(i,j)=prefac * (C(i,j) ** dd_here) * BK(fl+1)
                              
  
              endif
  
negative.code: 
        DOUBLE PRECISION C(nx,ny), Gt(nx,ny)
        DOUBLE PRECISION origin_val
        DOUBLE PRECISION rem, dd_here, far
        DOUBLE PRECISION GA, prefac, snu
        INTEGER nx, ny, i, j, fl, N, nd, cmin, cmax
        DOUBLE PRECISION BK(50), DGAMMA
  