% Illustration of Bergeron process.
% Brian Hornbuckle, January 22, 2007.

clear;

T = linspace(-10,30,1000);

a = 0.611;  b = 17.502;  c = 240.97;

es = a.*exp(b.*T./(T + c));

b = 21.87;  c = 265.5;

es_ice = a.*exp(b.*T./(T + c));

plot(T,es,'b-',T,es_ice,'g-');
axis([-10 5 0 1]);  grid on;
legend('over water','over ice');