Home > atmlab > mie > epsalwater.m

epsalwater

PURPOSE ^

Dielectric permittivity of saline water according

SYNOPSIS ^

function result = epsalwater(fGHz, TK, S)

DESCRIPTION ^

 Dielectric permittivity of saline water according
 Meissner and Wentz (2004)
 Frequency range: 1 to 1000 GHz, 
 Input: fGHz: frequency in GHz, TK: temperature in K, 
 Salinity S in promille
 Program by Matzler, Aug. 2004

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

epsalwater.m

SOURCE CODE ^

0001 function result = epsalwater(fGHz, TK, S)
0002 
0003 % Dielectric permittivity of saline water according
0004 % Meissner and Wentz (2004)
0005 % Frequency range: 1 to 1000 GHz,
0006 % Input: fGHz: frequency in GHz, TK: temperature in K,
0007 % Salinity S in promille
0008 % Program by Matzler, Aug. 2004
0009 
0010 
0011     a = dbstack;
0012     %
0013     if length(a)==1 | ~strncmp(a(2).file,'mie',3) | ~strncmp(a(2).file,'eps',3)
0014       error('This function can just be used by the Mie functions.');
0015     end
0016 
0017 
0018 T=TK-273.16; % temp conversion to C
0019 
0020 a0=5.723;
0021 a1=0.022379;
0022 a2=-7.1237e-04;
0023 a3=5.0478;
0024 a4=-7.0315e-02;
0025 a5=6.0059e-04;
0026 a6=3.6143;
0027 a7=2.8841e-02;
0028 a8=0.13652;
0029 a9=1.4825e-03;
0030 a10=2.4166e-04;
0031 
0032 b0=-3.56417e-03;
0033 b1= 4.74868e-06;
0034 b2= 1.15574e-05;
0035 b3= 2.39357e-03;
0036 b4=-3.13530e-05;
0037 b5= 2.52477e-07;
0038 b6=-6.28908e-03;
0039 b7= 1.76032e-04;
0040 b8=-9.22144e-05;
0041 b9=-1.99723e-02;
0042 b10=1.81176e-04;
0043 b11=-2.04265e-03;
0044 b12=1.57883e-04;
0045 
0046 T2=T.*T;
0047 T3=T2.*T;
0048 T4=T3.*T;
0049 S2=S.*S;
0050 alfa0=(6.9431+3.2841*S-0.099486*S2)./(84.850+69.024*S+S2);
0051 alfa1= 49.843-0.2276*S+0.00198*S2;
0052 RTQ=1+alfa0.*(T-15)./(alfa1+T);
0053 R15=S.*(37.5109+5.45216*S+1.4409e-02*S2)./(1004.75+182.283*S+S2);
0054 sigma35=2.903602+8.607e-02*T+4.738817e-04*T2-2.991e-06*T3+4.3047e-09*T4;
0055 sigma=sigma35.*RTQ.*R15;
0056 es0=(3.70886e04-82.168*T)./(421.854+T);
0057 es=es0.*exp(b0.*S+b1.*S2+b2.*T.*S);
0058 v10=(45+T)./(a3+a4.*T+a5.*T2);
0059 v1=v10.*(1+S.*(b3+b4.*T+b5.*T2))
0060 e10=a0+a1*T+a2*T2;
0061 e1=e10.*exp(b6.*S+b7.*S2+b8.*T.*S);
0062 v20=(45+T)./(a8+a9.*T+a10.*T2);
0063 v2=v20.*(1+S.*(b9+b10.*T))
0064 einf0=a6+a7*T;
0065 einf=einf0.*(1+S.*(b11+b12.*T));
0066 e0=8.85418782e-03;
0067 eps=(es-e1)./(1-i*fGHz./v1)+(e1-einf)./(1-i*fGHz./v2)+einf+i*sigma./(2*pi*e0.*fGHz);
0068 
0069 result=eps;

Generated on Mon 15-Sep-2014 13:31:28 by m2html © 2005