# make figures for MOH93

IDDAT=ida1.dat idb1.dat ida2.dat idb2.dat ida3.dat idb3.dat \
    ida4.dat idb4.dat ida5.dat idb5.dat

figs : fig4 fig5 fig6 pwrflo

fig4 : $(IDDAT) 
fig5 : 4kio0.dat 4kio1.dat 
fig6 : fdpswp.dat

pwrflo : bmvel.dat pwrflo.dat 
tdm.cal : moh93.par
	tdm -C moh93
bmvel.dat pwrflo.dat : tdm.cal moh93.par
	tdm -c moh93 f1=2002 f2=2490 l1=60 l2=55 ifrmt=2 jfrmt=4 ntsgat=10000 ncyc=4
	pwrflo

$(IDDAT) :
	tdm -i moh93
	reor
	datman fig4

4kio0.dat: # passive, linear case (gam=0)
	tdm -i -e 4kio.lst moh93 gam=0 
	mv 4kio.dat 4kio0.dat 
4kio1.dat: # active, non-linear case (gam=1)
	tdm -i -e 4kio.lst moh93 gam=1
	mv 4kio.dat 4kio1.dat 
fdpswp.dat : fdpswp.lst 
	tdm -i -e fdpswp.lst moh93

clean:
	rm -f *.bak *.log tdm.* *.stm

empty:
	rm -f *.bak *.log tdm.* *.stm *.tsr *.dat

