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

[surendra at geodynamics.org]

Author: surendra
Date: 2012-10-23 22:15:45 -0700 (Tue, 23 Oct 2012)
New Revision: 20895

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90
Log:
Fixed errors related to zero-faults

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	2012-10-24 04:08:16 UTC (rev 20894)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90	2012-10-24 05:15:45 UTC (rev 20895)
@@ -169,6 +169,9 @@
     return 
   endif
 
+  read(IIN_PAR,*) nbfaults
+  if (nbfaults==0) return
+
   filename = prname(1:len_trim(prname))//'fault_db.bin'
   open(unit=IIN_BIN,file=trim(filename),status='old',action='read',form='unformatted',iostat=ier)
   if( ier /= 0 ) then
@@ -177,8 +180,6 @@
   endif
   ! WARNING TO DO: should be an MPI abort
 
-  read(IIN_PAR,*) nbfaults
-  if (nbfaults==0) return
   ! Reading etas of each fault
   do iflt = 1,nbfaults
     read(IIN_PAR,*) ! etas

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90	2012-10-24 04:08:16 UTC (rev 20894)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_kinematic.f90	2012-10-24 05:15:45 UTC (rev 20895)
@@ -113,13 +113,14 @@
     return 
   endif
 
+  read(IIN_PAR,*) nbfaults
+  if (nbfaults==0) return
+
   filename = prname(1:len_trim(prname))//'fault_db.bin'
   open(unit=IIN_BIN,file=trim(filename),status='old',action='read',form='unformatted',iostat=ier)
   if( ier /= 0 ) stop 'Have not found proc*_fault_db.bin'
   ! WARNING TO DO: should be an MPI abort
 
-  read(IIN_PAR,*) nbfaults
-  if (nbfaults==0) return
   read(IIN_PAR,*)  ! eta
   read(IIN_PAR,*) SIMULATION_TYPE
   if ( SIMULATION_TYPE == 2 ) then