# 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)
