[cig-commits] r18557 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: decompose_mesh_SCOTCH src

[percygalvez at geodynamics.org]

Author: percygalvez
Date: 2011-06-08 15:41:42 -0700 (Wed, 08 Jun 2011)
New Revision: 18557

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90
Log:
fixing kelvin_voigt.bin files

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90	2011-06-08 22:38:08 UTC (rev 18556)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90	2011-06-08 22:41:42 UTC (rev 18557)
@@ -330,9 +330,9 @@
     endif
   enddo
 
- ! list of elements per node
- !   nnodes_elmnts(i) = number of elements containing node #i (i=0:nnodes-1)     
- !   nodes_elmnts(nsize*i:nsize*i+nnodes_elmnts(i)-1) = index of elements (starting at 0) containing node #i  
+! list of elements per node
+!   nnodes_elmnts(i) = number of elements containing node #i (i=0:nnodes-1)     
+!   nodes_elmnts(nsize*i:nsize*i+nnodes_elmnts(i)-1) = index of elements (starting at 0) containing node #i  
 !  nnodes_elmnts(:) = 0
 !  nodes_elmnts(:) = 0
 !  do i = 0, esize*nelmnts-1

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90	2011-06-08 22:38:08 UTC (rev 18556)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90	2011-06-08 22:41:42 UTC (rev 18557)
@@ -177,7 +177,8 @@
   if (.not. ANY_FAULT_IN_THIS_PROC ) return
 
   do iflt=1,size(fault_db)
-
+! TO DO : Percy, iface is computed with fault open , we have to change this and close  *_dummy
+!                at this point.
     call setup_iface(fault_db(iflt),nnodes_ext_mesh,nodes_coords_ext_mesh,nspec,nglob,ibool)
 
     ! saving gll indices for each fault face, needed for ibulks
@@ -610,21 +611,27 @@
 
   integer :: nbfaults,iflt,ier
   character(len=256) :: filename
+  integer :: size_Kelvin_Voigt
 
   if (.not.ANY_FAULT) return
-
+! opening Kelvin_voig_eta.bin (Necessary for all processors , if number of fault elements = 0 
+! then the file will be empty with 
+    filename = prname(1:len_trim(prname))//'Kelvin_voigt_eta.bin'
+  open(unit=IOUT,file=trim(filename),status='unknown',action='write',form='unformatted',iostat=ier)
+  if( ier /= 0 ) then
+    write(IOUT,*) 'error opening file ',trim(filename)
+    stop 
+  endif
+   
 ! saves mesh file proc***_Kelvin_voigt_eta.bin
   if (allocated(Kelvin_Voigt_eta)) then
-    filename = prname(1:len_trim(prname))//'Kelvin_voigt_eta.bin'
-    open(unit=IOUT,file=trim(filename),status='unknown',action='write',form='unformatted',iostat=ier)
-    if( ier /= 0 ) then
-      write(IOUT,*) 'error opening file ',trim(filename)
-      stop 
-    endif
     write(IOUT) size(Kelvin_Voigt_eta)
-    write(IOUT) Kelvin_Voigt_eta  
-    close(IOUT)
+  else 
+    size_Kelvin_Voigt = 0
   endif
+  write(IOUT) size_Kelvin_Voigt
+  if (size_Kelvin_Voigt /= 0) Write(IOUT) Kelvin_Voigt_eta  
+  Close(IOUT)
 
 ! saves mesh file proc***_fault_db.bin
   filename = prname(1:len_trim(prname))//'fault_db.bin'