sphere = TriangulatedSphere(model=basedMDwarfModel, times=times, nTheta=5, nPhi=5, noiseScale=0.001)
sphere.evolve(0.05, 1000*60, 0.25, 30, 5*60, 1, 10)
#sphere.evolve(0.05, 100*60, 0.01, 30, 5*60, 0.1, 10)
#flareProb, flareMeanDuration, flareMeanAmplitude, flareMeanLat, flareDurationScale, flareAmplitudeScale, flareLatScale):

# Assuming that the components are tidally locked we only need to subtend one hemisphere to get the incident luminosity at any given time
hemiLuminosity = sphere.get_luminosity_over_solid_angle(0, 0, omega=2*np.pi, time='*')