diff options
| -rw-r--r-- | mkcuts_dvpT.py | 22 |
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(). |