Commit d4ee553f authored by Florian Wittkamp's avatar Florian Wittkamp

Bugfix: Acoustic free surface (merge issue)

parent 744ec46b
...@@ -36,6 +36,7 @@ void surface_acoustic_PML(int ndepth, float ** sp){ ...@@ -36,6 +36,7 @@ void surface_acoustic_PML(int ndepth, float ** sp){
j=ndepth; /* The free surface is located exactly in y=1/2*dh !! */ j=ndepth; /* The free surface is located exactly in y=1/2*dh !! */
for (i=1;i<=NX;i++){ for (i=1;i<=NX;i++){
sp[j][i] = 0;
for (m=1; m<=fdoh; m++) { for (m=1; m<=fdoh; m++) {
sp[j-m][i] = -sp[j+m][i]; sp[j-m][i] = -sp[j+m][i];
} }
......
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