[cig-commits] r21313 - in seismo/2D/SPECFEM2D/trunk/src: meshfem2D specfem2D

[cristini at geodynamics.org]

Author: cristini
Date: 2013-02-01 01:42:24 -0800 (Fri, 01 Feb 2013)
New Revision: 21313

Modified:
   seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_materials.f90
   seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90
Log:
change reading of Par_file for anisotropic media

Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_materials.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_materials.f90	2013-02-01 04:55:00 UTC (rev 21312)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_materials.f90	2013-02-01 09:42:24 UTC (rev 21313)
@@ -131,16 +131,16 @@
         ! anisotropic materials
 
         rho_s(i) = val0read
-        cp(i) = val1read
-        cs(i) = val2read
-        aniso3(i) = val3read
-        aniso4(i) = val4read
-        aniso5(i) = val5read
-        aniso6(i) = val6read
-        aniso7(i) = val7read
-        aniso8(i) = val8read
-        QKappa(i) = val9read
-        Qmu(i) = val10read
+        aniso3(i) = val1read
+        aniso4(i) = val2read
+        aniso5(i) = val3read
+        aniso6(i) = val4read
+        aniso7(i) = val5read
+        aniso8(i) = val6read
+        QKappa(i) = val7read
+        Qmu(i) = val8read
+        cp(i) = sqrt(val4read/val0read)
+        cs(i) = sqrt(val6read/val0read)
      else
 
         ! poroelastic materials

Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90	2013-02-01 04:55:00 UTC (rev 21312)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90	2013-02-01 09:42:24 UTC (rev 21313)
@@ -278,9 +278,9 @@
                     permxx(i),permxz(i),permzz(i),kappa_s(i),&
                     kappa_f(i),kappa_fr(i),eta_f(i),mu_fr(i),Qmu(i)
       else
-         write(15,*) i,icodemat(i),rho_s(i),cp(i),cs(i), &
+         write(15,*) i,icodemat(i),rho_s(i), &
                     aniso3(i),aniso4(i),aniso5(i),aniso6(i),&
-                    aniso7(i),aniso8(i),QKappa(i),Qmu(i),0,0
+                    aniso7(i),aniso8(i),QKappa(i),Qmu(i),0,0,0,0
       endif
     enddo
 

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90	2013-02-01 04:55:00 UTC (rev 21312)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/gmat01.f90	2013-02-01 09:42:24 UTC (rev 21313)
@@ -133,10 +133,6 @@
 
         density(1) =val0
 
-        ! P and S velocity
-        cp = 20
-        cs = 10
-
         ! Anisotropy parameters
         c11 = val1
         c13 = val2
@@ -145,6 +141,10 @@
         c35 = val5
         c55 = val6
 
+        ! P and S velocity
+        cp = sqrt(c33/density(1))
+        cs = sqrt(c55/density(1))
+
         ! QKappa and Qmu values
         !QKappa = val9
         !Qmu = val10