function qsat = qs(p,Temp)

%  qs(p,Temp) is saturation mixing ratio based on Wexler's formula for es.
%   From Bolton, 1980, MWR, 108, 1046-1053.

  Rd    = 287.04;
  Rv    = 461.50;
  esat = 611.2*exp(17.67*(Temp-273)./(Temp-29.5));
  qsat = (Rd/Rv)*esat./(p-esat);