[MSNoise] "keep_all" in s001configurator.py

[Kohtaro R. Araragi]

Hi Thomas,

Thank you very much, Thomas. I should have downloaded from link at
https://github.com/ROBelgium/MSNoise/releases NOT from
https://github.com/ROBelgium/MSNoise. I will check the release version
and let you know if you find anything around "keep_all" options or other
things.

Sincerely,
Run

(11/11/14 1:20 PM), Thomas Lecocq wrote:
> Hello Run,
>
> Please send your emails to the mailing list, your remarks could be
> useful for many users !
>
>
>
> Le 22/10/2014 11:47, Kohtaro R. Araragi a e'crit :
>> Hi Thomas,
>>
>> Thank you very much for your help the other day. I had to change a
>> little bit after that, but your e-mail with other MSNoise users in
>> the past can help me to resolve all prolbem and now I can use all
>> function until s06compute_dtt.py.
>>
>> Today, I would like to you let you know about what I find in MSNoise
>> function.
>>
>> In "s001configurator.py", we can set "keep_all" as Y if we want to
>> keep each corr_duration CCF (or all 30 seconds cross-corr). However,
>> when I set Y in "keep_all", "s03compute_cc.py" did not work. When I
>> check the code and error message, it seems that current version of
>> s03compute_cc.py use "thisdate" for processing(see bold part in the
>> source below) of "keep_all" before value is assigned. When I set
>> "keep_all" as "N", everything is working and I can compute C.C. I may
>> make some mistakes, and I would be happy if you let me know where I
>> can change.
> You are using the "master" version of MSNoise, which contains this bug
> ! use "releases" on github to be sure to use the bug-free versions !
>
> Thomas
>
>>
>> You may have updated this part and if so, I appreciate if you let me
>> know.
>>
>> Sincerely,
>> Run
>>
>> ----------------------------------------------------------------
>> *if keep_all:**
>> **add_corr(db, station1.replace('.', '_'), station2.replace(**
>> **'.', '_'), filter.ref, thisdate, thistime, min30 / fe, "ZZ", corr,
>> fe)**
>> * if keep_days:
>> if not np.any(np.isnan(corr)) and not np.any(np.isinf(corr)):
>> daycorr += corr
>> ndaycorr += 1
>>
>> if keep_days:
>> thisdate = time.strftime(
>> "%Y-%m-%d", time.gmtime(basetime))
>> thistime = time.strftime(
>> "%H_%M", time.gmtime(basetime))
>> add_corr(
>> db, station1.replace(
>> '.', '_'), station2.replace('.', '_'), filter.ref,
>> thisdate, thistime, min30 / fe, 'ZZ', daycorr, fe, day=True,
>> ncorr=ndaycorr)
>> update_job(db, goal_day, orig_pair, 'CC', 'D')
>> logging.info("Job Finished. It took %.2f seconds" % (time.time() - jt))
>> ----------------------------------------------------------------
>


-- 

____________________________________________________
K. Run Araragi
PhD studnet/M.Sc.
ISTerre, France & 
Earthquake Research Institute
1-1-1 Yayoi, Bunkyo-ku
Tokyo, Japan 113-0032

Skype name: kalessinlord
Tel: 81-3-5841-5697
Cel: 81-90-1166-3091
     33-970-461-805 
mail:kohtaro.araragi at ujf-grenoble.fr
____________________________________________________