*******************************************************************************
*                             SoundTrack, v. 1.00                             *
*                                                                             *
*                                     by                                      *
*                                                                             *
*             L. Padilla (e-mail: padilla at domain "gae ucm es")             *
*                                                                             *
*                            Madrid, January 2003                             *
*                                                                             *
* Latest version at http://www.gae.ucm.es/~padilla/extrawork/soundtrack.kumac *
*******************************************************************************

macro soundtrack

close 0
hi/del 0
opt *
set *
zone 1 2

nt/cre 1 soundtrack1 2 ! ! t s
nt/rea 1 soundtrack1.dat

nt/cre 2 soundtrack2 2 ! ! t s
nt/rea 2 soundtrack2.dat

1d 11 'High density track' 1000 0. 10.
1d 12 'Low density track' 1000 0. 10.

2d 21 'High density track' 500 0. 5. 1000 0. 10.
2d 22 'Low density track' 500 0. 5. 1000 0. 10.

nt/plot 1.s%t ! ! ! ! l
atit 'Time (s)' 'Amplitude'
ms1 = ($HINFO(1000000,ymax) + abs($HINFO(1000000,ymin)))/12.

nt/plot 2.s%t ! ! ! ! l
atit 'Time (s)' 'Amplitude'
ms2 = ($HINFO(1000000,ymax) + abs($HINFO(1000000,ymin)))/12.

mess Thresholds for high density track: [ms1] and low density track: [ms2]

nt/pro 11 1.1000.*soundtrack.f([ms1],60)
nt/pro 12 2.1000.*soundtrack.f([ms2],20)

nt/pro 21 1.1000.*soundtrack.f([ms1],60)%t
nt/pro 22 2.1000.*soundtrack.f([ms2],20)%t

zone 2 2

wait

hi/plot 11(:1.)
atit 'Flux reversal duration (ms)' 'Number of flux reversals'
hi/plot 12(:1.5)
atit 'Flux reversal duration (ms)' 'Number of flux reversals'

hi/plot 21(1.5:1.85,:1.)
atit 'Time (s)' 'Flux reversal duration (ms)'
hi/plot 22(1.7:1.9,:1.5)
atit 'Time (s)' 'Flux reversal duration (ms)'

return