Skip to content

Commit

Permalink
Add files via upload
Browse files Browse the repository at this point in the history 
I changed the directions 2 and 3 each other.
  • Loading branch information
@ohjuwon
ohjuwon authored Sep 8, 2016
1 parent 6c09c75 commit 115bd1b
Showing 1 changed file with 10 additions and 10 deletions.
20 changes: 10 additions & 10 deletions SOFI3D_org/src/hh_elastic.c
View file
Original file line number Diff line number Diff line change
Expand Up @@ -54,8 +54,8 @@ void model_elastic(float *** rho, float *** pi, float *** u,
/* x=1, y=2 in Tsvankin [1997] (e.g.) epsx=epsion1 & epsy=epsilon2 */

/* parameters for layer 1 */
const float vpv1=3500.0, poi1=0.25, epsx1=0.1, epsy1=0.2, delx1=0.1, dely1=0.3, delxy1=0.1,
gamx1=0.3, gamy1=0.1, rho1=2000.0, h=100000.0;
const float vpv1=3500.0, poi1=0.25, epsx1=0., epsy1=0.2, delx1=0., dely1=0., delxy1=0.,
gamx1=0., gamy1=0., rho1=2000.0, h=100000.0;

/* parameters for layer 2 */
//const float vp2=5700.0, vs2=3400.0, rho2=2500.0;
Expand Down Expand Up @@ -99,7 +99,7 @@ void model_elastic(float *** rho, float *** pi, float *** u,

else{
Vpv=vpv2; Poi=poi2; Epsx=epsx2; Epsy=epsy2; Delx=delx2; Dely=dely2;
Delxy=delxy2; Gamx=gamx2; Gamy=gamy2; Rho=rho2; }
Delxy=delxy1; Gamx=gamx2; Gamy=gamy2; Rho=rho2; }

/*=========================================================
* modify up to this point for ELASTIC model definition
Expand All @@ -122,17 +122,17 @@ void model_elastic(float *** rho, float *** pi, float *** u,
u[jj][ii][kk]=muv;
pi[jj][ii][kk]=piv;

C33[jj][ii][kk]=Rho*Vpv*Vpv;
C22[jj][ii][kk]=Rho*Vpv*Vpv;
C55[jj][ii][kk]=Rho*Vsv*Vsv;
C11[jj][ii][kk]=(1+2*Epsx)*Rho*Vpv*Vpv;
C22[jj][ii][kk]=(1+2*Epsy)*Rho*Vpv*Vpv;
C66[jj][ii][kk]=(1+2*Gamx)*Rho*Vsv*Vsv;
C44[jj][ii][kk]=((1+2*Gamx)*Rho*Vsv*Vsv)/(1+2*Gamy);
C33[jj][ii][kk]=(1+2*Epsy)*Rho*Vpv*Vpv;
C44[jj][ii][kk]=(1+2*Gamx)*Rho*Vsv*Vsv;
C66[jj][ii][kk]=((1+2*Gamx)*Rho*Vsv*Vsv)/(1+2*Gamy);
C13[jj][ii][kk]=Rho*sqrt((Vpv*Vpv-Vsv*Vsv)*((1+2*Delx)*Vpv*Vpv-Vsv*Vsv))-Rho*Vsv*Vsv ;
C23[jj][ii][kk]=Rho*sqrt((Vpv*Vpv-(1+2*Gamx)*Vsv*Vsv/(1+2*Gamy))*((1+2*Dely)*Vpv*Vpv-(1+2*Gamx)*Vsv*Vsv/(1+2*Gamy)))-Rho*(1+2*Gamx)*Vsv*Vsv/(1+2*Gamy);
C12[jj][ii][kk]=Rho*sqrt(((1+2*Epsx)*Vpv*Vpv-(1+2*Gamx)*Vsv*Vsv)*((1+2*Delxy)*(1+2*Epsx)*Vpv*Vpv-(1+2*Gamx)*Vsv*Vsv))-Rho*(1+2*Gamx)*Vsv*Vsv ;
C12[jj][ii][kk]=Rho*sqrt((Vpv*Vpv-(1+2*Gamx)*Vsv*Vsv/(1+2*Gamy))*((1+2*Dely)*Vpv*Vpv-(1+2*Gamx)*Vsv*Vsv/(1+2*Gamy)))-Rho*(1+2*Gamx)*Vsv*Vsv/(1+2*Gamy);
C23[jj][ii][kk]=Rho*sqrt(((1+2*Epsx)*Vpv*Vpv-(1+2*Gamx)*Vsv*Vsv)*((1+2*Delxy)*(1+2*Epsx)*Vpv*Vpv-(1+2*Gamx)*Vsv*Vsv))-Rho*(1+2*Gamx)*Vsv*Vsv ;
rho[jj][ii][kk]=Rho;


if (WRITE_MODELFILES==1) {
vpv[jj][ii][kk]=Vpv;
Expand Down

0 comments on commit 115bd1b

Please sign in to comment.