Added LuminosityCycle figure.

VirtualPlanetaryLaboratory · Feb 14, 2024 · 5d83759 · 5d83759
1 parent 4b89543
commit 5d83759
Show file tree

Hide file tree

Showing 5 changed files with 46 additions and 3 deletions.
diff --git a/examples/LuminosityCycle/LuminosityCycle.pdf b/examples/LuminosityCycle/LuminosityCycle.pdf
diff --git a/examples/LuminosityCycle/makeplot.py b/examples/LuminosityCycle/makeplot.py
@@ -0,0 +1,43 @@
+import pathlib
+import sys
+
+import matplotlib.pyplot as plt
+import vplot
+
+import vplanet
+
+# Path hacks
+path = pathlib.Path(__file__).parents[0].absolute()
+sys.path.insert(1, str(path.parents[0]))
+from get_args import get_args
+
+# Run vplanet
+out = vplanet.run(path / "vpl.in", units=False)
+
+fig = plt.figure(figsize=(6.5, 8))
+plt.subplot(2, 1, 1)
+plt.plot(
+    out.sinewave.Time,
+    out.sinewave.Luminosity,
+    color=vplot.colors.red,
+    label="SineWave",
+)
+plt.plot(
+    out.none.Time,
+    out.none.Luminosity,
+    color=vplot.colors.pale_blue,
+    label="None",
+)
+plt.legend(loc="best", title="Stellar Model")
+plt.ylabel(r"Luminosity (L$_\odot$)")
+plt.xlabel("Time (yr)")
+
+plt.subplot(2, 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)")
+
+# Save the figure
+ext = get_args().ext
+fig.savefig(path / f"LuminosityCycle.{ext}", bbox_inches="tight", dpi=600)
diff --git a/examples/LuminosityCycle/none.in b/examples/LuminosityCycle/none.in
@@ -3,6 +3,6 @@ sName	                  none
 saModules	              stellar
 sStellarModel             none
 dMass                     1
-dRadius                   0.00454
+dRadius                   0.00465
 dLuminosity               -1
 saOutputOrder Time -Luminosity -Radius Temperature -RotPer -LXUVTot RadGyra
diff --git a/examples/LuminosityCycle/sinewave.in b/examples/LuminosityCycle/sinewave.in
@@ -3,7 +3,7 @@ sName	                  sinewave
 saModules	              stellar
 sStellarModel             sinewave
 dMass                     1
-dRadius                   0.00454
+dRadius                   0.00465
 dLuminosity               -1
 dLuminosityAmplitude      -0.1
 dLuminosityPeriod         -100

diff --git a/examples/RadHeat/makeplot.py b/examples/RadHeat/makeplot.py
@@ -16,7 +16,7 @@
 # Run vplanet
 out = vplanet.run(path / "vpl.in", units=False)
 
-fig = plt.figure(figsize=(8.5, 8))
+fig = plt.figure(figsize=(6.5, 8))
 plt.subplot(2, 2, 1)
 plt.plot(
     out.earth.Time,