# compute tdm data for madison MOH90

CFLAGS= -g -I..
PGMS = mkfig2 mkfig3 mkfig4 mkdm4
TSRH = madp5a.tsr madp5b.tsr madp5c.tsr madp5d.tsr
TSR1 = mad1a.tsr mad1b.tsr mad1c.tsr mad1d.tsr
TSR2 = mad2a.tsr mad2b.tsr mad2c.tsr mad2d.tsr
TSR4 = mad4a.tsr mad4b.tsr mad4c.tsr mad4d.tsr
TSRS =  mad0.tsr $(TSRH) $(TSR1) $(TSR2) $(TSR4)
FIG1 = fig1a1.dat fig1a2.dat fig1a3.dat fig1a4.dat \
       fig1b1.dat fig1b2.dat fig1b3.dat fig1b4.dat 
FIG2 = fig2a1.dat fig2a2.dat fig2a3.dat fig2a4.dat fig2a5.dat
FIG3 = fig3a1.dat fig3a2.dat fig3a3.dat fig3a4.dat fig3a5.dat fig3a6.dat
FIG4 = fig4a1.dat fig4a2.dat
DM1 = dm1a.dat dm1b.dat
DM2 = dm2a.dat dm2b.dat
DM3 = dm3a1.dat dm3a2.dat dm3a3.dat dm3a4.dat dm3a5.dat \
       dm3b1.dat dm3b2.dat dm3b3.dat dm3b4.dat dm3b5.dat
DM4 = dm4a0.dat dm4a1.dat dm4a2.dat
DM10 = dm10a1.dat dm10a2.dat dm10a3.dat dm10a4.dat dm10a5.dat \
        dm10b1.dat dm10b2.dat dm10b3.dat dm10b4.dat dm10b5.dat
DM13 = dm13a.dat dm13b.dat dm13c.dat dm13d.dat

all : figs dm

figs: fig1 fig2 fig3 fig4 

dm: dm1 dm2 dm3 dm4 dm10 dm13

fig1 : fig1.plt $(FIG1) 
$(FIG1) : mad0.tsr
	datman fig1

dm1 : dm1.plt $(DM1) 
$(DM1): mad0.tsr
	datman dm1
dm2 : dm2.plt $(DM2) 
$(DM2): mad0.tsr
	datman dm2
dm3 : dm3.plt $(DM3) 
$(DM3): mad0.tsr
	datman dm3
dm4 : dm4.plt $(DM4) 
$(DM4): mkdm4
	mkdm4 mad0.tsr
fig2 : fig2.plt $(FIG2) 
$(FIG2): mkfig2 $(TSRS)
	mkfig2 mad1a.tsr 
fig3 : fig3.plt $(FIG3) 
$(FIG3): mkfig3 $(TSRS)
	mkfig3 mad1a.tsr 
fig4 : fig4.plt fig4a.inc $(FIG4) 
$(FIG4): mkfig4 $(TSRS)
	mkfig4 mad1a.tsr 
dm10 : dm10.plt $(DM10) 
$(DM10): $(TSRS)
	datman dm10
dm13 : dm13.plt $(DM13) 
$(DM13): $(TSRS)
	datman dm13

#
# click
#
mad0.tsr : mad0.par
	tdm -i -f mad0
	reor mad0
#
# f=0.5kHz
#
madp5a.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbrp5.stm stgain=1.897e-10
	reor mad0 madp5a
madp5b.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbrp5.stm stgain=1.897e-9
	reor mad0 madp5b
madp5c.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbrp5.stm stgain=1.897e-8
	reor mad0 madp5c
madp5d.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbrp5.stm stgain=1.897e-7
	reor mad0 madp5d
#
# f=1kHz
#
mad1a.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr1.stm stgain=1.826e-10
	reor mad0 mad1a
mad1b.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr1.stm stgain=1.826e-9 
	reor mad0 mad1b
mad1c.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr1.stm stgain=1.826e-8 
	reor mad0 mad1c
mad1d.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr1.stm stgain=1.826e-7 
	reor mad0 mad1d
#
# f=2kHz
#
mad2a.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr2.stm stgain=2.151e-10
	reor mad0 mad2a
mad2b.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr2.stm stgain=2.151e-9 
	reor mad0 mad2b
mad2c.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr2.stm stgain=2.151e-8 
	reor mad0 mad2c
mad2d.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr2.stm stgain=2.164e-7 
	reor mad0 mad2d
#
# f=4kHz
#
mad4a.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr4.stm stgain=1.182e-9
	reor mad0 mad4a
mad4b.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr4.stm stgain=1.182e-8 
	reor mad0 mad4b
mad4c.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr4.stm stgain=1.182e-7 
	reor mad0 mad4c
mad4d.tsr : mad0.par
	tdm -i -f mad0 stimfn=tnbr4.stm stgain=1.182e-6 
	reor mad0 mad4d

mkfig2 : mkfig2.o
	$(CC) -o mkfig2 mkfig2.o -lm 
mkfig3 : mkfig3.o
	$(CC) -o mkfig3 mkfig3.o -lm 
mkfig4 : mkfig4.o
	$(CC) -o mkfig4 mkfig4.o -lm
mkdm4 : mkdm4.o
	$(CC) -o mkdm4 mkdm4.o -lm

clean :
	rm -f *.tsd *.sta $(PGMS) *.exe *.log

empty :
	rm -f *.tsd *.sta $(PGMS) *.exe *.log *.tsr *.dat

