[cig-commits] r20934 - seismo/3D/FAULT_SOURCE/branches/new_fault_db/src

[surendra at geodynamics.org]

Author: surendra
Date: 2012-10-25 19:13:20 -0700 (Thu, 25 Oct 2012)
New Revision: 20934

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_common.f90
Log:
More compile-time errors fixed

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_common.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_common.f90	2012-10-26 01:05:03 UTC (rev 20933)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_common.f90	2012-10-26 02:13:20 UTC (rev 20934)
@@ -30,19 +30,19 @@
 
 !---------------------------------------------------------------------
 
-subroutine initialize_fault (bc,IIN_BIN,Minv,dt)
+subroutine initialize_fault (bc,IIN_BIN,Minv,dt_tmp)
 
   use specfem_par
 
   class(fault_type), intent(inout) :: bc
   real(kind=CUSTOM_REAL), intent(in)  :: Minv(:)
   integer, intent(in)                 :: IIN_BIN
-  real(kind=CUSTOM_REAL), intent(in)  :: dt
+  real(kind=CUSTOM_REAL), intent(in)  :: dt_tmp
 
   real(kind=CUSTOM_REAL) :: tmp_vec(3,NGLOB_AB)
   real(kind=CUSTOM_REAL), dimension(:,:), allocatable   :: jacobian2Dw
   real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: normal
-  real(kind=CUSTOM_REAL), dimension(:), allocatable :: nxyz
+  real(kind=CUSTOM_REAL), dimension(:,:), allocatable :: nxyz
   integer, dimension(:,:), allocatable :: ibool1
   integer :: ij,k,e
 
@@ -72,7 +72,7 @@
     read(IIN_BIN) bc%coord(2,:)
     read(IIN_BIN) bc%coord(3,:)
 
-    bc%dt = dt
+    bc%dt = dt_tmp
 
     bc%B = 0e0_CUSTOM_REAL
     allocate(nxyz(3,bc%nglob))
@@ -120,7 +120,7 @@
     !   Z = 1/( B1/M1 + B2/M2 ) / (0.5*dt)
     ! T_stick = Z*Vfree traction as if the fault was stuck (no displ discontinuity) 
     ! NOTE: same Bi on both sides, see note above
-    bc%Z = 1.e0_CUSTOM_REAL/(0.5e0_CUSTOM_REAL*dt * bc%B *( bc%invM1 + bc%invM2 ))
+    bc%Z = 1.e0_CUSTOM_REAL/(0.5e0_CUSTOM_REAL*dt_tmp * bc%B *( bc%invM1 + bc%invM2 ))
     ! WARNING: In non-split nodes at fault edges M is assembled across the fault.
     ! hence invM1+invM2=2/(M1+M2) instead of 1/M1+1/M2
     ! In a symmetric mesh (M1=M2) Z will be twice its intended value