use inpolygin for private flux, fix R_Thomson=0.99