Constrain circle fit to have center on axis

* read-shapes/read-region-file-one-shock.py:  New option --on-axis
  forces yc=0 so that center of circle that is fitted to the bowshock
  shape has its center on the star-proplyd axis.