function result = epswater93(fGHz, TK)

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

    a = dbstack;
    %
    if length(a)==1 | ~strncmp(a(2).file,'mie',3) | ~strncmp(a(2).file,'eps',3)
      error('This function can just be used by the Mie functions.');
    end

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;