Commit 5944d0e4 authored by Florian Wittkamp's avatar Florian Wittkamp

Bugfix: Free surface viscoelastic

parent d4ee553f
...@@ -36,18 +36,17 @@ void surface_PML(int ndepth, float ** vx, float ** vy, float ** sxx, float ** sy ...@@ -36,18 +36,17 @@ void surface_PML(int ndepth, float ** vx, float ** vy, float ** sxx, float ** sy
float vxx, vyy, sump=0.0; float vxx, vyy, sump=0.0;
float dh24, dthalbe; float dh24, dthalbe;
float *pts, ws, sumu, sumpi, mu, pi; float *pts, ws, sumu, sumpi, mu, pi;
extern float DT, DH, *FL, WAVETYPE; extern float DT, DH, *FL;
extern int NX, INVMAT1, L; extern int NX, INVMAT1, L;
extern int FW, BOUNDARY; extern int FW, BOUNDARY;
extern int NPROCX, NPROCY, POS[3], MYID; extern int NPROCX, NPROCY, POS[3], MYID;
extern int FDORDER; extern int FDORDER,WAVETYPE;
extern float F_REF; extern float F_REF;
fdoh = FDORDER/2; fdoh = FDORDER/2;
dthalbe=DT/2.0; dthalbe=DT/2.0;
dh24=1.0/DH; dh24=1.0/DH;
if (WAVETYPE==1||WAVETYPE==3){ if (WAVETYPE==1||WAVETYPE==3){
/* vector for maxwellbodies */ /* vector for maxwellbodies */
pts=vector(1,L); pts=vector(1,L);
......
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