Adjust time-window such that observed and synthetic data can be expressed as:...

Adjust time-window such that observed and synthetic  data can be expressed as:  d_windowed=d exp(-gamma t).
parent c82a0201
...@@ -139,10 +139,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc, ...@@ -139,10 +139,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
time = (float)(j * DT); time = (float)(j * DT);
dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]); dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]);
taper = exp(-GAMMA*dumpa*dumpa); taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]); dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]);
taper1 = exp(-GAMMA*dumpb*dumpb); taper1 = exp(-GAMMA*dumpb);
if(time>=picked_times_m[1][i]+picked_times_m[2][i]){ if(time>=picked_times_m[1][i]+picked_times_m[2][i]){
sectiondata[i][j] = sectiondata[i][j] * taper;} sectiondata[i][j] = sectiondata[i][j] * taper;}
...@@ -160,10 +160,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc, ...@@ -160,10 +160,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
time = (float)(j * DT); time = (float)(j * DT);
dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]); dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]);
taper = exp(-GAMMA*dumpa*dumpa); taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]); dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]);
taper1 = exp(-GAMMA*dumpb*dumpb); taper1 = exp(-GAMMA*dumpb);
dummysection[i][j] = sectiondata[i][j]; dummysection[i][j] = sectiondata[i][j];
...@@ -174,10 +174,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc, ...@@ -174,10 +174,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
dummysection[i][j] = sectiondata[i][j] * taper1;} dummysection[i][j] = sectiondata[i][j] * taper1;}
dumpc = (time-picked_times_m[4][i]-picked_times_m[5][i]); dumpc = (time-picked_times_m[4][i]-picked_times_m[5][i]);
taper2 = exp(-GAMMA*dumpc*dumpc); taper2 = exp(-GAMMA*dumpc);
dumpd = (time-picked_times_m[4][i]+picked_times_m[6][i]); dumpd = (time-picked_times_m[4][i]+picked_times_m[6][i]);
taper3 = exp(-GAMMA*dumpd*dumpd); taper3 = exp(-GAMMA*dumpd);
if(time>=picked_times_m[4][i]+picked_times_m[5][i]){ if(time>=picked_times_m[4][i]+picked_times_m[5][i]){
sectiondata[i][j] = sectiondata[i][j] * taper2;} sectiondata[i][j] = sectiondata[i][j] * taper2;}
...@@ -196,10 +196,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc, ...@@ -196,10 +196,10 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
time = (float)(j * DT); time = (float)(j * DT);
dumpa = (time-picked_times[i]-TWLENGTH_PLUS); dumpa = (time-picked_times[i]-TWLENGTH_PLUS);
taper = exp(-GAMMA*dumpa*dumpa); taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times[i]+TWLENGTH_MINUS); dumpb = (time-picked_times[i]+TWLENGTH_MINUS);
taper1 = exp(-GAMMA*dumpb*dumpb); taper1 = exp(-GAMMA*dumpb);
if(time>=picked_times[i]+TWLENGTH_PLUS){ if(time>=picked_times[i]+TWLENGTH_PLUS){
sectiondata[i][j] = sectiondata[i][j] * taper;} sectiondata[i][j] = sectiondata[i][j] * taper;}
...@@ -221,4 +221,4 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc, ...@@ -221,4 +221,4 @@ void time_window(float **sectiondata, int iter, int ntr_glob, int **recpos_loc,
free_vector(picked_times,1,ntr_glob); free_vector(picked_times,1,ntr_glob);
} }
} /* end of function time_window.c */ } /* end of function time_window.c */
\ No newline at end of file
...@@ -101,10 +101,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i ...@@ -101,10 +101,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
time = (float)(j * DT); time = (float)(j * DT);
dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]); dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]);
taper = exp(-GAMMA*dumpa*dumpa); taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]); dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]);
taper1 = exp(-GAMMA*dumpb*dumpb); taper1 = exp(-GAMMA*dumpb);
if(time>=picked_times_m[1][i]+picked_times_m[2][i]){ if(time>=picked_times_m[1][i]+picked_times_m[2][i]){
sectiondata[i][j] = sectiondata[i][j] * taper;} sectiondata[i][j] = sectiondata[i][j] * taper;}
...@@ -123,10 +123,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i ...@@ -123,10 +123,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
time = (float)(j * DT); time = (float)(j * DT);
dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]); dumpa = (time-picked_times_m[1][i]-picked_times_m[2][i]);
taper = exp(-GAMMA*dumpa*dumpa); taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]); dumpb = (time-picked_times_m[1][i]+picked_times_m[3][i]);
taper1 = exp(-GAMMA*dumpb*dumpb); taper1 = exp(-GAMMA*dumpb);
dummysection[i][j] = sectiondata[i][j]; dummysection[i][j] = sectiondata[i][j];
...@@ -137,10 +137,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i ...@@ -137,10 +137,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
dummysection[i][j] = sectiondata[i][j] * taper1;} dummysection[i][j] = sectiondata[i][j] * taper1;}
dumpc = (time-picked_times_m[4][i]-picked_times_m[5][i]); dumpc = (time-picked_times_m[4][i]-picked_times_m[5][i]);
taper2 = exp(-GAMMA*dumpc*dumpc); taper2 = exp(-GAMMA*dumpc);
dumpd = (time-picked_times_m[4][i]+picked_times_m[6][i]); dumpd = (time-picked_times_m[4][i]+picked_times_m[6][i]);
taper3 = exp(-GAMMA*dumpd*dumpd); taper3 = exp(-GAMMA*dumpd);
if(time>=picked_times_m[4][i]+picked_times_m[5][i]){ if(time>=picked_times_m[4][i]+picked_times_m[5][i]){
sectiondata[i][j] = sectiondata[i][j] * taper2;} sectiondata[i][j] = sectiondata[i][j] * taper2;}
...@@ -159,10 +159,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i ...@@ -159,10 +159,10 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
time = (float)(j * DT); time = (float)(j * DT);
dumpa = (time-picked_times[i]-TWLENGTH_PLUS); dumpa = (time-picked_times[i]-TWLENGTH_PLUS);
taper = exp(-GAMMA*dumpa*dumpa); taper = exp(-GAMMA*dumpa);
dumpb = (time-picked_times[i]+TWLENGTH_MINUS); dumpb = (time-picked_times[i]+TWLENGTH_MINUS);
taper1 = exp(-GAMMA*dumpb*dumpb); taper1 = exp(-GAMMA*dumpb);
if(time>=picked_times[i]+TWLENGTH_PLUS){ if(time>=picked_times[i]+TWLENGTH_PLUS){
sectiondata[i][j] = sectiondata[i][j] * taper;} sectiondata[i][j] = sectiondata[i][j] * taper;}
...@@ -184,4 +184,4 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i ...@@ -184,4 +184,4 @@ void time_window_glob(float **sectiondata, int iter, int ntr_glob, int ns, int i
free_vector(picked_times,1,ntr_glob); free_vector(picked_times,1,ntr_glob);
} }
} /* end of function time_window_glob.c */ } /* end of function time_window_glob.c */
\ No newline at end of file
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment