#This script creates a waterfall plot of the time-based signal
##############################################################################################################
###
#What to do before?
#Extract sound selection (preserve times)
#Extract visible pulses
#for one column the number of periods should not be over 55
#otherwise step has to be ajusted 
#in this script the first period defines the horizontal edges, all other orient at it but dur will alternate
###############################################################################################################
fenster=1.9
step=0.17

Erase all

select PointProcess untitled

anzahlpunkte=Get number of points
anzahlperioden=anzahlpunkte-1
##########zahl vor / gibt an, bis wohin 
fraktion=12/anzahlperioden

pt1=Get time from index... 2
pt1$=fixed$ (pt1, 5)
pt12=Get time from index... 3
pt12$=fixed$ (pt12, 5)
pt2=Get time from index... 'anzahlpunkte'
pt2$=fixed$ (pt2, 5)

d1='pt12'-'pt1'

for i from 2 to anzahlpunkte-1
select PointProcess untitled
periodenanfang=Get time from index... i

periodenende=Get time from index... i+1

d2='periodenende'-'periodenanfang'
gf='d2' / 'd1'
x=6*'gf'

#echo 'x'
#pause ok?

x$=fixed$ (x, 2)

select Sound untitled


if i >1

oben=i*'step'-'step'
unten='oben'+'fenster'
else
oben=i*'fenster'-'fenster'
#oben=0
unten='fenster'
endif

#if i >1
Viewport... 0 'x$' 'oben' 'unten'
#else
#Viewport... 0 6 'oben' 'unten'
#endif

Draw... 'periodenanfang' 'periodenende' 0.5 -0.6 no curve
Text left... no 'i'
One mark right... 0 no yes no
One mark left... 0 no yes no
if i=2
Viewport... 0 6 0 12
p1=periodenanfang
p2=periodenanfang+(periodenende-periodenanfang)/4
p3=periodenanfang+(periodenende-periodenanfang)/2
p4=periodenanfang+(periodenende-periodenanfang)*0.75
p5=periodenende

t2=(periodenende-periodenanfang)/4
t2$=fixed$ (t2, 6)

t3=(periodenende-periodenanfang)/2
t3$=fixed$ (t3, 6)
t4=(periodenende-periodenanfang)*0.75
t4$=fixed$ (t4, 6)
t5=p5-periodenanfang
t5$=fixed$ (t5, 6)

One mark top... 'p1' no yes yes
One mark top... 'p2' no yes yes
One mark top... 'p3' no yes yes
One mark top... 'p4' no yes yes
One mark top... 'p5' no yes yes
One mark top... 'p1' no no no 0
One mark top... 'p2' no no no 't2$'
One mark top... 'p3' no no no 't3$'
One mark top... 'p4' no no no 't4$'
One mark top... 'p5' no no no 't5$'
endif
endfor


Viewport... 0 6 0 12
#Draw inner box
Text bottom... no 'anzahlperioden' cycles of selection from 'pt1$' to 'pt2$' sec
#Text left... yes VF_VFT F0/2 lost of 'subharmonic coupling'
Text left... yes short lowering of the fundamental in kargyraa2 (AEF-VF)

Viewport... 6 12 0 4
Draw... 0 0 0 0 yes curve
select PointProcess untitled
Draw... 0 0 no

select Sound untitled
To Spectrogram... 0.009 2500 0.002 20 Bartlett (triangular)
Viewport... 6 12 4 8
Paint... 0 0 0 0 100 yes 50 6 0 yes
Text top... no Bartlett (triangular)
select Sound untitled
To Spectrogram... 0.045 2500 0.001 20 Gaussian
Viewport... 6 12 8 12
Paint... 0 0 0 0 100 yes 50 6 0 yes
Text top... no Gaussian
select PointProcess untitled
Draw... 0 0 no

#################Sven Grawunder, November 2004, MPI EVA Leipzig######################################