Merge pull request #327 from VirtualPlanetaryLaboratory/SinewaveXUV

Sinewave xuv

Author: RoryBarnes

examples/LuminosityCycle/LuminosityCycle.png

@@ -15,7 +15,7 @@
 out = vplanet.run(path / "vpl.in", units=False)
 
 fig = plt.figure(figsize=(6.5, 8))
-plt.subplot(2, 1, 1)
+plt.subplot(3, 1, 1)
 plt.plot(
     out.sinewave.Time,
     out.sinewave.Luminosity,
@@ -28,16 +28,23 @@
     color=vplot.colors.pale_blue,
     label="None",
 )
-plt.legend(loc="best", title="Stellar Model")
+plt.legend(loc="best", title="STELLAR Model")
 plt.ylabel(r"Luminosity (L$_\odot$)")
 plt.xlabel("Time (yr)")
 
-plt.subplot(2, 1, 2)
+plt.subplot(3, 1, 2)
 plt.plot(out.sinewave.Time, out.sinewave.Temperature, color=vplot.colors.red)
 plt.plot(out.none.Time, out.none.Temperature, color=vplot.colors.pale_blue)
 plt.ylabel("Effective Temperature (K)")
 plt.xlabel("Time (yr)")
 
+plt.subplot(3, 1, 3)
+plt.plot(out.sinewave.Time, out.sinewave.LXUVTot, color=vplot.colors.red)
+plt.plot(out.none.Time, out.none.LXUVTot, color=vplot.colors.pale_blue)
+plt.ylabel(r"XUV Luminosity (L$_\odot$)")
+plt.xlabel("Time (yr)")
+
+
 # Save the figure
 ext = get_args().ext
 fig.savefig(path / f"LuminosityCycle.{ext}", bbox_inches="tight", dpi=300)

examples/LuminosityCycle/makeplot.py

@@ -6,7 +6,7 @@ dMass                     1
 dRadius                   0.00465
 dLuminosity               -1
 dLuminosityAmplitude      -0.1
-dLuminosityPeriod         -100
+dLuminosityPeriod         -10
 dLuminosityPhase          0
 
 saOutputOrder Time -Luminosity -Radius Temperature -RotPer -LXUVTot RadGyra

examples/LuminosityCycle/sinewave.in

@@ -8,10 +8,10 @@ sUnitLength               AU
 sUnitTime                 YEARS
 sUnitAngle                d
 bDoLog                    1
-iDigits                   6
+iDigits                   16
 dMinValue                 1e-10
 bDoForward                1
 bVarDt                    1
 dEta                      0.01
-dStopTime                 1000
-dOutputTime               1
+dStopTime                 100
+dOutputTime               0.1