aboutsummaryrefslogtreecommitdiffstats
path: root/mkcuts_dvpT.py
diff options
context:
space:
mode:
Diffstat (limited to 'mkcuts_dvpT.py')
-rw-r--r--mkcuts_dvpT.py22
1 files changed, 22 insertions, 0 deletions
diff --git a/mkcuts_dvpT.py b/mkcuts_dvpT.py
index ce6bed6..65317d9 100644
--- a/mkcuts_dvpT.py
+++ b/mkcuts_dvpT.py
@@ -120,6 +120,8 @@ def main():
rhoX = np.log10(rhoX)
preZ = np.log10(preZ)
preX = np.log10(preX)
+ lTX = np.log10(TX)
+ lTZ = np.log10(TZ)
# Mark Arrays vs. NaN
rhoZ = np.ma.masked_array(rhoZ, mask=np.isnan(rhoZ))
@@ -219,6 +221,22 @@ def main():
plt.xlabel('Y [AU]')
plt.ylabel('Z [AU]')
+ f9 = plt.figure(9)
+ plt.imshow(lTZ, cmap=tmp_cm, extent=ext, interpolation='none')
+ plt.colorbar()
+ plt.grid()
+ plt.title('Log10 Temperature Cut [K], t=%.2f' % output.info["time"])
+ plt.xlabel('X [AU]')
+ plt.ylabel('Y [AU]')
+
+ f10 = plt.figure(10)
+ plt.imshow(lTX, cmap=tmp_cm, extent=ext, interpolation='none')
+ plt.colorbar()
+ plt.grid()
+ plt.title('Log10 Temperature Cut [K], t=%.2f' % output.info["time"])
+ plt.xlabel('Y [AU]')
+ plt.ylabel('Z [AU]')
+
# Save Figures
if args.save:
if args.d:
@@ -233,6 +251,8 @@ def main():
if args.T:
f7.savefig('cut_Txy_%05d.png' % iout)
f8.savefig('cut_Tzy_%05d.png' % iout)
+ f9.savefig('cut_logTxy_%05d.png' % iout)
+ f10.savefig('cut_logTzy_%05d.png' % iout)
# Show Figures?
if args.show:
@@ -251,6 +271,8 @@ def main():
if args.T:
plt.close(f7)
plt.close(f8)
+ plt.close(f9)
+ plt.close(f10)
"""
Jump into Main().