in function [es, s] = satvap(T), the derivative of es to temperature s (kPa/C) needs a correction

[bobzsu]

in

a = 7.5; b = 237.3; %degrees C es = 6.10710.^(a.T./(b+T)); s = eslog(10)a*b./(b+T).^2;

Suggest to change this to eq. (6), Murray (1967) On the computation of water vapor pressure, J. Appl. Meteorol.,

es0 = 6.1078;

over water: a = 17.2693882; b = 35.86;

over ice: a = 21.8745584; b = 7.66;

es = es0exp(a.(T-273.16)./(T-b)); [ des/dT = es.((273.16-b)./(T-b).^2 ] s = es.((273.16-b)./(T-b).^2 ;

Somebody please double check