pl.legend(loc = 'upper right')
pl.axis([1e-9,1e-6,1e-24,1e-19])
pl.savefig('plots/skew_ret_water/skew_ret_water_A0_A2.pdf')
show()
ls4 = 1e9*ls[ 2]#2]
ls5 = 1e9*ls[12]#4]
ls6 = 1e9*ls[22]#6]
ls1 = 1e9*ls[32]#8]
ls2 = 1e9*ls[42]#12]
ls3 = 1e9*ls[52]#16]
fig = pl.figure()
ax = fig.add_axes([0.1,0.1,0.8,0.8])
#pl.semilogy(thetas, G_l_t_dt)
ax.semilogy(thetas, G_l_t_dt[ 2,:], label = r'$\ell$ = %1.2f nm' %ls4)
ax.semilogy(thetas, G_l_t_dt[12,:], label = r'$\ell$ = %1.2f nm' %ls5)
ax.semilogy(thetas, G_l_t_dt[22,:], label = r'$\ell$ = %1.2f nm' %ls6)
ax.semilogy(thetas, G_l_t_dt[32,:], label = r'$\ell$ = %1.2f nm' %ls1)
ax.semilogy(thetas, G_l_t_dt[42,:], label = r'$\ell$ = %1.2f nm' %ls2)
ax.semilogy(thetas, G_l_t_dt[52,:], label = r'$\ell$ = %1.2f nm' %ls3)
#ax.semilogy(0,0,'', label = r'$G_\theta = cos(2\theta)/2sin(\theta)$')
pl.xlabel(r'$Angle\,\,\mathrm{[radians]}$', size = 20)
pl.ylabel(r'$-G(\ell,\theta)\,\,\mathrm{[k_{B}T]}$', size = 20)
pl.axis([0,1.7,1e-10,1.0])
#pl.axis([0,1.7,1e-3,1e4])
#pl.title(r'$\mathrm{-G(\ell,\theta)\,vs.\,angle:\,skewed,\,retarded,\,water}$', size = 20)
pl.legend(loc = 'lower left')
#pl.savefig('plots/skew_ret_water/skew_ret_water_G_vs_theta.pdf')
#show()
pl.savefig('plots/skew_ret_water/G_vs_theta_fixed_l.pdf')
#ax = fig.add_axes([0.1,0.1,0.8,0.8])
##pl.semilogy(thetas, G_l_t_dt)
#pl.semilogy(thetas, G_l_t_dt[ 8,:], label = r'$\ell$ = %1.2f nm' %ls1)
#pl.semilogy(thetas, G_l_t_dt[12,:], label = r'$\ell$ = %1.2f nm' %ls2)
#pl.semilogy(thetas, G_l_t_dt[16,:], label = r'$\ell$ = %1.2f nm' %ls3)
#pl.xlabel(r'$Angle\,\,\mathrm{[radians]}$', size = 20)
#pl.ylabel(r'$G(\ell,\theta)\,\,\mathrm{[k_{B}T]}$', size = 20)
##pl.axis([(1./25)*np.pi,(3./4)*np.pi,105,135])
##pl.title(r'$\mathrm{-G(\ell,\theta)\,vs.\,angle:\,skewed,\,retarded,\,water}$', size = 20)
#pl.legend(loc = 'best')
#pl.savefig('plots/skew_ret_water/skew_ret_water_G_vs_theta.pdf')
#show()
fig = pl.figure()
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8])
#pl.semilogy(thetas, G_l_t_dt)
ax.semilogy(thetas, G_l_t_dt[8, :], label=r'$\ell$ = %1.2f nm' % ls1)
ax.semilogy(thetas, G_l_t_dt[12, :], label=r'$\ell$ = %1.2f nm' % ls2)
ax.semilogy(thetas, G_l_t_dt[16, :], label=r'$\ell$ = %1.2f nm' % ls3)
ax.semilogy(0, 0, '', label=r'$G_\theta = cos(2\theta)/2sin(\theta)$')
pl.xlabel(r'$Angle\,\,\mathrm{[radians]}$', size=20)
pl.ylabel(r'$-G(\ell,\theta)\,\,\mathrm{[k_{B}T]}$', size=20)
pl.axis([0, 1.7, 1e-3, 1e4])
#pl.title(r'$\mathrm{-G(\ell,\theta)\,vs.\,angle:\,skewed,\,retarded,\,water}$', size = 20)
pl.legend(loc='lower left')
#pl.savefig('plots/skew_ret_water/skew_ret_water_G_vs_theta.pdf')
#show()
#fig = pl.figure()
#ax = fig.add_axes([0.1,0.1,0.8,0.8])
#ax.plot(x_peak , y_peak, color = 'g',label = r'Synthesized')# $\epsilon$"($\omega$)', linestyle='-')
#ax.plot(x_nopeak, y_nopeak, color = 'b',label = r'Ab initio')# $\epsilon$"($\omega$)', linestyle='-')
pl.legend(loc='upper right')
pl.axis([1e-9, 1e-6, 1e-24, 1e-19])
pl.savefig('plots/compare/srw_thetas_A0_A2.pdf')
show()
ls4 = 1e9 * ls[2] #2]
ls5 = 1e9 * ls[12] #4]
ls6 = 1e9 * ls[22] #6]
ls1 = 1e9 * ls[32] #8]
ls2 = 1e9 * ls[42] #12]
ls3 = 1e9 * ls[52] #16]
fig = pl.figure()
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8])
#pl.semilogy(thetas, G_l_t_dt)
ax.semilogy(thetas, G_l_t_dt[2, :], label=r'$\ell$ = %1.2f nm' % ls4)
ax.semilogy(thetas, G_l_t_dt[12, :], label=r'$\ell$ = %1.2f nm' % ls5)
ax.semilogy(thetas, G_l_t_dt[22, :], label=r'$\ell$ = %1.2f nm' % ls6)
ax.semilogy(thetas, G_l_t_dt[32, :], label=r'$\ell$ = %1.2f nm' % ls1)
ax.semilogy(thetas, G_l_t_dt[42, :], label=r'$\ell$ = %1.2f nm' % ls2)
ax.semilogy(thetas, G_l_t_dt[52, :], label=r'$\ell$ = %1.2f nm' % ls3)
#ax.semilogy(0,0,'', label = r'$G_\theta = cos(2\theta)/2sin(\theta)$')
pl.xlabel(r'$Angle\,\,\mathrm{[radians]}$', size=20)
pl.ylabel(r'$-G(\ell,\theta)\,\,\mathrm{[k_{B}T]}$', size=20)
pl.axis([0, 1.7, 1e-10, 1.0])
#pl.title(r'$\mathrm{-G(\ell,\theta)\,vs.\,angle:\,skewed,\,retarded,\,water}$', size = 20)
pl.legend(loc='lower left')
#pl.savefig('plots/skew_ret_water/G_vs_theta_fixed_l.pdf')
show()
pl.figure()
#ax = fig.add_axes([0.1,0.1,0.8,0.8])
##pl.semilogy(thetas, G_l_t_dt)
#pl.semilogy(thetas, G_l_t_dt[ 8,:], label = r'$\ell$ = %1.2f nm' %ls1)
#pl.semilogy(thetas, G_l_t_dt[12,:], label = r'$\ell$ = %1.2f nm' %ls2)
#pl.semilogy(thetas, G_l_t_dt[16,:], label = r'$\ell$ = %1.2f nm' %ls3)
#pl.xlabel(r'$Angle\,\,\mathrm{[radians]}$', size = 20)
#pl.ylabel(r'$G(\ell,\theta)\,\,\mathrm{[k_{B}T]}$', size = 20)
##pl.axis([(1./25)*np.pi,(3./4)*np.pi,105,135])
##pl.title(r'$\mathrm{-G(\ell,\theta)\,vs.\,angle:\,skewed,\,retarded,\,water}$', size = 20)
#pl.legend(loc = 'best')
#pl.savefig('plots/skew_ret_water/skew_ret_water_G_vs_theta.pdf')
#show()
fig = pl.figure()
ax = fig.add_axes([0.1,0.1,0.8,0.8])
#pl.semilogy(thetas, G_l_t_dt)
ax.semilogy(thetas, G_l_t_dt[ 8,:], label = r'$\ell$ = %1.2f nm' %ls1)
ax.semilogy(thetas, G_l_t_dt[12,:], label = r'$\ell$ = %1.2f nm' %ls2)
ax.semilogy(thetas, G_l_t_dt[16,:], label = r'$\ell$ = %1.2f nm' %ls3)
ax.semilogy(0,0,'', label = r'$G_\theta = cos(2\theta)/2sin(\theta)$')
pl.xlabel(r'$Angle\,\,\mathrm{[radians]}$', size = 20)
pl.ylabel(r'$-G(\ell,\theta)\,\,\mathrm{[k_{B}T]}$', size = 20)
pl.axis([0,1.7,1e-3,1e4])
#pl.title(r'$\mathrm{-G(\ell,\theta)\,vs.\,angle:\,skewed,\,retarded,\,water}$', size = 20)
pl.legend(loc = 'lower left')
#pl.savefig('plots/skew_ret_water/skew_ret_water_G_vs_theta.pdf')
#show()
#fig = pl.figure()
#ax = fig.add_axes([0.1,0.1,0.8,0.8])
#ax.plot(x_peak , y_peak, color = 'g',label = r'Synthesized')# $\epsilon$"($\omega$)', linestyle='-')
#ax.plot(x_nopeak, y_nopeak, color = 'b',label = r'Ab initio')# $\epsilon$"($\omega$)', linestyle='-')