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

[ampuero at geodynamics.org]

Author: ampuero
Date: 2012-10-26 23:23:01 -0700 (Fri, 26 Oct 2012)
New Revision: 20956

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_generate_databases.f90
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90
Log:
finished removing FAULT directory; re-fixed rsf init

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_generate_databases.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_generate_databases.f90	2012-10-27 06:08:13 UTC (rev 20955)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_generate_databases.f90	2012-10-27 06:23:01 UTC (rev 20956)
@@ -91,7 +91,7 @@
  
  ! read fault input file
   nb = 0
-  open(unit=IIN,file='DATA/FAULT/Par_file_faults',status='old',action='read',iostat=ier)
+  open(unit=IIN,file='DATA/Par_file_faults',status='old',action='read',iostat=ier)
   if (ier==0) then    
     read(IIN,*) nb  
   else

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90	2012-10-27 06:08:13 UTC (rev 20955)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90	2012-10-27 06:23:01 UTC (rev 20956)
@@ -856,7 +856,7 @@
                         - f%f0 - f%a*log(f%V_init/f%V0) ) &
                       / f%b )
   else
-     f%theta =  f%a * log(TWO*f%V0/V_init * sinh(-sqrt(T0(1,:)**2+T0(2,:)**2)/T0(3,:)/f%a))
+     f%theta =  f%a * log(TWO*f%V0/f%V_init * sinh(-sqrt(T0(1,:)**2+T0(2,:)**2)/T0(3,:)/f%a))
   endif
 
  ! WARNING : ad hoc for SCEC benchmark TPV10x
@@ -868,7 +868,7 @@
   call init_2d_distribution(nucFload,coord,IIN_PAR,nFload)
 
  ! WARNING: the line below is only valid for pure strike-slip faulting
-  V(1,:) = bc%rsf%V_init
+  V(1,:) = f%V_init
 
 end subroutine rsf_init