In [ ]:

    

%matplotlib inline
import yt
import logging
logging.getLogger('yt').setLevel(logging.ERROR)
import numpy as np
import matplotlib
matplotlib.rcParams['figure.dpi'] = 150
import matplotlib.pyplot as plt
from tools import setup_cl


    

In [ ]:

    

dirs = [
        '/home/ychen/data/0only_1110_h0_rerun/',
        '/home/ychen/data/0only_0602_hydro/',
        '/home/ychen/data/0only_0529_h1/',
        '/home/ychen/data/0only_0605_hinf/',
        '/home/ychen/data/0only_1212_h0_10Myr_rerun/',
        '/home/ychen/data/0only_0413_hydro_10Myr/',
        '/home/ychen/data/0only_1022_h1_10Myr/',
        '/home/ychen/data/0only_0204_hinf_10Myr/']

colors, labels = setup_cl(dirs)

for dir in dirs:
    data = np.loadtxt(dir+'GridAnalysis_Xray_100kpc.txt')
    label = labels[dir]
    if '10Myr' in dir:
        plt.plot(data[:,1], data[:,2], color=colors[dir])
        #plt.fill_between(data[:,1], data[:,2], data[:,3],\
        #                 facecolor=colors[dir], lw=1, linestyle=':', edgecolor=colors[dir], alpha=0.3)
        #plt.plot(data[:,1], data[:,2], c=colors[dir])
    else:
        plt.plot(data[:,1], data[:,2], color=colors[dir], label=label)
    #plt.fill_between(data[:,1], data[:,4], -data[:,5],\
    #                 facecolor=colors[dir], lw=1, linestyle=':', edgecolor=colors[dir], alpha=0.3)

handles, labels = plt.axes().get_legend_handles_labels()
# reverse the order
plt.legend(handles[::-1], labels[::-1], loc=1)

#plt.ylabel('z (kpc)')
plt.xlabel('t (Myr)')
plt.xlim(0,500)
plt.axvline(10, ls=':', lw=1, alpha=0.5)
plt.ylim(2.30E44,2.55E44)
#plt.savefig('LobeSize_20Myr.pdf')
#plt.semilogx()


    

In [ ]:

    

colors, labels = setup_cl(dirs)
for dir in dirs:
    data = np.loadtxt(dir+'GridAnalysis_Xray_100kpc.txt')
    label = labels[dir]
    if '10Myr' in dir:
        plt.plot(data[:,1], data[:,3]/data[:,2], color=colors[dir])
    else:
        plt.plot(data[:,1], data[:,3]/data[:,2], color=colors[dir], label=label)
    #plt.fill_between(data[:,1], data[:,4], -data[:,5],\
    #                 facecolor=colors[dir], lw=1, linestyle=':', edgecolor=colors[dir], alpha=0.3)

handles, labels = plt.axes().get_legend_handles_labels()
# reverse the order
plt.legend(handles[::-1], labels[::-1], loc=1)

#plt.ylabel('z (kpc)')
plt.xlabel('t (Myr)')
plt.xlim(0,100)
plt.ylim(0.99995, 1)


    

In [ ]: