Re: [MSNoise] "keep_all" in s001configurator.py

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 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)) ----------------------------------------------------------------