ref_meltp.f90 File Reference

Functions/Subroutines

subroutine ref_meltp (F, fmelt, phw, phs, crefd)
 Compute the refractive index of a metling particle following the Fabry-Szyrmer (1999) core-shell model 5. The particles use a concentric core and shell with different dielectric properties. The particle core consists of ice inclusions in a liquid water matrix, which together serve as a matrix for air inclusions. The outer shell is modeled as ice inclusions in a water matrix, which together form inclusions in a matrix of air. The effective final dielectric constant is then approx. by the original Maxwell-Garnett theory More...
 

Function/Subroutine Documentation

subroutine ref_meltp ( real, intent(in)  F,
real, intent(in)  fmelt,
real, intent(in)  phw,
real, intent(in)  phs,
complex(8), intent(out)  crefd 
)

Compute the refractive index of a metling particle following the Fabry-Szyrmer (1999) core-shell model 5. The particles use a concentric core and shell with different dielectric properties. The particle core consists of ice inclusions in a liquid water matrix, which together serve as a matrix for air inclusions. The outer shell is modeled as ice inclusions in a water matrix, which together form inclusions in a matrix of air. The effective final dielectric constant is then approx. by the original Maxwell-Garnett theory

Parameters
F[in] Frequency [GHz]
fmelt[in] Melted fraction (0=snow, 1=water) Note that the input value must be >0 and <1, use ref_hydro.f90 otherwise
phw[in] Water density [kg/m3]
phs[in] Snow density [kg/m3]
crefd[out] Complex refractive index of liquid water
See Also
watoptic (water)
iceoptic (ice)
mg_ellips (model matrix/inclusions)
Author
Bernat P Treserras (berna.nosp@m.t.pt.nosp@m.reser.nosp@m.ras@.nosp@m.mcgil.nosp@m.l.ca) McGill University
Date
September 13th, 2016 - Initial version
From bbmodel.for May 1995 - Frederic Fabry