function y=qair3(y)
p=y(:,1);%pressure in mb
x=y(:,2);%temperature in C
h=y(:,3);%relative humidity in %
es=6.112.*exp(17.502.*x./(x+241.0)).*(1.0007+3.46e-6*p).*h/100;
y=es*622./(p-.378*es);