# Latest commands 2012-06-04:
#
cd ~/Seismo/gcf/ 
scpe -d onau '2012060?_????'   
# HARK! SNSN changed the stream name but 3u93 is still our default in Holt catalogs
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 -E guralp.tse,CMGT3AR -10Hz -c z -l 7200 $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