anchor.diff:  ! ================================
   ! = gamma(t) can be changed here =
   ! ================================
  -            dd_here=Gt(i,j)            
  +            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
                 C(i,j)=origin_val / dd_here
  
positive.code: ! ================================
              dd_here=Gt(i,j)
              
  
  
              if (C(i,j) .EQ. 0.0D0) then
  
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
  