croposp [WP]: confirm calibration of PSD by comparison with foutra

src/ts/croposp/Makefile

@@ -74,6 +74,7 @@ clean: ;
 	-find . -name \*.bak | xargs --no-run-if-empty /bin/rm -v
 	-find . -name \*.d | xargs --no-run-if-empty /bin/rm -v
 	-/bin/rm -vf flist *.bin *.o *.xxx *.ps *.pdf *~ $(PROGRAMS)
+	-/bin/rm -vf *.asc
 
 TESTCASEMAKE=$(filter-out %.bak,$(wildcard testcases/Makefile*))
 
@@ -243,9 +244,11 @@ CPSD_DIVISOR=100
 CPSD_PADFACTOR=1
 CPSD_NLOG=20
 CPSD_LOG=--log $(CPSD_NLOG)
+CPSD_OPT=--DEBUG
+CPSD_OPT=
 OUTPUTS=psd npsd coherency transfer phase
 $(addsuffix _out.xxx,$(OUTPUTS)): cpsd1.bin cpsd2.bin cpsd3.bin cpsd4.bin
-	croposp --DEBUG --verbose --overwrite --itype bin --trim \
+	croposp $(CPSD_OPT) --verbose --overwrite --itype bin --trim \
 	  $(CPSD_LOG) \
 	  --pattern="%I" \
 	  --nsegments=$(CPSD_NSEGMENTS) \
@@ -260,6 +263,19 @@ $(addsuffix _out.xxx,$(OUTPUTS)): cpsd1.bin cpsd2.bin cpsd3.bin cpsd4.bin
 	  --phase phase_out.xxx \
 	  $^
 
+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 \
+	  -nsegments=20 junk.xxx $^
+
+psd%foutra.xxx: psd%.001.asc
+	printf "# #1: %s\n" $< > $@
+	cat $< >> $@
+
+PSD.pdf: psd_out.xxx psd1foutra.xxx psd2foutra.xxx \
+  psd3foutra.xxx psd4foutra.xxx
+	$(CROPOSPLOT) -v --grid -o $@ $^
+
 %.pdf: %_out.xxx
 	$(CROPOSPLOT) -v --grid -o $@ $<