Fix error for velocity updates in elastic SH-waves

parent 6d2fbf22
...@@ -103,8 +103,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f ...@@ -103,8 +103,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx); sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy); syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT; uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT; uyz[j][i]=(f*vzy)/DT;
} }
} }
...@@ -159,8 +159,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f ...@@ -159,8 +159,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx); sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy); syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT; uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT; uyz[j][i]=(f*vzy)/DT;
} }
} }
...@@ -216,8 +216,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f ...@@ -216,8 +216,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx); sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy); syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT; uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT; uyz[j][i]=(f*vzy)/DT;
} }
} }
...@@ -275,8 +275,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f ...@@ -275,8 +275,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx); sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy); syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT; uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT; uyz[j][i]=(f*vzy)/DT;
} }
} }
...@@ -336,8 +336,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f ...@@ -336,8 +336,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx); sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy); syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT; uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT; uyz[j][i]=(f*vzy)/DT;
} }
} }
...@@ -399,8 +399,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f ...@@ -399,8 +399,8 @@ void update_s_elastic_PML_SH(int nx1, int nx2, int ny1, int ny2, float ** vz, f
sxz[j][i]+=(fipjp*vzx); sxz[j][i]+=(fipjp*vzx);
syz[j][i]+=(f*vzy); syz[j][i]+=(f*vzy);
uxz[j][i]+=(fipjp*vzx)/DT; uxz[j][i]=(fipjp*vzx)/DT;
uyz[j][i]+=(f*vzy)/DT; uyz[j][i]=(f*vzy)/DT;
} }
} }
......
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