Home > atmlab > h2o > parametrisations > eps_ice_Jiang_04.m

eps_ice_Jiang_04

PURPOSE ^

eps_ice_Jiang_04 Dielectric constant for pure ice according to Jiang 2004

SYNOPSIS ^

function e = eps_ice_Jiang_04( f, t )

DESCRIPTION ^

 eps_ice_Jiang_04   Dielectric constant for pure ice according to Jiang 2004

    Provides the complex dielectric constant following the Jiang 2004
    paper.

    Frequency limit is:
      f: [ 1  MHz, 3000 GHz] 
    The actual limits of temperature is not known. The following limit 
    is applied here
      t: [ 20 K, 280 K ] 


 FORMAT    e = eps_ice_Jiang_04( 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_Jiang_04.m

SOURCE CODE ^

0001 % eps_ice_Jiang_04   Dielectric constant for pure ice according to Jiang 2004
0002 %
0003 %    Provides the complex dielectric constant following the Jiang 2004
0004 %    paper.
0005 %
0006 %    Frequency limit is:
0007 %      f: [ 1  MHz, 3000 GHz]
0008 %    The actual limits of temperature is not known. The following limit
0009 %    is applied here
0010 %      t: [ 20 K, 280 K ]
0011 %
0012 %
0013 % FORMAT    e = eps_ice_Jiang_04( f, t )
0014 %
0015 % OUT   e   Complex dielectric constant
0016 % IN    f   Frequency                     [Hz]
0017 %       t   Temperature                   [K]
0018 
0019 % 2014-06-19   Created by Maryam Jamali
0020 
0021 
0022 function e = eps_ice_Jiang_04( f, t )
0023 
0024 if t<20  |  t>280
0025   error('Valid range for temperature is 20-280 K'); 
0026 end
0027 
0028 t = t - 273.15;
0029 
0030 theta = ( 300 ./ (273.15 + t) ) - 1;
0031 fghz  = f/1e9;
0032 
0033 if fghz<0.001  |  fghz>3000
0034   error('Valid range for frequency is 0.001-3000 GHz'); 
0035 end
0036 
0037 
0038 a     = 1e-4 .* (50.4 + 62 .* theta) .* exp(-22.1 .* theta);
0039 b     = 1e-4 .* ((0.502 - 0.131.*theta) ./ (1 + theta)) + 0.542e-6 .* ((1+theta) ./ (theta+0.0073)).^2;
0040 
0041 e     = 3.15 + i * ( a./fghz + b*fghz + 1.16e-11*fghz.^3);

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