# Set other parameters ShowPlanetMoonEclipses = True # True: the reality; False would be no mutual # eclipses. Of course unphysical, but useful for tests and comparisons) ShowPlanet = False # True: Planet+Moon; False: Moon only Noise = 0 # [ppm per minute]; 0 = no noise is added NumberOfTransits = 0 # How many (randomly chosen) transits are observed; # if 0 then all available are sampled (and their number is 10*Quality). PhaseToHighlight = 0.2 # If no highlighting is desired, choose value < 0 Quality = 250 # Radius of star in pixels --> size of numerical sampling grid NumberOfSamples = 250 # How many transits are to be sampled # Curve MyNewCurve = PyOSE.curve(StellarRadius, limb1, limb2, PlanetRadius, PlanetAxis, PlanetImpact, PlanetPeriod, MoonRadius, MoonAxis, MoonEccentricity, MoonAscendingNode, MoonLongitudePeriastron, MoonInclination, ShowPlanetMoonEclipses, ShowPlanet, Quality, NumberOfSamples, Noise, NumberOfTransits) Time = MyNewCurve[0][1:] Flux = MyNewCurve[1][1:] ax = plt.axes() plt.plot(Time, Flux, color='k') plt.rc('text', usetex=True) plt.rc('font', family='serif') plt.tick_params(axis='both', which='major', labelsize=16) plt.xlabel('time around planetary mid-transit [days]', fontsize=16) plt.ylabel('normalized stellar brightness [ppm]', fontsize=16) plt.axis([-0.2, +0.2, -100, 1], set_aspect='equal', fontsize=16) ax.tick_params(direction='out') MoonAxis = 61162
# eclipses. Of course unphysical, but useful for tests and comparisons) ShowPlanet = False # True: Planet+Moon; False: Moon only Noise = 0 # [ppm per minute]; 0 = no noise is added NumberOfTransits = 0 # How many (randomly chosen) transits are observed; # if 0 then all available are sampled (and their number is 10*Quality). PhaseToHighlight = 0.2 # If no highlighting is desired, choose value < 0 Quality = 250 # Radius of star in pixels --> size of numerical sampling grid NumberOfSamples = 250 # How many transits are to be sampled # Curve PlanetImpact = 0.0# 0.25 # [0..1.x]; central transit is 0. MoonInclination = 90.0 # 0..90 in degrees. 0 is the reference plain (no incl). MyNewCurve = PyOSE.curve(StellarRadius, limb1, limb2, PlanetRadius, PlanetAxis, PlanetImpact, PlanetPeriod, MoonRadius, MoonAxis, MoonEccentricity, MoonAscendingNode, MoonLongitudePeriastron, MoonInclination, ShowPlanetMoonEclipses, ShowPlanet, Quality, NumberOfSamples, Noise, NumberOfTransits) Time = MyNewCurve[0][1:] Flux = MyNewCurve[1][1:] ax = plt.axes() plt.plot(Time, Flux, color = 'k') plt.rc('text', usetex=True) plt.rc('font', family='serif') plt.tick_params(axis='both', which='major', labelsize=16) plt.xlabel('time around planetary mid-transit [days]',fontsize=16) plt.ylabel('normalized stellar brightness [ppm]',fontsize=16) plt.axis([-0.2, +0.2, -100, 1], set_aspect='equal', fontsize=16) ax.tick_params(direction='out') """
ShowPlanetMoonEclipses = True # True: the reality; False would be no mutual # eclipses. Of course unphysical, but useful for tests and comparisons) ShowPlanet = False # True: Planet+Moon; False: Moon only Noise = 0 # [ppm per minute]; 0 = no noise is added NumberOfTransits = 0 # How many (randomly chosen) transits are observed; # if 0 then all available are sampled (and their number is 10*Quality). PhaseToHighlight = 0.09 # If no highlighting is desired, choose value < 0 Quality = 250 # Radius of star in pixels --> size of numerical sampling grid NumberOfSamples = 250 # How many transits are to be sampled # 3D model plt.rc('text', usetex=True) plt.rc('font', family='serif') MyModelview = PyOSE.modelview( StellarRadius, limb1, limb2, PlanetRadius, PlanetImpact, MoonRadius, MoonAxis, MoonEccentricity, MoonAscendingNode, MoonLongitudePeriastron, MoonInclination, PhaseToHighlight, Quality) ax = plt.axes() ax.arrow(0, 0.4, 0.364, 0, head_width=0.0, head_length=0.0, fc='k', ec='k', zorder = 10) # direction #ax.arrow(0.364, 0.4, 0, -0.1, head_width=0.0, head_length=0.0, fc='k', ec='k') ax.arrow(0.356, 0.0, 0.0, 0.205, head_width=0.0, head_length=0.0, fc='k', ec='k', zorder = 10) # b_S ax.tick_params(direction='out') plt.tick_params(axis='both', which='major', labelsize=16) plt.annotate(r"direction", xy=(0.17, 0.42), size=16) plt.annotate(r"$b_S$", xy=(0.4, 0.05), size=16) #plt.annotate(r"$b_P$ ", xy=(-0.12, 0.1), size=16)
# Set other parameters ShowPlanetMoonEclipses = True # True: the reality; False would be no mutual # eclipses. Of course unphysical, but useful for tests and comparisons ShowPlanet = False # True: Planet+Moon; False: Moon only Noise = 0 # [ppm per minute]; 0 = no noise is added NumberOfTransits = 0 # How many (randomly chosen) transits are observed; # if 0 then all available are sampled (and their number is 10*Quality). PhaseToHighlight = 0.2 # If no highlighting is desired, choose value < 0 Quality = 5000 # Radius of star in pixels --> size of numerical sampling grid NumberOfSamples = 250 # How many transits are to be sampled # River MyRiverKepler = PyOSE.river(StellarRadius, limb1, limb2, PlanetRadius, PlanetAxis, PlanetImpact, PlanetPeriod, MoonRadius, MoonAxis, MoonEccentricity, MoonAscendingNode, MoonLongitudePeriastron, MoonInclination, ShowPlanetMoonEclipses, Quality, NumberOfSamples, Noise) # River function returns pixel map. To plot time axis, call function timeaxis MyTime = PyOSE.timeaxis(PlanetPeriod, PlanetAxis, MoonRadius, StellarRadius, Quality) plt.imshow(MyRiverKepler, cmap=cm.gray, interpolation='none', aspect='auto', extent=[MyTime[0], -MyTime[0], 1, 0]) plt.rc('text', usetex=True) plt.rc('font', family='serif') ax = plt.axes()
ShowPlanetMoonEclipses = True # True: the reality; False would be no mutual # eclipses. Of course unphysical, but useful for tests and comparisons ShowPlanet = False # True: Planet+Moon; False: Moon only Noise = 0 # [ppm per minute]; 0 = no noise is added NumberOfTransits = 0 # How many (randomly chosen) transits are observed; # if 0 then all available are sampled (and their number is 10*Quality). PhaseToHighlight = 0.2 # If no highlighting is desired, choose value < 0 Quality = 5000 # Radius of star in pixels --> size of numerical sampling grid NumberOfSamples = 250 # How many transits are to be sampled # River MyRiverKepler = PyOSE.river( StellarRadius, limb1, limb2, PlanetRadius, PlanetAxis, PlanetImpact, PlanetPeriod, MoonRadius, MoonAxis, MoonEccentricity, MoonAscendingNode, MoonLongitudePeriastron, MoonInclination, ShowPlanetMoonEclipses, Quality, NumberOfSamples, Noise) # River function returns pixel map. To plot time axis, call function timeaxis MyTime = PyOSE.timeaxis( PlanetPeriod, PlanetAxis, MoonRadius, StellarRadius, Quality) plt.imshow(MyRiverKepler, cmap=cm.gray, interpolation='none', aspect='auto', extent=[MyTime[0], -MyTime[0], 1, 0]) plt.rc('text', usetex=True) plt.rc('font', family='serif') ax = plt.axes() T_dur_P = 0.1477 / 2 plt.plot([T_dur_P, T_dur_P], [0, 1], 'k--', linewidth = 1.5)
ShowPlanetMoonEclipses = True # True: the reality; False would be no mutual # eclipses. Of course unphysical, but useful for tests and comparisons ShowPlanet = False # True: Planet+Moon; False: Moon only Noise = 0 # [ppm per minute]; 0 = no noise is added NumberOfTransits = 0 # How many (randomly chosen) transits are observed; # if 0 then all available are sampled (and their number is 10*Quality). PhaseToHighlight = 0.2 # If no highlighting is desired, choose value < 0 Quality = 250 # Radius of star in pixels --> size of numerical sampling grid NumberOfSamples = 250 # How many transits are to be sampled # 3D model plt.rc('text', usetex=True) plt.rc('font', family='serif') MyModelview = PyOSE.modelview( StellarRadius, limb1, limb2, PlanetRadius, PlanetImpact, MoonRadius, MoonAxis, MoonEccentricity, MoonAscendingNode, MoonLongitudePeriastron, MoonInclination, PhaseToHighlight, Quality) ax = plt.axes() ax.arrow(0, 0, 0, 0.14, width=0.001, head_width=0.03, head_length=0.03, fc='k', ec='k', zorder = 10) ax.tick_params(direction='out') plt.tick_params(axis='both', which='major', labelsize=16) plt.annotate(r"$b_P$ ", xy=(0.04, 0.02), size=16) plt.xlabel('distance [stellar radii]',fontsize=16) plt.ylabel('distance [stellar radii]',fontsize=16) ax.set_aspect('equal') plt.savefig("fig_7a.eps", bbox_inches='tight') plt.savefig("fig_7a.pdf", bbox_inches='tight') MyModelview.show()
ShowPlanetMoonEclipses = True # True: the reality; False would be no mutual # eclipses. Of course unphysical, but useful for tests and comparisons ShowPlanet = False # True: Planet+Moon; False: Moon only Noise = 0 # [ppm per minute]; 0 = no noise is added NumberOfTransits = 0 # How many (randomly chosen) transits are observed; # if 0 then all available are sampled (and their number is 10*Quality). PhaseToHighlight = 0.2 # If no highlighting is desired, choose value < 0 Quality = 250 # Radius of star in pixels --> size of numerical sampling grid NumberOfSamples = 250 # How many transits are to be sampled # 3D model plt.rc('text', usetex=True) plt.rc('font', family='serif') MyModelview = PyOSE.modelview(StellarRadius, limb1, limb2, PlanetRadius, PlanetImpact, MoonRadius, MoonAxis, MoonEccentricity, MoonAscendingNode, MoonLongitudePeriastron, MoonInclination, PhaseToHighlight, Quality) ax = plt.axes() ax.arrow(0, 0, 0, 0.14, width=0.001, head_width=0.03, head_length=0.03, fc='k', ec='k', zorder=10) ax.tick_params(direction='out') plt.tick_params(axis='both', which='major', labelsize=16)