Home > atmlab > mie > epsoil.m

epsoil

PURPOSE ^

Dielectric permittivity soil according to Dobson et al. 1985,

SYNOPSIS ^

function result = epsoil(fGHz, TK, rob, mv, S, C)

DESCRIPTION ^

 Dielectric permittivity soil according to Dobson et al. 1985,
 as formulated by Ulaby et al. 1986, Monograph on
 "Microwave Ramote Sensing", Vol. 3, p. 2102.
 Input: 
 fGHz: frequency, GHz 
 TK: temperature, K 
 rob: bulk density, g/cm^3
 mv: volumetric soil moisture
 S, C: Sand and Clay fractions, respectively
 M�zler, June 2002

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

epsoil.m

SOURCE CODE ^

0001 function result = epsoil(fGHz, TK, rob, mv, S, C)
0002 
0003 % Dielectric permittivity soil according to Dobson et al. 1985,
0004 % as formulated by Ulaby et al. 1986, Monograph on
0005 % "Microwave Ramote Sensing", Vol. 3, p. 2102.
0006 % Input:
0007 % fGHz: frequency, GHz
0008 % TK: temperature, K
0009 % rob: bulk density, g/cm^3
0010 % mv: volumetric soil moisture
0011 % S, C: Sand and Clay fractions, respectively
0012 % M�zler, June 2002
0013 
0014 
0015     a = dbstack;
0016     %
0017     if length(a)==1 | ~strncmp(a(2).file,'mie',3) | ~strncmp(a(2).file,'eps',3)
0018       error('This function can just be used by the Mie functions.');
0019     end
0020 
0021 ess=4.7;              % permittivity of solid material
0022 ross=2.65;            % density of solid material, g/cm^3
0023 alfa=0.65;
0024 beta=1.09-0.11*S+0.18*C;
0025 po=(ross-rob)./ross;  %  Porosity
0026 ew=epswater(fGHz, TK);
0027 vss=1-po;
0028 va=po-mv;
0029 ealfa=(1-po)*ess.^alfa+po-mv+(mv.^beta).*(ew.^alfa);
0030 eps=ealfa.^(1/alfa);
0031 result=eps;

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