diff --git a/extrap/main/migr.c b/extrap/main/migr.c
index 0cebce0d396d0db8808bf9505f90135e434f7997..a1b86384b4c9d0914ab83804d2eae494539ec8fd 100644
--- a/extrap/main/migr.c
+++ b/extrap/main/migr.c
@@ -376,7 +376,7 @@ int main(int argc, char *argv[])
 			vwarn("number of depth steps set to maximum");
 			ndepth = nzm;
 		}
-		if(nx > nxm) verr("number of x position in model to small");
+		if(nx > nxm) verr("number of x position in model(%d) to small (data has %d)", nxm, nx);
 
 		velmod1 = (float *)malloc(nxm*ndepth*sizeof(float));
 		velmod2 = (float *)malloc(nxm*ndepth*sizeof(float));
diff --git a/fdelmodc/decomposition.c b/fdelmodc/decomposition.c
index a47b41e1cc7292522fbc553d9fdc2a461a3d4260..9fb26cd3ce9c55353c88cbaf0c2e7c1950880997 100644
--- a/fdelmodc/decomposition.c
+++ b/fdelmodc/decomposition.c
@@ -429,6 +429,7 @@ void decudF(float om, float rho, float cp, float dx, int nkx, float kangle, floa
 	return;
 }
 
+/*
 void findMatch(complex *A, complex *B, int nx, float om, float trange, float arange, float *topt, float *aopt) 
 {
 	int ix, it, ia;
@@ -463,6 +464,7 @@ void findMatch(complex *A, complex *B, int nx, float om, float trange, float ara
 	free(shifted);
 	return;
 }
+*/
 
 complex froot(float x)
 {
diff --git a/marchenko/Makefile b/marchenko/Makefile
index f60dc90c6a7bf57855e730656e13683f8dfe908a..03f9a5b0428649ab0e7eac4e0dab74fa24aead7a 100644
--- a/marchenko/Makefile
+++ b/marchenko/Makefile
@@ -61,7 +61,7 @@ fmute:	$(OBJJ)
 OBJH	= $(SRCH:%.c=%.o)
 
 marchenko:	$(OBJH) 
-	$(CC) $(LDFLAGS) $(OPTC) $(CFLAGS) -o marchenko $(OBJH) $(LIBS)
+	$(CC) $(LDFLAGS) $(OPTC) $(CFLAGS) -o marchenko1 $(OBJH) $(LIBS)
 
 OBJP	= $(SRCP:%.c=%.o)
 
diff --git a/marchenko/marchenko.c b/marchenko/marchenko.c
index 3fd1f9bca86f7c5267bead3e0da3d75f8997b3cf..ca54d5a1c47c97c966bfa923008420d34461eb16 100644
--- a/marchenko/marchenko.c
+++ b/marchenko/marchenko.c
@@ -581,21 +581,24 @@ int main (int argc, char **argv)
         }
         /* Apply mute with window for Gmin */
 		if ( plane_wave==1 ) {
-            applyMute_tshift(Gmin, muteW, smooth, 4, Nfoc, nxs, nts, ixpos, npos, shift, 0, tsynW);
-			itmin = 0;
+            /* for plane wave with angle shift itmin downward */
+			int *tsynSW;
+    		tsynSW = (int *)malloc(Nfoc*nxs*sizeof(int));
+			itmin = NINT(tshift/dt);
+            for (i = 0; i < nxs; i++) tsynSW[i] = tsynW[i] + itmin;
+            //timeShift(Gmin, nts, npos, dt, tshift, fmin, fmax);
+
+            /* applyMute_tshift(Gmin, muteW, smooth, 4, Nfoc, nxs, nts, ixpos, npos, shift, 0, tsynSW);
             for (l = 0; l < Nfoc; l++) {
                 for (i = 0; i < npos; i++) {
                     memcpy(&trace[0],&Gmin[l*nxs*nts+i*nts],nts*sizeof(float));
-                    for (j = 0; j < itmin; j++) {
-                        Gmin[l*nxs*nts+i*nts+j] = 0.0;
-                    }
                     for (j = 0; j < nts; j++) {
                         Gmin[l*nxs*nts+i*nts+j] = trace[j];
                     }
                 }
-            }
+            } */
             /* for plane wave with angle shift itmin downward */
-            timeShift(Gmin, nts, npos, dt, tshift, fmin, fmax);
+            //timeShift(Gmin, nts, npos, dt, tshift, fmin, fmax);
 		}
 		else {
             applyMute(Gmin, muteW, smooth, 4, Nfoc, nxs, nts, ixpos, npos, shift, tsynW);