changes to manual plots
@@ -154,25 +154,31 @@ plt.close()
out2 = vplot.GetOutput('womoon')
# plt.rcParams["mathtext."]
+tm, oblm, pAm, precf = np.loadtxt('marshnb/solarsys.Mars.backward',unpack = True)
+tm1, em1, incm1, argpm1, longam1, oblm1, pAm1, precf1 = np.loadtxt('marsvpl/solarsys.Mars.backward',unpack = True)
fig = plt.figure(figsize=(8.5,4))
fig.subplots_adjust(wspace=0.25)
plt.subplot(1,2,1)
-plt.plot(out.Earth.Time/1e6,out.Earth.Obliquity,'k-',label='Earth with Moon')
-plt.plot(out2.Earth.Time/1e6,out2.Earth.Obliquity,'k--',label='Earth without Moon')
+plt.plot(out.Earth.Time/1e6,out.Earth.Obliquity,'-',c=vplot.colors.pale_blue,label='Earth with Moon')
+plt.plot(out2.Earth.Time/1e6,out2.Earth.Obliquity,'-',c=vplot.colors.dark_blue,label='Earth without Moon')
plt.ylabel(r'$\varepsilon$ ($^{\circ}$)')
plt.xlabel('Time (Myr)')
plt.xlim(0,1)
plt.legend(loc='upper right',fontsize=10)
ax=plt.subplot(1,2,2)
-plt.plot(out.Mars.Time/1e6,np.ones_like((out.Mars.Obliquity))*np.nan,'k-',label='Mars')
+plt.plot(tm1/1e6,oblm1,'-',color=vplot.colors.orange,zorder=100,lw=1,label='Mars (DistOrb+DistRot)')
+plt.plot(tm/1e6,oblm,'-',c=vplot.colors.red,lw=1,label='Mars (HNBody+DistRot)')
+#plt.plot(out4.Mars.Time/1e6,out4.Mars.Obliquity,'k-',label='Mars')
plt.ylabel(r'$\varepsilon$ ($^{\circ}$)')
plt.xlabel('Time (Myr)')
-plt.xlim(0,1)
-plt.ylim(0,10)
-plt.text(0.2,0.5,'Mars evolution here?',transform=ax.transAxes)
+plt.xlim(-10,0)
+plt.ylim(5,55)
+plt.legend(loc='upper right', fontsize=10)
+# plt.ylim(0,10)
+#plt.text(0.2,0.5,'Mars evolution here?',transform=ax.transAxes)
plt.savefig('ObliqComp.pdf')
plt.close()
