Commit a3da04c0 authored by niklas.thiel's avatar niklas.thiel

BUGFIX calculating receiver positions

parent 2e14b8db
...@@ -114,7 +114,6 @@ int **receiver(int *ntr, float **srcpos, int shotno){ ...@@ -114,7 +114,6 @@ int **receiver(int *ntr, float **srcpos, int shotno){
nxrec2=iround(XREC1/DH); /* in gridpoints */ nxrec2=iround(XREC1/DH); /* in gridpoints */
nyrec2=iround(YREC1/DH); nyrec2=iround(YREC1/DH);
offset=iround((XREC2-srcpos[1][1])/DH); offset=iround((XREC2-srcpos[1][1])/DH);
} }
if (nyrec1 != nyrec2){ if (nyrec1 != nyrec2){
fprintf(FP,"\n\n"); fprintf(FP,"\n\n");
...@@ -126,7 +125,7 @@ int **receiver(int *ntr, float **srcpos, int shotno){ ...@@ -126,7 +125,7 @@ int **receiver(int *ntr, float **srcpos, int shotno){
if (offset<0) recdist=-NGEOPH; if (offset<0) recdist=-NGEOPH;
else recdist=NGEOPH; else recdist=NGEOPH;
*ntr=iround((nxrec2-nxrec1)/recdist)+1; *ntr=iround((nxrec2-nxrec1)/recdist)+1;
if (shotno>0) nxrec1=nxrec1+(shotdist*(shotno-1)); if (shotno>0) nxrec1=nxrec1+((srcpos[1][shotno]-srcpos[1][1])/DH);
recpos=imatrix(1,3,1,*ntr); recpos=imatrix(1,3,1,*ntr);
n=0; n=0;
Markdown is supported
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment