function es=svp(t,p)
%computes the saturation vapor pressure in mb (Buck, 1981)
%above sea water, esea=es(1-0.000537*Salinity)
es=6.1121.*exp(17.502.*t./(t+240.97)).*(1.0007+3.46e-6*p);  %saturation vapor pressure (mb)
