# Make file for Tomography


#############################
##  Artman, February 2004  ##
#############################

# Calculate SNR
Artman_%_SNR.H: Artman_SI_C_%_env.H
	SNRcalc.x < Artman_SI_C_$*_env.H header=Artman_SI_C_$*_env.H@@ par=${P}/pickpars.p unit=0.01 > $@

# Eikonal Tomography
Artman_%_Etom.H: Artman_%_SNR.H
	EikonalTomo_omp.x < Artman_$*_SNR.H  par=${P}/tomogrid.p  par=${P}/statpars.p  > $@
Artman_%_Etom1.H: Artman_%_SNR.H
	EikonalTomo_omp.x < Artman_$*_SNR.H  par=${P}/tomogrid.p  par=${P}/statpars1.p > $@
Artman_%_Etom2.H: Artman_%_SNR.H
	EikonalTomo_omp.x < Artman_$*_SNR.H  par=${P}/tomogrid.p  par=${P}/statpars2.p > $@
Artman_%_Etom3.H: Artman_%_SNR.H
	EikonalTomo_omp.x < Artman_$*_SNR.H  par=${P}/tomogrid.p  par=${P}/statpars3.p > $@
Artman_%_Etom4.H: Artman_%_SNR.H
	EikonalTomo_omp.x < Artman_$*_SNR.H  par=${P}/tomogrid.p  par=${P}/statpars4.p > $@

Artman_V_Etom: Artman_V_Etom1.H Artman_V_Etom2.H Artman_V_Etom3.H Artman_V_Etom4.H

# Mask Map
Artman_%_Mask.H: Artman_%_Etom.H Artman_%_xy_alive.H
	Window < Artman_$*_Etom.H n3=1 | sdrMaskValhall.x value=0 header=Artman_$*_xy_alive.H unit=1000 rmax=125 angmin=180 angrmax=1500 getmask=1 > $@


#############################
##  Jianhua, January 2005  ##
#############################

# Compute S/N and pick
Jianhua_PICKPARS = v=350 dtau=2
Jianhua_STATPARS = snr=2 minr=0 maxr=1750

# Calculate SNR
Jianhua_V_SNR.H: Jianhua_SI_C_V_env.H
	SNRcalc.x < Jianhua_SI_C_V_env.H header=Jianhua_SI_C_V_env.H@@ par=${P}/pickpars.p unit=0.01 > $@
Jianhua_V3_SNR.H: Jianhua_SI3_C_V_env.H
	SNRcalc.x < Jianhua_SI3_C_V_env.H header=Jianhua_SI3_C_V_env.H@@ par=${P}/pickpars.p unit=0.01 > $@

# Eikonal Tomography
Jianhua_V_Etom.H: Jianhua_V_SNR.H
	EikonalTomo_omp.x < Jianhua_V_SNR.H   par=${P}/tomogrid.p par=${P}/statpars.p > $@
Jianhua_V3_Etom.H: Jianhua_V3_SNR.H
	EikonalTomo_omp.x < Jianhua_V3_SNR.H   par=${P}/tomogrid.p par=${P}/statpars.p > $@

# Mask Map
Jianhua_V_Mask.H: Jianhua_V_Etom.H Jianhua_V_xy_alive.H
	Window < Jianhua_V_Etom.H n3=1 | sdrMaskValhall.x value=0 header=Jianhua_V_xy_alive.H unit=1000 rmax=125 angmin=180 angrmax=1500 getmask=1 > $@
Jianhua_V3_Mask.H: Jianhua_V3_Etom.H Jianhua_V_xy_alive.H
	Window < Jianhua_V3_Etom.H n3=1 | sdrMaskValhall.x value=0 header=Jianhua_V_xy_alive.H unit=1000 rmax=125 angmin=180 angrmax=1500 getmask=1 > $@

###########################
##  Laura, October 2008  ##
###########################

# Compute S/N and pick
Laura_PICKPARS = v=350 dtau=2
Laura_STATPARS = snr=2 minr=0 maxr=1750

# Calculate SNR
Laura_V_SNR.H: Laura_SI_C_V_env.H
	SNRcalc.x < $< header=$<@@ par=${P}/pickpars.p  > $@

# Eikonal Tomography
Laura_V_Etom.H: Laura_V_SNR.H
	EikonalTomo_omp.x < Laura_V_SNR.H  par=${P}/tomogrid.p par=${P}/statpars.p > $@ estimates=estimates.H velocity=1

# Mask Map
Laura_V_Mask.H: Laura_V_Etom.H Laura_V_xy_alive.H
	Window < Laura_V_Etom.H n3=1 | sdrMaskValhall.x value=0 header=Laura_V_xy_alive.H unit=1000 rmax=125 angmin=180 angrmax=1500 getmask=1 > $@


##############################
##  Example shot for Laura  ##
##############################

# Calculate SNR
Laura_SHOT_V_SNR.H: Laura_SHOT_C_V_env.H
	SNRcalc.x < $< header=$<@@ par=${P}/pickpars.p  > $@

# Eikonal Tomography
Laura_SHOT_V_Etom.H Laura_SHOT_V_tsurf.H: Laura_SHOT_V_SNR.H
	EikonalTomo_omp.x < $<  par=${P}/tomogrid.p par=${P}/statpars.p tsurf=Laura_SHOT_V_tsurf.H > $@


