def Phi(mm):
return -1.44296504*log(mm)
phi = array([0, .5, 1., 1.5, 2., 2.5, 3., 3.5, 4]);
phima = Phi(ma)
phisa = Phi(sa)
phimg = Phi(mg)
figure()
plot(phi,phi,'-k')
plot(phi-0.5,phi,'--k') # one sieve larger
plot(phi[0]*ones_like(d00.geom_mean),Phi(d00.geom_mean),'o',color='0.8')
plot(phi[1]*ones_like(d05.geom_mean),Phi(d05.geom_mean),'o',color='0.8')
plot(phi[2]*ones_like(d10.geom_mean),Phi(d10.geom_mean),'o',color='0.8')
plot(phi[3]*ones_like(d15.geom_mean),Phi(d15.geom_mean),'o',color='0.8')
plot(phi[4]*ones_like(d20.geom_mean),Phi(d20.geom_mean),'o',color='0.8')
plot(phi[5]*ones_like(d25.geom_mean),Phi(d25.geom_mean),'o',color='0.8')
plot(phi[6]*ones_like(d30.geom_mean),Phi(d30.geom_mean),'o',color='0.8')
plot(phi[7]*ones_like(d35.geom_mean),Phi(d35.geom_mean),'o',color='0.8')
plot(phi[8]*ones_like(d40.geom_mean),Phi(d40.geom_mean),'o',color='0.8')
plot(phi,phimg,'or')
xlim([-.5,4.5])
xlabel('Sieve size (phi)')
ylabel('DGS geometric mean (phi)')
title('MOF Poking Eyeball Tests, May 3, 2013')
#plt.savefig('D:\crs\proj\poking_eyeball_tests\MOF_test_050313\compare.png')