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

[percygalvez at geodynamics.org]

Author: percygalvez
Date: 2011-06-29 00:41:57 -0700 (Wed, 29 Jun 2011)
New Revision: 18673

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/compute_forces_elastic_noDev.f90
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/iterate_time.f90
Log:
loading .f90 files from brutus splay_fault run

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/compute_forces_elastic_noDev.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/compute_forces_elastic_noDev.f90	2011-06-29 07:15:20 UTC (rev 18672)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/compute_forces_elastic_noDev.f90	2011-06-29 07:41:57 UTC (rev 18673)
@@ -234,8 +234,9 @@
        endif
 
       ! Kelvin Voigt damping: artificial viscosity around dynamic faults
+
        if (allocated(Kelvin_Voigt_eta)) then
-         eta = Kelvin_Voigt_eta(ispec)
+         eta = Kelvin_Voigt_eta(ispec)   
          do k=1,NGLLZ
            do j=1,NGLLY
              do i=1,NGLLX

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90	2011-06-29 07:15:20 UTC (rev 18672)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90	2011-06-29 07:41:57 UTC (rev 18673)
@@ -524,8 +524,6 @@
 
 ! Solve for shear stress
   tnorm = sqrt( T(1,:)*T(1,:) + T(2,:)*T(2,:))
-! TO DO :
-! To avoid division by cero (temporal)
   tnorm = max(tnorm,1e0_CUSTOM_REAL)
   t1 = T(1,:)/tnorm
   t2 = T(2,:)/tnorm

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/iterate_time.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/iterate_time.f90	2011-06-29 07:15:20 UTC (rev 18672)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/iterate_time.f90	2011-06-29 07:41:57 UTC (rev 18673)
@@ -64,6 +64,7 @@
 
   do it = 1,NSTEP
   
+    write(IMAIN,*) 'it:',it !test 
     ! simulation status output and stability check
     if(mod(it,NTSTEP_BETWEEN_OUTPUT_INFO) == 0 .or. it == 5) then
       call it_check_stability()    
@@ -79,6 +80,7 @@
     ! elastic solver
     if( ELASTIC_SIMULATION ) call compute_forces_elastic()
 
+    write(IMAIN,*) 'compute_forces_elastic_done' !test 
     ! poroelastic solver
     if( POROELASTIC_SIMULATION ) stop 'poroelastic simulation not implemented yet'