t_start, t_end = -20, 20
n = 1000
time = jnp.linspace(t_start, t_end, n)
flux = light_curve(system)(time).T[0]
n_plots = 5
times = jnp.linspace(t_start, t_end, n_plots) # times when the surface maps are shown
and produces the following plot where the phase of the surface maps are (more or less) aligned to the light curve:
@lgrcia: Let me know how you feel about this change!
This addresses #214 by making the changes:
and produces the following plot where the phase of the surface maps are (more or less) aligned to the light curve:
@lgrcia: Let me know how you feel about this change!