Home > atmlab > mie > epsice.m

epsice

PURPOSE ^

Function for calculating the relative permittivity of pure ice in the

SYNOPSIS ^

function z = epsice(f,TK)

DESCRIPTION ^

 Function for calculating the relative permittivity of pure ice in the 
 microwave region, according to C. M�zler, "Microwave properties of ice 
 and snow", in B. Schmitt et al. (eds.) Solar System Ices, Astrophys. 
 and Space Sci. Library, Vol. 227, Kluwer Academic Publishers, 
 Dordrecht, pp. 241-257 (1998). Input:
 f = frequency in GHz, range 0.01 to 3000
 TK = temperature (K), range 20 to 273.15

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

epsice.m

SOURCE CODE ^

0001 function z = epsice(f,TK)
0002 % Function for calculating the relative permittivity of pure ice in the
0003 % microwave region, according to C. M�zler, "Microwave properties of ice
0004 % and snow", in B. Schmitt et al. (eds.) Solar System Ices, Astrophys.
0005 % and Space Sci. Library, Vol. 227, Kluwer Academic Publishers,
0006 % Dordrecht, pp. 241-257 (1998). Input:
0007 % f = frequency in GHz, range 0.01 to 3000
0008 % TK = temperature (K), range 20 to 273.15
0009 
0010     a = dbstack;
0011     %
0012     if length(a)==1 | ~strncmp(a(2).file,'mie',3) | ~strncmp(a(2).file,'eps',3)
0013       error('This function can just be used by the Mie functions.');
0014     end
0015     
0016     B1 = 0.0207;
0017     B2 = 1.16e-11;
0018     b = 335;
0019     deltabeta = exp(-10.02 + 0.0364*(TK-273));
0020     betam = (B1/TK) * ( exp(b/TK) / ((exp(b/TK)-1)^2) ) + B2*f^2;
0021     beta = betam + deltabeta;
0022     theta = 300 / TK - 1;
0023     alfa = (0.00504 + 0.0062*theta)*exp(-22.1*theta);
0024     z = 3.1884 + 9.1e-4*(TK-273);
0025     z = z + i*(alfa/f + beta*f);

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