function Ph_plotter(P,h)
open Ph_diagram.fig
hold on
semilogy(h,P)
axis([-4.1e5 -3.92e5 1e6 2e7])
text(-4e5,3.8e6,'$0\unit{\dcu}$')
hold off
