current_t.py

#!/usr/bin/env python

import matplotlib.pyplot as plt
import numpy as np

fig = plt.figure()
#fig.suptitle('Asymptotic current as a function of the temperature')

plt.xlabel('$\kappa_b T/\hbar\Delta$')
plt.ylabel('$I/I_0$')

plt.xscale('log')
plt.grid('on')

#plt.text(1, 0.5, '$\Omega=2\Delta$, $\\alpha=0.005$, $\\rho_0=\\vert z;-\\rangle$', bbox={'facecolor':'white','alpha':0.5})

R_CURRT = (np.loadtxt('RED-STAT-CURR-T.dat')).T
C_CURRT = (np.loadtxt('CP-STAT-CURR-T.dat')).T


rfig = plt.plot(R_CURRT[0],R_CURRT[1],color='red', label='Redfield dynamics')
cfig = plt.plot(C_CURRT[0],C_CURRT[1],color='blue', label='CP dynamics')

plt.legend(('Redfield dynamics', 'CP dynamics'), loc='upper right',
		bbox_to_anchor=(0.95, 0.95))

plt.show()