Difference between revisions of "Rayleigh scattering"

Revision as of 14:58, 28 September 2022

in optcv.F90 :

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

Let's write it :

TRAY \( \displaystyle = \tau_{RAY}(\lambda) dP\)

In exo_k we have :

TRAY \( \displaystyle = \sigma_{exok} dN \) with dm in kg/m2

TRAY \( \displaystyle = \sigma_{exok} \frac{dm}{m_{molecule}} \) with dm in kg/m2

TRAY \( \displaystyle = \frac{\sigma_{exok}}{g * m_{molecule}} dP\)

so

\( \displaystyle \tau_{RAY}(\lambda) = \frac{\sigma_{exok}}{g * m_{molecule}} \)

To be noticed :

We write \( \displaystyle TRAY = \tau_{RAY} dP\)

TAURAY(NW) is calculated in calc_rayleigh.F90

TAURAY(NW) \[ \displaystyle = \frac{\int_{\lambda'} TAUVAR(\lambda') B_{\lambda} \, \mathrm{d}\lambda'}{\int B_{\lambda} \, \mathrm{d}\lambda'} \]

\( \displaystyle \tau(\lambda) \) is called TAUVAR

@@ Line 41: / Line 41: @@
 TAURAY(NW) is calculated in calc_rayleigh.F90
-TAURAY(NW) : <math> \displaystyle \tau_{RAY} = \frac{\int_{\lambda'} \tau(\lambda) B_{\lambda} \, \mathrm{d}\lambda'}{\int B_{\lambda} \, \mathrm{d}\lambda'} </math>
+TAURAY(NW) : <math> \displaystyle = \frac{\int_{\lambda'} TAUVAR(\lambda') B_{\lambda} \, \mathrm{d}\lambda'}{\int B_{\lambda} \, \mathrm{d}\lambda'} </math>
 <math> \displaystyle \tau(\lambda) </math> is called TAUVAR