[MSNoise] Weights used in the MWCS
Thomas Lecocq
thomas.lecocq at seismology.be
Wed Nov 25 07:27:51 UTC 2015
You're welcome
Thomas
Le 23/11/2015 15:19, Thomas Lecocq a écrit :
> 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
>> __________________________________________
>
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
More information about the MSNoise
mailing list