mexpvat_bowshock

from irfpy.mexpvat import domain
from irfpy.mexpvat import orbnum

if __name__ == "__main__":

    for orb in range(100, 20000):
        peri = orbnum.get_pericenter(orb)
        try:
            bsx = domain.bowshock_crossings(orb)
        except KeyboardInterrupt:
            raise
        except:
            print(orb, '{:%FT%T}'.format(peri), 'EXCEPTION')
            continue
        print(orb, '{:%FT%T}'.format(peri), end=' ')
        for bs in bsx:
            print("{:%FT%T}".format(bs), end=' ')
        print()