# Latest commands 2012-06-04: # cd ~/Seismo/gcf/ scpe '2012060?_????' gcf2ts `ls 3u93z2/2012060*`
SORRY, gcf2ts is completely revised.
New version includes the ts-append stage # # # ##################################################################### # # HOW TO process gcf data until a low-pass filtered and decimated # time series is obtained # # scpe YYYYMMDD_hhmm # or scpe YYYYMMDD_hh'??' # # or scpe YYYYMMDD. # # scp from elder 3-d gcf into the three folders 3u93?2/ # We always retrieve a full hour # # gcf2ts 3u93e2/20110611_2200e.gcf # or for 10-minute bunches from 1-minute gcf: # # foreach f ( 3u93e2/20110611_22?0e.gcf ) # gcf2ts $f # end # # converts all files of that hour into *.sac and then into *.ts # cutting them up into 10-minute pieces. # It will reprocess all *.sac files in the directory. # Remove them if you want to speed up the process. # (A remove-sac option in gcf2ts would be helpful) # # appgcfts -1Hz G2011.163/GCF.2011:163:23:40:00.3U93N2.ts appgcfts G2011.163/GCF.2011:163:23:40:00.3U93N2.ts appgcfts -100Hz G2011.163/GCF.2011:163:23:40:00.3U93N2.ts # # You can do this with only one file at a time; the main purpose # is filtering, warranting segments that can be abutted later on. # The first call creates a 10-minute file with 1 Hz sampling, # the second a 10-minute file with 10 Hz sampling. # Output to bin ts file where the sampling rate is added to the # input name: G2011.163/GCF.2011:163:23:40:00.3U93N2-10Hz.ts # for the second example. # (appgcfts has no help text yet) # # You may prefer to append ts-files unaltered/unfiltered using tslist-app # THIS IS FAST! foreach c ( E N Z ) tslist-app -O:`label SV,$c` d/tmp.mc -I + \ SAC/G2013.093/GCF.2013:093:03:*.3U93${c}2.ts end # # If you are convinced that you must filter / resample, # you have to accept a rather long processing time # appgcfts -1Hz -A OS-2011-06-12-22-40-E-1Hz.tsf \ G2011.163/GCF.2011:163:22:[4-5]0:00.3U93E2.ts \ G2011.163/GCF.2011:163:23:[0-1]0:00.3U93E2.ts # # appends several 10-min sections and produces an ascii file. # The collecting ascii file must be specified after the -A # Option -a continues to collect into a pre-existent file, # while -A requests to start all over again from the head. # # Instead of manually specifying all files, you can process an # unbroken sequence of files using # foreach c ( E N Z ) appgcfts -1Hz -A OS-2011-06-12-22-00-${c}-1Hz.tsf \ `chain-of-gcf G2011.162/GCF.2011:162:22:00:00.3U93${c}2.ts` end # # A tslist command to read such tsf files and creating an mc-file would be # tslist OS-2011-06-12-22-00-N-1Hz.tsf -k4 -gi4,2i3,i4,2i3,i4,t43,f12.0 \ # -I -O:`label SEIS,N` OS-2011-06-12-22-00-1Hz.mc # # This operation is integrated in appgcfts: # # Either after the collecting step appgcfts -O:`label SV,E` d/sv-20110616-000020-120959.mc -I d/app.tsf # ... or foreach c ( E N Z ) appgcfts -1Hz -O`label SEIS,$c` d/OS-2011-06-12-22-00-1Hz.mc \ -A d/OS-2011-06-12-22-00-${c}-1Hz.tsf \ `chain-of-gcf G2011.162/GCF.2011:162:22:00:00.3U93${c}2.ts` end # # For Ludger: foreach c ( E N Z ) appgcfts -100Hz -E del.tse,DN -A OS-2011-06-11-22-20-${c}-100Hz.tsf \ `chain-of-gcf G2011.162/GCF.2011:162:22:20:00.3U93${c}2.ts` end foreach c ( E N Z ) appgcfts -100Hz -E del.tse,DS -A OS-2011-06-15-00-20-${c}-100Hz.tsf \ `chain-of-gcf G2011.166/GCF.2011:166:00:20:00.3U93${c}2.ts` end gcf2ts uu/jun/29/3u93e2/201106291320003u93e2.gcf