Module roche_shadow 

Functions

roche_shadow

roche_shadow_py

Compute roche shadow region in equatorial plane x,y,s = shadow(q, iangle, phi, n=200, dist=5., acc=1.e-4),