Home > atmlab > h2o > parametrisations > eps_ice_liebe93.m

eps_ice_liebe93

PURPOSE ^

EPS_ICE_LIEBE93 Dielectric constant for pure ice according to Liebe 93

SYNOPSIS ^

function e = eps_ice_liebe93( f, t )

DESCRIPTION ^

 EPS_ICE_LIEBE93   Dielectric constant for pure ice according to Liebe 93

    Provides the complex dielectric constant following the Liebe 1993
    paper.

    The actual limits of the parameterisation are not known. The following
    limits are here applied:
      f: [ 10 MHz, 1000 GHz] 
      t: [ 20 K, 280 K ] 

    Another parameterisation is provided by *epsice*.

 FORMAT    e = eps_ice_liebe93( f, t )
        
 OUT   e   Complex dielectric constant
 IN    f   Frequency [Hz]
       t   Temperature [K]

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

eps_ice_liebe93.m

SOURCE CODE ^

0001 % EPS_ICE_LIEBE93   Dielectric constant for pure ice according to Liebe 93
0002 %
0003 %    Provides the complex dielectric constant following the Liebe 1993
0004 %    paper.
0005 %
0006 %    The actual limits of the parameterisation are not known. The following
0007 %    limits are here applied:
0008 %      f: [ 10 MHz, 1000 GHz]
0009 %      t: [ 20 K, 280 K ]
0010 %
0011 %    Another parameterisation is provided by *epsice*.
0012 %
0013 % FORMAT    e = eps_ice_liebe93( f, t )
0014 %
0015 % OUT   e   Complex dielectric constant
0016 % IN    f   Frequency [Hz]
0017 %       t   Temperature [K]
0018 
0019 % 2004-10-22   Created by Patrick Eriksson
0020 
0021 
0022 function e = eps_ice_liebe93( f, t )
0023 
0024 
0025 theta = 300 ./ t;
0026 fghz  = f/1e9;
0027 
0028 if fghz<0.01  |  fghz>1000
0029   error('Valid range for frequency is 0.01-1000 GHz'); 
0030 end
0031 if t<20  |  t>280
0032   error('Valid range for temperature is 20-280 K'); 
0033 end
0034 
0035 
0036 a     = ( theta - 1 ) * exp( 17 - 22.1*theta );
0037 b     = 1e-5 * ( (0.233/(1-0.993/theta)).^2 + 6.33/theta - 1.31 ); 
0038 
0039 e     = 3.15 + i * ( a./fghz + b*fghz );

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