include NorneMakefile
D=Data/
DD=$(datapath)
DD2=$(datapath2)
M=$B/match4.x  
debug1: $(M)
	$(M) base=base.H moni=moni.H nodes=1  irepeat=10 gate=80 nhalf_trk=0  nhalf_bin=1 verb=y  time_u=0. time_d=0 time_r=200 filt_len=500 hori=imap1.txt hori2=imap3.txt optimize=y\
	len_var=.5 r_var=0.5 u_var=0.5 d_var=0.5 perc_var=.95 niter=15  \
	time_r=200 time_r1=160. time_r2=240. \
	time_u1=0   time_u2=100. time_u=200  \
	time_d1=100 time_d2=0.   time_d=200 \
	npop=15 fixed=n lrepeat=y  allnodes=y \
	filter_smp=$@.filter_smp  filter_opt=$@.filter_opt \
	moni_smp=$@.moni_smp        moni_opt=$@.moni_opt \
	rrr_smp=$@.rrr_smp           rrr_opt=$@.rrr_opt    rrr_org=$@.rrr_org \
	dif_smp=$@.dif_smp           dif_opt=$@.dif_opt    dif_org=$@.dif_org \
	 smpfit=$@.smpfit             optfit=$@.optfit      orgfit=$@.orgfit\
	 smp_up=./imap1.txt  \
	 smp_dwn=./imap3.txt \
	 opt_up1=./imap1.txt \
	 opt_dwn1=./imap3.txt\
	 opt_up2=./imap4.txt \
	 opt_dwn2=./imap5.txt
	 mv $@.filter_smp $@.moni_smp $@.rrr_smp $@.rrr_org $@.dif_smp $@.dif_org $@.smpfit $@.orgfit $D/
	 mv $@.filter_opt $@.moni_opt $@.rrr_opt            $@.dif_opt                      $@.optfit $D/


one = ../Oneway/SData
D = Data


Get-%:
	Window3d < $(one)/chirp-$*-1-s.H | Transp > $D/$*-1-s.H
	Window3d < $(one)/chirp-$*-2-s.H | Transp > $D/$*-2-s.H
	Window3d < $(one)/chirp-$*-3-s.H | Transp > $D/$*-3-s.H
	Window3d < $(one)/chirp-$*-4-s.H | Transp > $D/$*-4-s.H
	Window3d < $(one)/chirp-$*-5-s.H | Transp > $D/$*-5-s.H
	touch $@

get-%:
	Window3d < $(one)/shift-$*-1-s.H | Transp > $D/$*-1-s.H
	Window3d < $(one)/shift-$*-2-s.H | Transp > $D/$*-2-s.H
	Window3d < $(one)/shift-$*-3-s.H | Transp > $D/$*-3-s.H
	Window3d < $(one)/shift-$*-4-s.H | Transp > $D/$*-4-s.H
	Window3d < $(one)/shift-$*-5-s.H | Transp > $D/$*-5-s.H
	touch $@

med-get-%:
	Window3d < $(one)/shift-$*-1-s.H | Transp | Median.x wind1=0 wind2=2> $D/med-$*-1-s.H
	Window3d < $(one)/shift-$*-2-s.H | Transp | Median.x wind1=0 wind2=2> $D/med-$*-2-s.H
	Window3d < $(one)/shift-$*-3-s.H | Transp | Median.x wind1=0 wind2=2> $D/med-$*-3-s.H
	Window3d < $(one)/shift-$*-4-s.H | Transp | Median.x wind1=0 wind2=2> $D/med-$*-4-s.H
	Window3d < $(one)/shift-$*-5-s.H | Transp | Median.x wind1=0 wind2=2> $D/med-$*-5-s.H
	touch $@

match-mig-51-5mins-%: $(M) get-mig-51-5mins 
	$(M) base=$D/mig-51-5mins-1-s.H moni=$D/mig-51-5mins-$*-s.H $(match-par)
match-mig-51-rnd-5mins-%: $(M) get-mig-51-rnd-5mins 
	$(M) base=$D/mig-51-rnd-5mins-1-s.H moni=$D/mig-51-rnd-5mins-$*-s.H $(match-par)
match-mig-nocode-%: $(M) get-mig-nocode  
	$(M) base=$D/mig-nocode-1-s.H moni=$D/mig-nocode-$*-s.H $(match-par)
Match-mig-noise-02-5mins-%: $(M) get-mig-noise-02-5mins 
	$(M) base=$D/mig-noise-02-5mins-1-s.H moni=$D/mig-noise-02-5mins-$*-s.H $(match-par)
Match-mig-noise-05-5mins-%: $(M) get-mig-noise-05-5mins 
	$(M) base=$D/mig-noise-05-5mins-1-s.H moni=$D/mig-noise-05-5mins-$*-s.H $(match-par)
Match-mig-noise-10-5mins-%: $(M) get-mig-noise-10-5mins 
	$(M) base=$D/mig-noise-10-5mins-1-s.H moni=$D/mig-noise-10-5mins-$*-s.H $(match-par)
Match-mig-noise-10-6secs-%: $(M) get-mig-noise-10-6secs 
	$(M) base=$D/mig-noise-10-6secs-1-s.H moni=$D/mig-noise-10-6secs-$*-s.H $(match-par)
Match-mig-noise-05-6secs-%: $(M) Get-mig-noise-05-6secs 
	$(M) base=$D/mig-noise-05-6secs-1-s.H moni=$D/mig-noise-05-6secs-$*-s.H $(match-par)
Match-mig-noise-05-36secs-%: $(M) Get-mig-noise-05-36secs 
	$(M) base=$D/mig-noise-05-36secs-1-s.H moni=$D/mig-noise-05-36secs-$*-s.H $(match-par)
Match-mig-noise-05-72secs-%: $(M) Get-mig-noise-05-72secs 
	$(M) base=$D/mig-noise-05-72secs-1-s.H moni=$D/mig-noise-05-72secs-$*-s.H $(match-par)
Match-mig-noise-05-144secs-%: $(M) Get-mig-noise-05-144secs 
	$(M) base=$D/mig-noise-05-144secs-1-s.H moni=$D/mig-noise-05-144secs-$*-s.H $(match-par)
Match-mig-noise-05-288secs-%: $(M) Get-mig-noise-05-288secs 
	$(M) base=$D/mig-noise-05-288secs-1-s.H moni=$D/mig-noise-05-288secs-$*-s.H $(match-par)


match-05-%:
	make Match-mig-noise-05-$*-2 &
	make Match-mig-noise-05-$*-3 & 
	make Match-mig-noise-05-$*-4 & 
	make Match-mig-noise-05-$*-5 & 

match-par= \
        nodes=1  irepeat=10  nhalf_trk=0  nhalf_bin=1 verb=y  time_u=0. time_d=0  \
	hori=imap1.txt hori2=imap3.txt optimize=y \
	len_var=.4 ramp_var=0.4 up_var=0.4 dwn_var=0.4 perc_var=.95 niter=20  \
	time_r=200   time_r1=160. time_r2=240. \
	time_u1=-16  time_u2=-100. time_u=0     \
	time_d1=100  time_d2=0.    time_d=0     \
	npop=15 fixed=n lrepeat=n  allnodes=y \
	percent1=10e-8 percent2=10e-4 percent=10e-6 \
	npop=15 fixed=n lrepeat=y  allnodes=n nodes=16 \
	filter_smp=$@-filter-smp.H  filter_opt=$@-filter-opt.H  \
	moni_smp=$@-moni-smp.H       moni_opt=$@-moni-opt.H 	\
	rrr_smp=$@-rrr-smp.H           rrr_opt=$@-rrr-opt.H    rrr_org=$@-rrr-org.H \
	dif_smp=$@-dif-smp.H           dif_opt=$@-dif-opt.H    dif_org=$@-dif-org.H \
	smpfit=$@-fit-smp.H             optfit=$@-fit-opt.H      orgfit=$@-fit-org.H \
	fitmap=$@-fit-map.H  \
	smp_up=./imap1.txt  \
	smp_dwn=./imap3.txt \
	opt_up1=./imap1.txt \
	opt_dwn1=./imap3.txt\
	opt_up2=./imap4.txt \
	opt_dwn2=./imap5.txt\
	wind1=5  wind2=5  \
	rect1=5  rect2=5  \
	median=y gauss=y \
	median=y gauss=y \
	irepeat=3 niter=5 npop=5 nrepeat_out=3 nodes=16 ; \
	mv $@-filter-smp.H $@-moni-smp.H $@-rrr-smp.H $@-rrr-org.H $@-dif-smp.H $@-dif-org.H $@-fit-smp.H $@-fit-org.H             $D/ ; \
	mv $@-filter-opt.H $@-moni-opt.H $@-rrr-opt.H              $@-dif-opt.H                           $@-fit-opt.H $@-fit-map.H $D/ ; \
	touch $@ 




match-med-mig-51-5mins-%: $(M) med-get-mig-51-5mins 
	$(M) base=$D/med-mig-51-5mins-1-s.H moni=$D/med-mig-51-5mins-$*-s.H nodes=1  irepeat=10 gate=80 nhalf_trk=0  nhalf_bin=1 verb=y  time_u=0. time_d=0 time_r=200 filt_len=500 \
	hori=imap1.txt hori2=imap3.txt optimize=y \
	len_var=.4 ramp_var=0.4 up_var=0.4 dwn_var=0.4 perc_var=.95 niter=15  \
	time_r=200  time_r1=160. time_r2=240. \
	time_u1=-16  time_u2=-24. time_u=0    \
	time_d1=100  time_d2=0.    time_d=0   \
	npop=15 fixed=n lrepeat=n  allnodes=y \
	filter_smp=$@-filter-smp.H  filter_opt=$@-filter-opt.H  \
	moni_smp=$@-moni-smp.H       moni_opt=$@-moni-opt.H 	\
	rrr_smp=$@-rrr-smp.H           rrr_opt=$@-rrr-opt.H    rrr_org=$@-rrr-org.H \
	dif_smp=$@-dif-smp.H           dif_opt=$@-dif-opt.H    dif_org=$@-dif-org.H \
	smpfit=$@-fit-smp.H             optfit=$@-fit-opt.H      orgfit=$@-fit-org.H \
	fitmap=$@-fit-map.H  \
	smp_up=./imap1.txt  \
	smp_dwn=./imap3.txt \
	opt_up1=./imap1.txt \
	opt_dwn1=./imap3.txt\
	opt_up2=./imap4.txt \
	opt_dwn2=./imap5.txt\
	wind1=9  wind2=9  \
	rect1=5  rect2=5  \
	median=y gauss=y \
	irepeat=5 npop=10 
	mv $@-filter-smp.H $@-moni-smp.H $@-rrr-smp.H $@-rrr-org.H $@-dif-smp.H $@-dif-org.H $@-fit-smp.H $@-fit-org.H             $D/
	mv $@-filter-opt.H $@-moni-opt.H $@-rrr-opt.H              $@-dif-opt.H                           $@-fit-opt.H $@-fit-map.H $D/
	touch $@ 



debug-%: $(M)
	$(M) base=base.H moni=moni-$*.H nodes=1  irepeat=10 gate=80 nhalf_trk=0  nhalf_bin=1 verb=y  time_u=0. time_d=0 time_r=200 filt_len=500 hori=imap1.txt hori2=imap3.txt optimize=y\
	len_var=.4 ramp_var=0.4 up_var=0.4 dwn_var=0.4 perc_var=.95 niter=15  \
	time_r=200  time_r1=160. time_r2=240. \
	time_u1=16  time_u2=100. time_u=0     \
	time_d1=100 time_d2=0.   time_d=0     \
	npop=15 fixed=n lrepeat=y  allnodes=y \
	filter_smp=$@.filter_smp.H  filter_opt=$@.filter_opt.H  \
	moni_smp=$@.moni_smp.H       moni_opt=$@.moni_opt.H \
	rrr_smp=$@.rrr_smp.H           rrr_opt=$@.rrr_opt.H    rrr_org=$@.rrr_org.H \
	dif_smp=$@.dif_smp.H           dif_opt=$@.dif_opt.H    dif_org=$@.dif_org.H \
	smpfit=$@.smpfit.H             optfit=$@.optfit.H      orgfit=$@.orgfit.H \
	fitmap=$@.fitmap.H  \
	smp_up=./imap1.txt  \
	smp_dwn=./imap3.txt \
	opt_up1=./imap1.txt \
	opt_dwn1=./imap3.txt\
	opt_up2=./imap4.txt \
	opt_dwn2=./imap5.txt\
	wind1=5  wind2=5  \
	rect1=10 rect2=10 \
	irepeat=5 npop=15 
	mv $@.filter_smp.H $@.moni_smp.H $@.rrr_smp.H $@.rrr_org.H $@.dif_smp.H $@.dif_org.H $@.smpfit.H $@.orgfit.H             $D/
	mv $@.filter_opt.H $@.moni_opt.H $@.rrr_opt.H              $@.dif_opt.H                          $@.optfit.H $@.fitmap.H $D/
	touch $@ 
        



debug_new: $(M)
	#Window3d j1=1 j2=1 j3=1 min3= 1801 n3= 10 < ../Base_near_psdm_cgg_crop.rfc.vt.H $(DD2) > base.H
	#Window3d j1=1 j2=1 j3=1 min3= 1801 n3= 10 < ../Moni_near_psdm_cgg_crop.rfc.vt.H $(DD2) > moni.H
	#$(M) base=base.H moni=moni.H nodes=1  irepeat=7 gate=80 nhalf_trk=1  nhalf_bin=1 verb=y $(DD2)  time_u=2000. time_d=0 time_r=200 filt_len=300 hori=match_filter_event_smooth2.dat hori2=gboyega.dat optimize=y 
	$(M) base=../base.para.H moni=../moni.para.H nodes=1  irepeat=8 gate=80 nhalf_trk=1  nhalf_bin=1 verb=y $(DD2)  time_u=0. time_d=0 time_r=200 filt_len=300 hori=match_filter_event_smooth2.dat hori2=gboyega.dat optimize=y\
        len_var=.5 r_var=0.5 u_var=0.5 d_var=0.5 perc_var=.95 niter=15  \
	npop=32 fixed=n lrepeat=y  allnodes=y \
	filter_smp=$@.filter_smp  filter_opt=$@.filter_opt \
	moni_smp=$@.moni_smp        moni_opt=$@.moni_opt \
	rrr_smp=$@.rrr_smp           rrr_opt=$@.rrr_opt    rrr_org=$@.rrr_org \
	dif_smp=$@.dif_smp           dif_opt=$@.dif_opt    dif_org=$@.dif_org \
	 smpfit=$@.smpfit             optfit=$@.optfit      orgfit=$@.orgfit\
         smp_up=../imap2.dat  \
         smp_dwn=../imap3.dat \
         opt_up1=../imap1.dat \
         opt_dwn1=../imap3.dat\
         opt_up2=../imap4.dat \
         opt_dwn2=../imap5.dat
	 mv $@.filter_smp $@.moni_smp $@.rrr_smp $@.rrr_org $@.dif_smp $@.dif_org $@.smpfit $@.orgfit $D/
	 mv $@.filter_opt $@.moni_opt $@.rrr_opt            $@.dif_opt                      $@.optfit $D/

debug_new2: $(M)
	#Window3d j1=1 j2=1 j3=1 min3= 1801 n3= 10 < ../Base_near_psdm_cgg_crop.rfc.vt.H $(DD2) > base.H
	#Window3d j1=1 j2=1 j3=1 min3= 1801 n3= 10 < ../Moni_near_psdm_cgg_crop.rfc.vt.H $(DD2) > moni.H
	#(M) base=base.H moni=moni.H nodes=1  irepeat=7 gate=80 nhalf_trk=1  nhalf_bin=1 verb=y $(DD2)  time_u=2000. time_d=0 time_r=200 filt_len=300 hori=match_filter_event_smooth2.dat hori2=gboyega.dat optimize=y
	$(M) base=base.H moni=moni.H nodes=1  irepeat=8 gate=80 nhalf_trk=1  nhalf_bin=1 verb=y $(DD2)  time_u=0. time_d=0 time_r=200 filt_len=400 hori=match_filter_event_smooth2.dat hori2=gboyega.dat optimize=y\
        len_var=.5 r_var=0.5 u_var=0.5 d_var=0.5 perc_var=.95 niter=15  \
	time_u1=0. time_d1=250 \
	time_filt1=100 time_filt2=500 \
        npop=32 fixed=n lrepeat=y  allnodes=y \
        filter_smp=$@.filter_smp  filter_opt=$@.filter_opt \
        moni_smp=$@.moni_smp        moni_opt=$@.moni_opt \
        rrr_smp=$@.rrr_smp           rrr_opt=$@.rrr_opt    rrr_org=$@.rrr_org \
        dif_smp=$@.dif_smp           dif_opt=$@.dif_opt    dif_org=$@.dif_org \
         smpfit=$@.smpfit             optfit=$@.optfit      orgfit=$@.orgfit\
         smp_up=../imap1.dat  \
         smp_dwn=../imap3.dat \
         opt_up1=../imap1.dat \
         opt_dwn1=../imap3.dat\
         opt_up2=../imap4.dat \
         opt_dwn2=../imap5.da.
	 mv $@.filter_smp $@.moni_smp $@.rrr_smp $@.rrr_org $@.dif_smp $@.dif_org $@.smpfit $@.orgfit $D/
	 mv $@.filter_opt $@.moni_opt $@.rrr_opt            $@.dif_opt                      $@.optfit $D/
debug_new3: $(M)
	#Window3d j1=1 j2=1 j3=1 min3= 1801 n3= 10 < ../Base_near_psdm_cgg_crop.rfc.vt.H $(DD2) > base.H
	#Window3d j1=1 j2=1 j3=1 min3= 1801 n3= 10 < ../Moni_near_psdm_cgg_crop.rfc.vt.H $(DD2) > moni.H
	#$(M) base=base.H moni=moni.H nodes=1  irepeat=7 gate=80 nhalf_trk=1  nhalf_bin=1 verb=y $(DD2)  time_u=2000. time_d=0 time_r=200 filt_len=300 hori=match_filter_event_smooth2.dat hori2=gboyega.dat optimize=y 
	$(M) base=../base.para.H moni=../moni.para.H nodes=1  irepeat=8 gate=80 nhalf_trk=1 nhalf_bin=1 verb=y $(DD2)  time_u=0. time_d=0 time_r=200 filt_len=300 hori=match_filter_event_smooth2.dat hori2=gboyega.dat optimize=y\
        len_var=.5 r_var=0.5 u_var=0.5 d_var=0.5 perc_var=.95 niter=12 \
	npop=24 fixed=n lrepeat=y  allnodes=y \
	time_u1=0.    time_d1=250.\
	time_u2=-100. time_d2=0.  \
	time_filt1=100 time_filt2=400 \
	filter_smp=$@.filter_smp  filter_opt=$@.filter_opt \
	moni_smp=$@.moni_smp        moni_opt=$@.moni_opt \
	rrr_smp=$@.rrr_smp           rrr_opt=$@.rrr_opt    rrr_org=$@.rrr_org \
	dif_smp=$@.dif_smp           dif_opt=$@.dif_opt    dif_org=$@.dif_org \
	 smpfit=$@.smpfit             optfit=$@.optfit      orgfit=$@.orgfit\
         smp_up=../imap1.dat  \
         smp_dwn=../imap3.dat \
         opt_up1=../imap2.dat \
         opt_dwn1=../imap3.dat\
         opt_up2=../imap4.dat \
         opt_dwn2=../imap5.dat
	 mv $@.filter_smp $@.moni_smp $@.rrr_smp $@.rrr_org $@.dif_smp $@.dif_org $@.smpfit $@.orgfit $D/
	 mv $@.filter_opt $@.moni_opt $@.rrr_opt            $@.dif_opt                      $@.optfit $D/

debug_new4: $(M)
	#Window3d j1=1 j2=1 j3=1 min3= 2158 n3= 9 < ../Base_near_psdm_cgg_crop.rfc.vt.H $(DD2) > base.H
	#Window3d j1=1 j2=1 j3=1 min3= 2158 n3= 9 < ../Moni_near_psdm_cgg_crop.rfc.vt.H $(DD2) > moni.H
	#$(M) base=base.H moni=moni.H nodes=1  irepeat=7 gate=80 nhalf_trk=1  nhalf_bin=1 verb=y $(DD2)  time_u=2000. time_d=0 time_r=200 filt_len=300 hori=match_filter_event_smooth2.dat hori2=gboyega.dat optimize=y 
	$(M) base=base.H moni=moni.H nodes=1  irepeat=7 gate=80 nhalf_trk=1 nhalf_bin=1 verb=y $(DD2)  time_u=0. time_d=0 time_r=200 filt_len=300 hori=match_filter_event_smooth2.dat hori2=gboyega.dat optimize=y\
        len_var=.5 r_var=0.5 u_var=0.5 d_var=0.5 perc_var=.95 niter=25\
	npop=14 fixed=n lrepeat=y  allnodes=y \
	time_u1=-250. time_d1=250.\
	time_u2=-250. time_d2=0.  \
	time_filt1=100 time_filt2=400 \
	filter_smp=$@.filter_smp  filter_opt=$@.filter_opt \
	moni_smp=$@.moni_smp        moni_opt=$@.moni_opt \
	rrr_smp=$@.rrr_smp           rrr_opt=$@.rrr_opt    rrr_org=$@.rrr_org \
	dif_smp=$@.dif_smp           dif_opt=$@.dif_opt    dif_org=$@.dif_org \
	 smpfit=$@.smpfit             optfit=$@.optfit      orgfit=$@.orgfit\
         smp_up=../imap2.dat  \
         smp_dwn=../imap3.dat \
         opt_up1=../imap1.dat \
         opt_dwn1=../imap3.5.dat\
         opt_up2=../imap4.5.dat \
         opt_dwn2=../imap4.dat
	 mv $@.filter_smp $@.moni_smp $@.rrr_smp $@.rrr_org $@.dif_smp $@.dif_org $@.smpfit $@.orgfit $D/
	 mv $@.filter_opt $@.moni_opt $@.rrr_opt            $@.dif_opt                      $@.optfit $D/

1:
	touch 1
#debug%vt:1
#	$(S2V2) $(DD) vtfile=$@ sepfile=debug$*H clipmin=-10 clipmax=10 
#	mv  $@ lvt/

%rrr_smp.vt: %rrr_smp ; Window3d n3=554 f3=0   < $*rrr_smp > 1; Window3d n3=    f3=554 < $*rrr_smp > 2; cp ../match_near_psdm_cgg_crop.rrr2.vt $@ ; $(S2V) $(DD2) vtfile=$@ sepfile=1 sepfile2=2 clipmin=-100000 clipmax=100000 break=y
%rrr_opt.vt: %rrr_opt ; Window3d n3=554 f3=0   < $*rrr_opt > 1; Window3d n3=    f3=554 < $*rrr_opt > 2; cp ../match_near_psdm_cgg_crop.rrr2.vt $@ ; $(S2V) $(DD2) vtfile=$@ sepfile=1 sepfile2=2 clipmin=-100000 clipmax=100000 break=y
%rrr_org.vt: %rrr_org ; Window3d n3=554 f3=0   < $*rrr_org > 1; Window3d n3=    f3=554 < $*rrr_org > 2; cp ../match_near_psdm_cgg_crop.rrr2.vt $@ ; $(S2V) $(DD2) vtfile=$@ sepfile=1 sepfile2=2 clipmin=-100000 clipmax=100000 break=y
%fit.vt: %fit; Reverse which=3 < $< > a; Window3d n3=554 f3=0 <a > 1; Window3d f3=554 <a> 2; cp ../oneslice.vt $@ ; $(S2V) $(DD2) vtfile=$@ sepfile=1 sepfile2=2 clipmin=-100000 clipmax=100000 break=y
%.vt: %; Window3d n3=554 f3=0   < $* > $*1; Window3d n3=    f3=554 < $* > $*2; cp ../match_near_psdm_cgg_crop.dif2.vt $@ ; $(S2V) $(DD2) vtfile=$@ sepfile=$*1 sepfile2=$*2 clipmin=-100000 clipmax=100000 break=y;Rm $*1 $*2





oney = ../Oneway/Yemi/
get-yemi:
	Window3d min1=3 max1=12 min3=0 max3=4.5 < $(oney)/smooth_img.survey_1.line_19-t.H  | Transp > yemi-match-base.H; mv yemi-match-base.H $(oney)/
	Window3d min1=3 max1=12 min3=0 max3=4.5 < $(oney)/smooth_img.survey_11.line_19-t.H | Transp > yemi-match-moni-raw.H; mv yemi-match-moni-raw.H $(oney)/
	Window3d min1=3 max1=12 min3=0 max3=4.5 < $(oney)/smooth_img.survey_11.line_19-s.H | Transp > yemi-match-moni.H; mv yemi-match-moni.H $(oney)/

#=====
yemi-match: $(M) get-yemi
	$(M) base=$(oney)/yemi-match-base.H moni=$(oney)/yemi-match-moni.H $(match-yemi-par)
	Cat3d Data/yemi-match-dif-org.H Data/yemi-match-dif-opt.H | Window3d min2=5 max2=9 | Grey color=j pclip=100| Tube


match-yemi-par= \
        nodes=15  irepeat=10  nhalf_trk=0  nhalf_bin=1 verb=y  time_u=0. time_d=0  \
	hori=imap1.txt hori2=imap3.txt optimize=y \
	len_var=.4 ramp_var=0.4 up_var=0.4 dwn_var=0.4 perc_var=.95 niter=20  \
	gate=0.064 \
	time_filt=.2   time_filt1=.2  time_filt2=.25 \
	time_r=.05   time_r1=.02 time_r2=.05 \
	time_u1=-.05   time_u2=-.01 \
	time_d1=.05    time_d2=0.  \
	npop=15 fixed=n lrepeat=n  allnodes=y \
	percent1=10e-8 percent2=10e-4 percent=10e-6 \
	npop=15 fixed=n lrepeat=y  allnodes=n nodes=16 \
	filter_smp=$@-filter-smp.H  filter_opt=$@-filter-opt.H  \
	moni_smp=$@-moni-smp.H       moni_opt=$@-moni-opt.H 	\
	rrr_smp=$@-rrr-smp.H           rrr_opt=$@-rrr-opt.H    rrr_org=$@-rrr-org.H \
	dif_smp=$@-dif-smp.H           dif_opt=$@-dif-opt.H    dif_org=$@-dif-org.H \
	smpfit=$@-fit-smp.H             optfit=$@-fit-opt.H      orgfit=$@-fit-org.H \
	fitmap=$@-fit-map.H  \
	smp_up=./reflector_1.txt  \
	smp_dwn=./reflector_4.txt \
	opt_up1=./reflector_1.txt \
	opt_dwn1=./reflector_4.txt\
	opt_up2=./reflector_5.txt \
	opt_dwn2=./reflector_6.txt\
	wind1=15  wind2=10  \
	rect1=15  rect2=10  \
	median=y gauss=y \
	irepeat=5 npop=10 \
	median=y gauss=y \
	fact=2 gap=300 lgap=y \
	irepeat=5 niter=10 npop=10 nrepeat_out=2 nodes=16 verb=10 ; \
	mv $@-filter-smp.H $@-moni-smp.H $@-rrr-smp.H $@-rrr-org.H $@-dif-smp.H $@-dif-org.H $@-fit-smp.H $@-fit-org.H             $D/ ; \
	mv $@-filter-opt.H $@-moni-opt.H $@-rrr-opt.H              $@-dif-opt.H                           $@-fit-opt.H $@-fit-map.H $D/  ;  \
	touch $@ 


