Commit a666c082 authored by thomas.forbriger's avatar thomas.forbriger
Browse files

croposp [WP]: set raw data for tests

parent 761702bc
......@@ -40,11 +40,9 @@ clean: ;
#
# set time series parameters
# --------------------------
CPSD_NSAMPLESRAW=100000
#CPSD_NSAMPLESRAW=1000
CPSD_DT=0.1
CPSD_DT=0.01
CPSD_DUR=8640.
#CPSD_DUR=864.
CPSD_DUR=864.
# ----------------------------------------------------------------------
# create raw time series
......@@ -54,8 +52,7 @@ CPSD_DUR=8640.
# own specific seed value for the random number generator. We use the
# nanosecond field of the system clock.
#
CPSD_SIGGENOPT=-v -o -n $(CPSD_NSAMPLESRAW) -d $(CPSD_DT) \
-T $(CPSD_DUR) -ot bin
CPSD_SIGGENOPT=-v -o -d $(CPSD_DT) -T $(CPSD_DUR) -ot bin
cpsd1raw.bin:
/bin/bash -c 'export GSL_RNG_SEED=$$(date +"%N"); \
......@@ -94,27 +91,27 @@ rawsig.ps: cpsd1raw.bin cpsd2raw.bin cpsd3raw.bin cpsd4raw.bin cpsdbgraw.bin
# common in all reacordings.
#
# set amplitudes
CPSD_AMP1=1.e-3
CPSD_AMP2=1.e-4
CPSD_AMP3=3.e-3
CPSD_AMP4=5.e-5
CPSD_AMPBG=2.
CPSD_AMP1=0.2
CPSD_AMP2=4.
CPSD_AMP3=100.
CPSD_AMP4=1.e-2
CPSD_AMPBG=10.
# modify white noise to channel specific frequency characteristic
cpsd1mod.bin: cpsd1raw.bin
printf "fac $(CPSD_AMP1)\nlpb 1.,2.\ndif 0\nend\n" \
printf "fac $(CPSD_AMP1)\nhpb 1.,2\nend\n" \
| tidofi -type bin -Type bin -v -o -cs $@ $<
cpsd2mod.bin: cpsd2raw.bin
printf "fac $(CPSD_AMP2)\nint 0\navg 0end\n" \
printf "fac $(CPSD_AMP2)\nhpb 4.,4\nlpb 2.,4\nend\n" \
| tidofi -type bin -Type bin -v -o -cs $@ $<
cpsd3mod.bin: cpsd3raw.bin
printf "fac $(CPSD_AMP3)\nlpb 10.,2\nend\n" \
printf "fac $(CPSD_AMP3)\nlp2 100.,3\nend\n" \
| tidofi -type bin -Type bin -v -o -cs $@ $<
cpsd4mod.bin: cpsd4raw.bin
printf "fac $(CPSD_AMP4)\nhpb 100.0,2\nend\n" \
printf "fac $(CPSD_AMP4)\nend\n" \
| tidofi -type bin -Type bin -v -o -cs $@ $<
cpsdbgmod.bin: cpsdbgraw.bin
printf "fac $(CPSD_AMPBG)\nlpb 2.,4\navg 0\nhpb 200.,2\nend\n" \
printf "fac $(CPSD_AMPBG)\nlp1 2.\navg 0\nhpb 100.,4\nend\n" \
| tidofi -type bin -Type bin -v -o -cs $@ $<
# add common noise signal to channel specific noise floor
......@@ -266,6 +263,7 @@ CPSD_RAWREFFILES=$(addprefix psd,$(addsuffix modfoutra.xxx,$(CPSD_SET) bg))
# background noise separately)
PSDraw.pdf: psdraw.xxx $(CPSD_RAWREFFILES)
$(CROPOSPLOT) --title="$@ \n $(CPSD_CROPOSP_OPT)" \
--ylim=1.e-14:100. \
--titfontsize='small' -v --grid -o $@ $^
# display PSD of apparent seismometer output in comparison with foutra
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment