Difference between revisions of "Rayleigh scattering"

From Planets
Jump to: navigation, search
Line 20: Line 20:
 
TAURAY(NW) is calculated in calc_rayleigh.F90
 
TAURAY(NW) is calculated in calc_rayleigh.F90
  
TAURAY(NW) = <math> \displaystyle \frac{\int \tau(\lambda) B_{\lambda} \, \mathrm{d}\lambda}{} </math>
+
TAURAY(NW) = <math> \displaystyle \frac{\int \tau(\lambda) B_{\lambda} \, \mathrm{d}\lambda}{\frac{\int B_{\lambda} \, \mathrm{d}\lambda} </math>

Revision as of 14:28, 28 September 2022

About Rayleigh scattering in LMDZ Generic

Formalism

References

Hansen (1974) : https://ui.adsabs.harvard.edu/link_gateway/1974SSRv...16..527H/ADS_PDF

Rayleigh routine in exo_k : http://perso.astrophy.u-bordeaux.fr/~jleconte/exo_k-doc/_modules/exo_k/rayleigh.html#Rayleigh.sigma_mol

Exo_k uses formalism from : Caldas (2019) : https://hal.archives-ouvertes.fr/hal-02005332/document

Equation

in optcv.F90 :

TRAY(K,NW) = TAURAY(NW) * DPR(K)

TAURAY(NW) is calculated in calc_rayleigh.F90

TAURAY(NW) = \( \displaystyle \frac{\int \tau(\lambda) B_{\lambda} \, \mathrm{d}\lambda}{\frac{\int B_{\lambda} \, \mathrm{d}\lambda} \)