Commit 859ffe34 authored by thomas.forbriger's avatar thomas.forbriger Committed by thomas.forbriger
Browse files

boxcar

This is a legacy commit from before 2015-03-01.
It may be incomplete as well as inconsistent.
See COPYING.legacy and README.history for details.


SVN Path:     http://gpitrsvn.gpi.uni-karlsruhe.de/repos/TFSoftware/trunk
SVN Revision: 524
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 5dbf97ab
c this is <siggen.f>
c------------------------------------------------------------------------------
c ($Id: siggen.f,v 1.6 2001-10-25 12:03:11 forbrig Exp $)
c ($Id: siggen.f,v 1.7 2001-11-02 15:14:05 forbrig Exp $)
c
c 24/10/2001 by Thomas Forbriger (IMGF Frankfurt)
c
......@@ -8,15 +8,16 @@ c SIGnal GENerator
c
c REVISIONS and CHANGES
c 24/10/2001 V1.0 Thomas Forbriger
c 02/11/2001 V1.1 boxcar
c
c==============================================================================
c
program siggen
c
character*(*) version
parameter(version='SIGGEN V1.0 SIGnal GENerator')
parameter(version='SIGGEN V1.1 SIGnal GENerator')
character*(*) SIGGEN_CVS_ID
parameter(SIGGEN_CVS_ID='$Id: siggen.f,v 1.6 2001-10-25 12:03:11 forbrig Exp $')
parameter(SIGGEN_CVS_ID='$Id: siggen.f,v 1.7 2001-11-02 15:14:05 forbrig Exp $')
c
c parameters
integer nsig,ncyc
......@@ -177,6 +178,11 @@ c
print *,' single zero-crossing at Tm and terminates at Te'
print *,' with value 0 and frequency 0Hz.'
print *,' '
print *,'8 boxcar'
print *,' '
print *,' All samples are zero. Only ''n'' samples after'
print *,' ''Ta'' seconds will be ''a''.'
print *,' '
print *,SIGGEN_CVS_ID
stop
endif
......@@ -407,6 +413,32 @@ c
data(i)=0.
endif
enddo
elseif (nsig.eq.8) then
if (verbose) then
print *,'boxcar'
print *,'parameters:'
print 50,'amplitude','a',a,' '
print 50,'sampling interval','d',d,'s'
print 50,'time series lentgh','T',t,'s'
print 50,'onset of wavelet','Ta',ta,'s'
print 50,'number of non-zero samples','n',ncyc,' '
endif
nsamples=int(t/d)
if (verbose) then
print *,'derived parameters:'
print 52,'number of samples',nsamples,' '
endif
do i=1,nsamples
ti=d*float(i-1)
if (ti.lt.ta) then
data(i)=0.
elseif (ncyc.gt.0) then
data(i)=a
ncyc=ncyc-1
else
data(i)=0.
endif
enddo
else
stop 'ERROR: unknown signal'
endif
......
Supports Markdown
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