From 253bd10a3a4114be6b448d2d911be5dc31a0ad61 Mon Sep 17 00:00:00 2001
From: Pawan Bharadwaj <bharadwaj.pawan@gmail.com>
Date: Thu, 23 Aug 2018 13:03:11 -0400
Subject: [PATCH] added downhole gallery problem; testing multicomponent
 inversion

---
 src/Data/misfits.jl   |  1 -
 src/Gallery/fwi.jl    | 18 +++++++++++++++++-
 test/Fdtd/backprop.jl |  2 +-
 3 files changed, 18 insertions(+), 3 deletions(-)

diff --git a/src/Data/misfits.jl b/src/Data/misfits.jl
index 7b7fea6b..4dde5c63 100644
--- a/src/Data/misfits.jl
+++ b/src/Data/misfits.jl
@@ -54,7 +54,6 @@ function func_grad!(pa::P_misfit, grad=nothing)
 		else
 			error("invalid grad")
 		end
-
 		J += JJ 
 	end
 	return J
diff --git a/src/Gallery/fwi.jl b/src/Gallery/fwi.jl
index 7706c2f6..ba547ce6 100644
--- a/src/Gallery/fwi.jl
+++ b/src/Gallery/fwi.jl
@@ -1,7 +1,7 @@
 
 # create pizza problem
 
-function xfwi_problem(attrib::Symbol; born_flag=false, rfields=[:Vx,:Vz,:P])
+function xfwi_problem(attrib::Symbol; born_flag=false, rfields=[:Vx,:Vz,:P], α=0.0)
 
 	if(attrib==:pizza)
 		# starting model
@@ -27,6 +27,22 @@ function xfwi_problem(attrib::Symbol; born_flag=false, rfields=[:Vx,:Vz,:P])
 		igrid=Grid.M2D_resamp(model.mgrid, 50.,50.,)
 		parameterization=[:χvp, :χρ, :null]
 		igrid_interp_scheme=:B2
+	elseif(attrib==:downhole)
+		gri = Grid.M2D(-30.,30.,-120.,30.,1.0,1.0)
+		model = Models.Seismic_zeros(gri)
+		Models.adjust_bounds!(model, [2500., 3500.], [2500., 3500.], [2500.,3500.])
+		model0=deepcopy(model)
+		print(model)
+		Models.Seismic_addon!(model, randn_perc=0.1, fields=[:χvp,:χρ])
+		Models.Seismic_addon!(model0, randn_perc=0.1, fields=[:χvp,:χρ])
+		Models.Seismic_addon!(model, ellip_rad=[2000., 5.], ellip_loc=[0.,0.],ellip_pert=0.1, ellip_α=α, fields=[:χvp],)
+
+		acqgeom=Acquisition.Geom_fixed(-75., -50., 0.0, -100., -50., 0.0, 2, 20, :vertical, :vertical)
+		acqsrc=Acquisition.Src_fixed_mod(acqgeom.nss,1,[:P],mod=model,nλ=4, tmaxfrac=0.8)
+		tgrid=acqsrc.tgrid 
+		igrid=Grid.M2D(-30.,30.,-40.,30.,40,100)
+		parameterization=[:χvp, :χρ, :null]
+		igrid_interp_scheme=:B2
 	else
 		error("invalid attrib")
 	end
diff --git a/test/Fdtd/backprop.jl b/test/Fdtd/backprop.jl
index 8b5eb810..ba16face 100644
--- a/test/Fdtd/backprop.jl
+++ b/test/Fdtd/backprop.jl
@@ -6,7 +6,7 @@ acqgeom =JuMIT.Acquisition.Geom_circ(nss=4,nr=100,rad=[990.,990.]);
 acqgeom =JuMIT.Acquisition.Geom_circ(nss=4,nr=100,rad=[0.,200.]);
 acqsrc=JuMIT.Acquisition.Src_fixed_mod(acqgeom.nss,1,[:P],mod=model, nλ=3, tmaxfrac=0.4)
 
-for sflags in [[1,-1],[2,3]]
+for sflags in [[1,-1],[2,-2]]
 	pa=JuMIT.Fdtd.Param(born_flag=false,npw=1, tgridmod=acqsrc.tgrid,
 	#	abs_trbl=[:null],
 		gmodel_flag=false,