croposp [WP]: add reference psd for tests

src/ts/croposp/tests/Makefile

@@ -156,14 +156,16 @@ CPSD_OPT=--DEBUG
 CPSD_OPT=
 OUTPUTS=psd npsd coherency transfer phase
 
-CPSD_SET=1 2 3 4
 CPSD_SET=1 2 3
+CPSD_SET=1 2 3 4
 CPSD_BINFILES=$(addprefix cpsd,$(addsuffix .bin,$(CPSD_SET)))
+CPSD_RAWFILES=$(addprefix cpsd,$(addsuffix mod.bin,$(CPSD_SET) bg)) \
 
 CPSD_CROPOSP_OPT=--nsegments=$(CPSD_NSEGMENTS) \
 		     --divisor=$(CPSD_DIVISOR) \
 		     --padfactor=$(CPSD_PADFACTOR) \
 		     --overlap=$(CPSD_OVERLAP) $(CPSD_LOG)
+
 $(addsuffix _out.xxx,$(OUTPUTS)): $(CPSD_BINFILES)
 	croposp $(CPSD_OPT) --verbose --overwrite --itype bin --trim \
 	  --pattern="%I" \
@@ -176,6 +178,13 @@ $(addsuffix _out.xxx,$(OUTPUTS)): $(CPSD_BINFILES)
 	  --phase phase_out.xxx \
 	  $^
 
+psdraw.xxx: $(CPSD_RAWFILES)
+	croposp $(CPSD_OPT) --verbose --overwrite --itype bin --trim \
+	  --pattern="%F" \
+	  $(CPSD_CROPOSP_OPT) \
+	  --datetolerance 0.1 \
+	  --psd $@ $^
+
 psd%.001.asc: cpsd%.bin
 	foutra -v -o -type bin -ASCII=$(patsubst cpsd%.bin,psd%,$<) \
 	  -power -rbw=0.1 -demean -detrend -logascii=0.1 -avgascii \
@@ -186,6 +195,11 @@ psd%foutra.xxx: psd%.001.asc
 	cat $< >> $@
 
 CPSD_REFFILES=$(addprefix psd,$(addsuffix foutra.xxx,$(CPSD_SET)))
+CPSD_RAWREFFILES=$(addprefix psd,$(addsuffix modfoutra.xxx,$(CPSD_SET) bg))
+
+PSDraw.pdf: psdraw.xxx $(CPSD_RAWREFFILES)
+	$(CROPOSPLOT) --title="$@ \n $(CPSD_CROPOSP_OPT)" \
+	  --titfontsize='small' -v --grid -o $@ $^
 
 PSD.pdf: psd_out.xxx $(CPSD_REFFILES)
 	$(CROPOSPLOT) --title="$@ \n $(CPSD_CROPOSP_OPT)" \