Hi Thomas,
Thank you very much, Thomas. I should have downloaded from link at
. 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@ujf-grenoble.fr
____________________________________________________