From 66ae7f86d6b4be11a152407ff8547af3b778f2e9 Mon Sep 17 00:00:00 2001
From: Volker Hoffmann <volker@cheleb.net>
Date: Mon, 23 Sep 2013 10:08:05 +0200
Subject: make and save log10(T) if T cuts are requested

---
 mkcuts_dvpT.py | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)

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().
-- 
cgit v1.1