anchor.diff:  cf2py double precision intent(in),check(origin_val>0)::origin_val
   
         DOUBLE PRECISION C(nx,ny), Gt(nx,ny)
  -      DOUBLE PRECISION ddx(nx), hx(nx), ddy(ny), hy(ny)
  +      DOUBLE PRECISION hx(nx), hy(ny)
         DOUBLE PRECISION origin_val
         DOUBLE PRECISION rem, dd_here, far
         DOUBLE PRECISION GA, prefac, snu
  
positive.code: 
  c
        SUBROUTINE nsst
       *(C,Gt,origin_val,ddx,ddy,hx,hy,
       * Bk,cmin,cmax,nx,ny,symm)
  
  cf2py threadsafe
  
negative.code: ! ================================
  ! = gamma(t) can be changed here =
  ! ================================
              dd_here=Gt(i,j)
              
              if (dd_here.GE. 0.01) then
                  far = dabs((dd_here+2.0D0)**2-0.25D0)*10.0D0
  