Hi Chris,

The design of MSNoise is that you do the moving stack of 60 days (=smoothing over the last 59 days) and compare that to a REF. This is achieved with "mov stack"= "msnoise stack -m". If you want to make the step stacking , i.e. only exporting the stack for each nth day, you need to run the "msnoise stack -s" command, but it's really not well documented and implemented, e.g. the dvv plots will probably fail.


The fact that the first 60 days are the same are a choice that goes back to the orginal codes I implemented, assuming the stack should be the same as long as the number of days of movstacks are not available. It's a bad implementation and I really would like to change this in the future.

The "-i 9999" command is only useful if you want to re-stack data computed more than 1 day ago.

Let me know how it goes

Thomas


On 17/09/2016 11:49, 蒋承孜 wrote:
Hello,


I came across a stack error when using MSNoise and could not find the solution 
in the documentation. 
Ihavefinishedthecomputationofthecross-correlationsformydataandnowIwouldliketo 
do60-daystack with10-dayintervals. However, itreturned 
wrongresultswhenItriedtorun 'msnoisestack -s -i999'. Forinstance, 
Iwouldliketostackcross-correlationfunctionscomputedfrom2009-09-10to2011-07-02. 
Theresulton2009-09-10wascorrect, buttheremainingstacks were exactthesame 
with2009-09-10. Theprocesswas pastedbelow. Iamnotsurewhatiswrong. HowcanI 
resolveit?


Thankyouinadvance,


Chris Chiang
Undergraduate Student
School of Geophysics and Information Technology
China University of Geosciences, Beijing
29 Xueyuan Road, Haidian District, Beijing, China








PROCESS:


(E:\Anaconda\Anaconda2) E:\Data\Mountain\filter0.1-0.2>msnoisestack -s -i999
LetsSTACK !
2016-09-1610:14:40 [DEBUG] Startingthestepstack
[]
2016-09-1610:14:40 [DEBUG] ProcessingYP_NE41:YP_NE64-ZZ-1
2016-09-1610:14:40 [DEBUG] Found650updateddays
2016-09-1610:14:40 [DEBUG] NewDataforYP_NE41:YP_NE64-ZZ-1
2016-09-1610:14:44 [DEBUG] YP_NE41_YP_NE642009-09-10 [2009-09-10 - 2009-11-08] 
(60daystack)
e:\anaconda\anaconda2\lib\site-packages\obspy\io\sac\util.py:371: UserWarning: 
Oldheaderhasinvalidreftime.
   warnings.warn(msg)
2016-09-1610:14:44 [DEBUG] YP_NE41_YP_NE642009-11-09 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:44 [DEBUG] YP_NE41_YP_NE642010-01-08 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:45 [DEBUG] YP_NE41_YP_NE642010-03-09 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:45 [DEBUG] YP_NE41_YP_NE642010-05-08 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:45 [DEBUG] YP_NE41_YP_NE642010-07-07 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:45 [DEBUG] YP_NE41_YP_NE642010-09-05 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642010-11-04 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642011-01-03 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642011-03-04 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642011-05-03 [2009-09-10 - 2009-11-08] 
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642011-07-02 [2009-09-10 - 2009-11-08] 
(60daystack
_______________________________________________
MSNoise mailing list
MSNoise@mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise

_______________________________________________
MSNoise mailing list
MSNoise@mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise

Reply via email to