Home > atmlab > h2o > parametrisations > eps_ice_matzler06.m

eps_ice_matzler06

PURPOSE ^

Function for calculating the complex relative permittivity

SYNOPSIS ^

function z = eps_ice_matzler06(f,TK)

DESCRIPTION ^

 Function for calculating the complex relative permittivity 
 of pure ice in the microwave region, 
 according to 
 Matzler, C., P.W. Rosenkranz, A. Battaglia and J.P. Wigneron (eds.), 
 "Thermal Microwave Radiation - Applications for Remote Sensing", 
 IEE Electromagnetic Waves Series 52, London, UK (2006), Chapter 5.

 FORMAT   z = eps_ice_matzler06(f,TK)
        
 OUT   z   Complex relative permittivity
 IN    f   Frequency [Hz]. 0.01-3000 GHz.
       TK  Temperature [K]. 20-273.15 K.

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

eps_ice_matzler06.m

SOURCE CODE ^

0001 % Function for calculating the complex relative permittivity
0002 % of pure ice in the microwave region,
0003 % according to
0004 % Matzler, C., P.W. Rosenkranz, A. Battaglia and J.P. Wigneron (eds.),
0005 % "Thermal Microwave Radiation - Applications for Remote Sensing",
0006 % IEE Electromagnetic Waves Series 52, London, UK (2006), Chapter 5.
0007 %
0008 % FORMAT   z = eps_ice_matzler06(f,TK)
0009 %
0010 % OUT   z   Complex relative permittivity
0011 % IN    f   Frequency [Hz]. 0.01-3000 GHz.
0012 %       TK  Temperature [K]. 20-273.15 K.
0013 
0014 % 2006-03-07   Function provided by C. Matzler (some adaptions by PE).
0015 
0016 function z = eps_ice_matzler06(f,TK)
0017 
0018     f = f/1e9;
0019 
0020     if f<0.01  |  f>3000
0021       error('Valid range for frequency is 0.01-3000 GHz'); 
0022     end
0023     if TK<20  |  TK>273.15
0024       error('Valid range for temperature is 20-273.15 K'); 
0025     end
0026     
0027     B1 = 0.0207;
0028     B2 = 1.16e-11;
0029     b = 335;
0030     deltabeta = exp(-9.963 + 0.0372.*(TK-273));
0031     betam = (B1./TK).* ( exp(b./TK)./ ((exp(b./TK)-1).^2) ) + B2*f.^2;
0032     beta = betam + deltabeta;
0033     theta = 300./TK - 1;
0034     alfa = (0.00504 + 0.0062*theta).*exp(-22.1*theta);
0035     z = 3.1884 + 9.1e-4*(TK-273);
0036     z = z + i*(alfa./f + beta.*f);

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