anchor.diff:                  
           grad = (aa - 1d0)/(1d0 - t1) * t1 * z(i) **(cc -1d0)
        & * cc / sigma
  -        grad = grad + -  z(i) **(cc -1d0) * cc/sigma -
  +        grad = grad + (-  z(i) **(cc -1d0) * cc/sigma) -
        & (cc -1d0)/(z(i) *sigma)
           grad = -grad
           if (nloc .NE. 1) then 
  
positive.code:         
          gradlike(i) = (aa - 1d0)/(1d0 - t1) * t1 * 
       & z(i) **(cc -1d0) * cc / sigma
          gradlike(i) = gradlike(i) + -  z(i) **(cc -1d0)
       & * cc/sigma - (cc -1d0)/(z(i) *sigma)
        enddo
        END SUBROUTINE 
        
  
negative.code:       
  	      DATA S /1.0e-5/, C /8.5/ S3 /8.333333333e-2/
  	      DATA S4 /8.333333333e-3/, S5 /3.968253968e-3/
  	      DATA D1 /-0.5772156649/
  	      
  	      psi = 0
  	
  