function y=qsea_p(x)
% computes saturation e at sea surface given x = [SST P]
ts=x(:,1); % sea temp
p=x(:,2);  % pressure mb
es=6.112.*exp(17.502.*ts./(ts+241.0))*.98.*(1.0007+3.46e-6.*p);
y=es*622./(p-.378*es);