function y=qice(y)
x=y(:,1);
p=y(:,2);
es=(1.0003+4.18e-6*p).*6.1115.*exp(22.452.*x./(x+272.55));
y=es*622./(p-.378*es);