diff --git a/marchenko/Makefile b/marchenko/Makefile index 03f9a5b0428649ab0e7eac4e0dab74fa24aead7a..f60dc90c6a7bf57855e730656e13683f8dfe908a 100644 --- a/marchenko/Makefile +++ b/marchenko/Makefile @@ -61,7 +61,7 @@ fmute: $(OBJJ) OBJH = $(SRCH:%.c=%.o) marchenko: $(OBJH) - $(CC) $(LDFLAGS) $(OPTC) $(CFLAGS) -o marchenko1 $(OBJH) $(LIBS) + $(CC) $(LDFLAGS) $(OPTC) $(CFLAGS) -o marchenko $(OBJH) $(LIBS) OBJP = $(SRCP:%.c=%.o) diff --git a/marchenko/marchenko.c b/marchenko/marchenko.c index ca54d5a1c47c97c966bfa923008420d34461eb16..18f36d295f6b3a6962b029f9969e4b39b49521c6 100644 --- a/marchenko/marchenko.c +++ b/marchenko/marchenko.c @@ -581,24 +581,9 @@ int main (int argc, char **argv) } /* Apply mute with window for Gmin */ if ( plane_wave==1 ) { - /* 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 < 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); + applyMute_tshift(Gmin, muteW, smooth, 4, Nfoc, nxs, nts, ixpos, npos, shift, 0, tsynW); + /* for plane wave with angle tshift downward */ + timeShift(Gmin, nts, npos, dt, tshift, fmin, fmax); } else { applyMute(Gmin, muteW, smooth, 4, Nfoc, nxs, nts, ixpos, npos, shift, tsynW);