Merge branch 'master' of ssh://gitserver//home/tforb/git.git/Seitosh

src/ts/croposp/tests/Makefile

@@ -20,6 +20,10 @@
 # along with this program. If not, see <http://www.gnu.org/licenses/>.
 # ----
 #
+# run with pairwise noise signal:
+#
+# make all PAIRWISE=yes
+#
 # REVISIONS and CHANGES 
 #    06/02/2019   V1.0   Thomas Forbriger
 # 
@@ -33,6 +37,9 @@ all: compilation.pdf
 CPSD_SET=1 2 3
 CPSD_SET=1 2 3 4
 
+# input data
+CPSD_RAWSET:=$(CPSD_SET)
+
 ALLPDF=PSD.pdf PSDraw.pdf \
   $(addprefix coherence_,$(addsuffix .pdf,$(CPSD_SET))) \
   $(addprefix transfer_,$(addsuffix .pdf,A B)) \
@@ -172,6 +179,15 @@ CPSD_AMP3=100.
 CPSD_AMP4=1.e-2
 CPSD_AMPBG=10.
 
+# test effect of pairwise coherent noise
+ADDITIONAL=
+FILTERBG2=dif  0
+ifdef PAIRWISE
+ADDITIONAL=cpsd2modX.bin
+FILTERBG2=rem
+CPSD_SET=1 3 4
+endif
+
 # modify white noise to channel specific frequency characteristic
 cpsd1mod.bin: cpsd1raw.bin
 	printf "fac  $(CPSD_AMP1)\nhpb  1.,2\nend\n" \
@@ -189,7 +205,7 @@ cpsdbg1mod.bin: cpsdbgraw.bin
 	printf "fac  $(CPSD_AMPBG)\nlp1  2.\navg  0\nhpb  100.,4\nend\n" \
 	  | tidofi -type bin -Type bin -v -o -cs $@ $<
 cpsdbg2mod.bin: cpsdbg1mod.bin
-	printf "dif  0.\nend\n" \
+	printf "$(FILTERBG2)\nend\n" \
 	  | tidofi -type bin -Type bin -v -o -cs $@ $<
 
 %.ascii: %.bin
@@ -202,13 +218,13 @@ DISCARDS=839.
 	  --secondsf $(DISCARDS) $@ $<
 
 # add common noise signal to channel specific noise floor
-cpsd1add.bin: cpsd1modX.bin cpsdbg1modX.bin
+cpsd1add.bin: cpsd1modX.bin cpsdbg1modX.bin $(ADDITIONAL)
 	teseco -v -o -a -type bin -Type bin $@ $^
-cpsd2add.bin: cpsd2modX.bin cpsdbg1modX.bin
+cpsd2add.bin: cpsd2modX.bin cpsdbg1modX.bin 
 	teseco -v -o -a -type bin -Type bin $@ $^
 cpsd3add.bin: cpsd3modX.bin cpsdbg1modX.bin
 	teseco -v -o -a -type bin -Type bin $@ $^
-cpsd4add.bin: cpsd4modX.bin cpsdbg2modX.bin
+cpsd4add.bin: cpsd4modX.bin cpsdbg2modX.bin $(ADDITIONAL)
 	teseco -v -o -a -type bin -Type bin $@ $^
 	
 # ----------------------------------------------------------------------
@@ -284,7 +300,7 @@ CPSD_CROPOSP_STDOPT=--verbose --overwrite --itype bin --trim \
 CPSD_BINFILES=$(addprefix cpsd,$(addsuffix .bin,$(CPSD_SET)))
 # compile list of raw files (apparent instrumental noise and seismic
 # background noise separately)
-CPSD_RAWFILES=$(addprefix cpsd,$(addsuffix mod.bin,$(CPSD_SET) bg1 bg2)) \
+CPSD_RAWFILES=$(addprefix cpsd,$(addsuffix mod.bin,$(CPSD_RAWSET) bg1 bg2)) \
 
 # run complete analysis on the apparent seismometer output
 $(addsuffix _out.xxx,$(OUTPUTS)): $(CPSD_BINFILES)