MAKE_DEPEND=yes
include ${SEPINC}/SEP.top
include Migmakefile

BIN=../Bin
PAR=../Par
FIG=../Fig
DATA=../Data
SIGDIR=/mnt/sepdata/1a/sigsbee_multiple

UF90LIBS=-lsep2df90 -lsepparf90 -lseppar -lsupersetf90 -lsuperset -lsepgeef90
UF90FLAGS=-O1
#UF90FLAGS =-parallel -openmp -openmp_report2  -O3 -axW -vec_report3
MAINLIBS=-lsep2df90 -lsep3df90 -lsep3d -lsepf90 -lsep -L${LIBDIR}  -lsepgeef90
#UF90LIBS=${MPILIBS} ${MAINLIBS} -lsepparf90 -lseppar_mpi -lsepgeef90 -lsep2df90 -static-libcxa -lguide#-Bstatic 
#UF77LIBS = ${BEILIB}


prim01.H: ${DATA}/cmp2667_muum.H ${BIN}/LSQpr.x
	${BIN}/LSQpr.x n1=650 n2=256 n3=1 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 par=${PAR}/muti.P <${DATA}/cmp2667_muum.H >${DATA}/$@

prim01.v: ${DATA}/prim01.H
	Grey <${DATA}/prim01.H >${FIG}/$@ 
	Tube <${FIG}/$@

cmpmut.H: ${DATA}/cmp2667_muum.H ${BIN}/Pmute.x
	${BIN}/Pmute.x n1=650 n2=256 n3=1 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 par=${PAR}/mute.P <${DATA}/cmp2667_muum.H >${DATA}/$@

cmpmut.v: ${DATA}/cmpmut.H
	Grey <${DATA}/cmpmut.H >${FIG}/$@ 
	Tube <${FIG}/$@

##################################
#
# The progs
#
mut:Pmute.x
	mv -f Pmute.x ../Bin/.
 
dst: DPtest.x
	mv -f DPtest.x ../Bin/.

avo:AVO.x
	mv -f AVO.x ../Bin/.

avo1:AVO1.x
	mv -f AVO1.x ../Bin/.

stk: Stack.x
	mv -f Stack.x ../Bin/.

LSQ: LSQpr.x
	mv -f LSQpr.x ../Bin/.

SUBT: Subtract.x
	mv -f Subtract.x ../Bin/.

NORM:norm.x
	mv -f norm.x ../Bin/.

#################################
#
# Now the sorted cigs
#
prim01sei.H: ${DATA}/prim_data1_h_sort.H ${BIN}/LSQpr.x
	${BIN}/LSQpr.x n1=650 n2=256 n3=10 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 par=${PAR}/mutei.P <${DATA}/prim_data1_h_sort.H >${DATA}/$@

prim01si.H: ${DATA}/prim_data1_h_sort.H ${BIN}/LSQpr.x
	${BIN}/LSQpr.x n1=650 n2=256 n3=10 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 par=${PAR}/muti.P <${DATA}/prim_data1_h_sort.H >${DATA}/$@

#######################
#
# Using Solver
#
#######################

xxx.H:
	<${DATA}/prim_data1_h.H Window f3=100 n3=1 >${DATA}/$@

prim_tiny.H: ${DATA}/prim_data2_h.H ${BIN}/AVO.x
	${BIN}/AVO.x par=${PAR}/mutei.P solv=1 <${DATA}/prim_data2_h.H >${DATA}/$@

prim_tiny1.H: ${DATA}/prim_data2_h.H ${BIN}/AVO.x
	${BIN}/AVO.x par=${PAR}/mutei.P solv=1 <${DATA}/prim_data2_h.H A_out=${DATA}/Atiny.H B_out=${DATA}/Btiny.H C_out=${DATA}/Ctiny.H >${DATA}/$@

prim_tiny200.H: ${DATA}/prim_data2_h.H ${BIN}/AVO.x
	${BIN}/AVO.x par=${PAR}/mutei.P solv=1 niter=200 <${DATA}/prim_data2_h.H >${DATA}/$@

prim_reg.H: ${DATA}/xxx.H ${BIN}/AVO.x
	${BIN}/AVO.x par=${PAR}/mutei.P solv=2 <${DATA}/xxx.H >${DATA}/$@

prim_prc_eps0.2.H: ${DATA}/prim_data2_h.H ${BIN}/AVO.x
	${BIN}/AVO.x par=${PAR}/mutei.P solv=3 eps=0. <${DATA}/prim_data2_h.H >${DATA}/$@

prim_prc1.H: ${DATA}/prim_data2_h.H ${BIN}/AVO.x
	${BIN}/AVO.x par=${PAR}/mutei.P solv=3 <${DATA}/prim_data2_h.H A_out=${DATA}/Atiny.H B_out=${DATA}/Btiny.H C_out=${DATA}/Ctiny.H > ${DATA}/$@

prim_new1.H: 
	${BIN}/AVO.x par=${PAR}/mutei.P solv=3 <junk1.H A_out=${DATA}/Anew.H B_out=${DATA}/Bnew.H C_out=${DATA}/Cnew.H > ${DATA}/$@
#	${BIN}/AVO.x par=${PAR}/mutei.P solv=3 <${DATA}/prims.all_2d.scale_mine.H A_out=${DATA}/Anew.H B_out=${DATA}/Bnew.H C_out=${DATA}/Cnew.H > ${DATA}/$@

prim_new3.H:
	${BIN}/AVO.x par=${PAR}/mutei.P solv=3 <${DATA}/prim_data2_h.H A_out=${DATA}/Anew3.H B_out=${DATA}/Bnew3.H C_out=${DATA}/Cnew3.H > ${DATA}/$@
#	${BIN}/AVO.x par=${PAR}/mutei.P solv=3 <${DATA}/prims.all_2d.scale_mine.H A_out=${DATA}/Anew.H B_out=${DATA}/Bnew.H C_out=${DATA}/Cnew.H > ${DATA}/$@

prim_new2.H: 
	${BIN}/AVO1.x par=${PAR}/mutei.P solv=3 <${DATA}/prims.all_2d.scale_mine.H A_out=${DATA}/Anew2.H B_out=${DATA}/Bnew2.H C_out=${DATA}/Cnew2.H > ${DATA}/$@
#
# Now the whole line
#
prim02w.H:
	${BIN}/LSQpr.x par=${PAR}/mutei.P <${DATA}/prims.all_2d.scale_mine.H A_out=${DATA}/Alsq.H B_out=${DATA}/Blsq.H C_out=${DATA}/Clsq.H >${DATA}/$@

prim01w.H: 
	${BIN}/LSQpr.x par=${PAR}/mutei.P <${DATA}/prims.all_2d.scale_mine.H >${DATA}/$@

prim01i.H: ${DATA}/prim_data1_h.H ${BIN}/LSQpr.x
	${BIN}/LSQpr.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=0.0 par=${PAR}/mutei.P <${DATA}/prim_data1_h.H >${DATA}/$@


#prim01.v: ${DATA}/prim01.H
#	Grey <${DATA}/prim01.H >${FIG}/$@ 
#	Tube <${FIG}/$@

cmpmutw.H: ${DATA}/prim_data1_h.H ${BIN}/Pmute.x
	${BIN}/Pmute.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 par=${PAR}/mute.P <${DATA}/prim_data1_h.H >${DATA}/$@

stkcmp01w.H: ${DATA}/cmpmutw.H ${BIN}/Stack.x
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 <${DATA}/cmpmutw.H >${DATA}/$@

stkprim01w.H: ${DATA}/prim01w.H ${BIN}/Stack.x
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 <${DATA}/prim01w.H >${DATA}/$@

stkprim04estw.H: ${DATA}/prim_est04w_f500.H ${BIN}/Stack.x
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 <${DATA}/prim_est04w_f500.H >${DATA}/$@



stkcmp01w_near.H: ${DATA}/cmpmutw.H ${BIN}/Stack.x
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 <${DATA}/cmpmutw.H >${DATA}/$@

stkprim01w_near.H: ${DATA}/prim01w.H ${BIN}/Stack.x
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 <${DATA}/prim01w.H >${DATA}/$@

stkcmp01w_nearmid.H: ${DATA}/cmpmutw.H ${BIN}/Stack.x
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 <${DATA}/cmpmutw.H >${DATA}/$@

stkprim01w_nearmid.H: ${DATA}/prim01w.H ${BIN}/Stack.x
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 <${DATA}/prim01w.H >${DATA}/$@

stkprim01i_near.H: ${DATA}/prim01i.H ${BIN}/Stack.x 
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 par=$(PAR)/ang1.P <${DATA}/prim01w.H >${DATA}/$@

stkprim01i_nearmid.H: ${DATA}/prim01i.H ${BIN}/Stack.x 
	${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 par=$(PAR)/ang2.P <${DATA}/prim01w.H >${DATA}/$@

stkprim01i.H: ${DATA}/prim01i.H ${BIN}/Stack.x ${BIN}/Pmute.x
	${BIN}/Pmute.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 par=$(PAR)/mute.P par=$(PAR)/ang0.P <${DATA}/prim01w.H | ${BIN}/Stack.x n1=650 n2=256 n3=900 d1=6.67 d2=0.46875 d3=26.67  o1=1000.5 o2=-60 o3=2667.0 >${DATA}/$@



#cmpmut.v: ${DATA}/cmpmut.H
#	Grey <${DATA}/cmpmut.H >${FIG}/$@ 
#	Tube <${FIG}/$@


clean:
	rm -f *.o *.fix.* *.mod *.x .make.dependencies.LINUX86_64

include ${SEPINC}/SEP.bottom