function result = epswater93(fGHz, TK)

% Dielectric permittivity of liquid water without salt according
% to Liebe et al. 1991 Int. J. IR+mm Waves 12(12), 659-675
% Mätzler, June 2002

TETA=1-300/TK;
e0=77.66-103.3*TETA;
e1=0.0671*e0;
f1=20.2+146.4*TETA+316*TETA.*TETA;
e2=3.52;  %+7.52*TETA, i.e. version of Liebe 1993
f2=39.8*f1;
eps=e2+(e1-e2)./(1-i*fGHz./f2)+(e0-e1)./(1-i*fGHz./f1);

result=eps;