[MSNoise] Weights used in the MWCS

Thomas Lecocq thomas.lecocq at oma.be
Mon Nov 23 14:19:18 UTC 2015


Hi,


Le 19/11/2015 07:41, Yoones Vaezi a écrit :
> Hi,
>
> Another problem I am facing is about the weights used for the MWCS 
> program. The weights are supposed to be calculated using equation A5 
> in Clarke et al. 2011 which is
> Inline image 2.
> But the following two lines in your Python code do not follow this 
> equation:
>
> w = 1.0 / (1.0 / (coh[indRange]**2) - 1.0)
> w[coh[indRange] >= 0.99] = 1.0 / (1.0 / 0.9801 - 1.0)
>
>
> Could you please let me know where the values 0.99 and 0.9801 are 
> coming from and what is the reason behind having these two lines in 
> the code?
from the original code of Clarke, so, no, no idea precisely. I assume it 
was to avoid infinite and zero division in the Fortran code.


Thomas

>
> Thank you very much for your time and consideration,
>
> Regards,
>
> Yoones
> __________________________________________



More information about the MSNoise mailing list