anchor.diff:  ! ================================
   ! = gamma(t) can be changed here =
   ! ================================
  -            dd_here=Gt(i,j)+(ddx(i)+ddy(j))*0.5D0
  +            dd_here=Gt(i,j)
               
               if (dd_here.GE. 0.01) then
                   far = dabs((dd_here+2.0D0)**2-0.25D0)*10.0D0
  
positive.code: ! ================================
  ! = gamma(t) can be changed here =
  ! ================================
              dd_here=Gt(i,j)+(ddx(i)+ddy(j))*0.5D0
              
              if (dd_here.GE. 0.01) then
                  far = dabs((dd_here+2.0D0)**2-0.25D0)*10.0D0
  
negative.code: 
  c
        SUBROUTINE nsst
       *(C,Gt,origin_val,hx,hy,
       * Bk,cmin,cmax,nx,ny,symm)
  
  cf2py threadsafe
  