[cig-commits] r17879 - seismo/2D/SPECFEM2D/trunk

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2011-02-15 16:36:16 -0800 (Tue, 15 Feb 2011)
New Revision: 17879

Modified:
   seismo/2D/SPECFEM2D/trunk/part_unstruct.F90
Log:
converted comments to print statements in subroutine read_abs_surface()


Modified: seismo/2D/SPECFEM2D/trunk/part_unstruct.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/part_unstruct.F90	2011-02-16 00:29:33 UTC (rev 17878)
+++ seismo/2D/SPECFEM2D/trunk/part_unstruct.F90	2011-02-16 00:36:16 UTC (rev 17879)
@@ -291,7 +291,7 @@
 
   open(unit=994, file=trim(filename), form='formatted' , status='old', action='read', iostat=ier)
   if( ier /= 0 ) then
-    print*,'error opening file: ',trim(filename)
+    print *,'error opening file: ',trim(filename)
     stop 'error read absorbing surface file'
   endif
 
@@ -300,15 +300,17 @@
   allocate(abs_surface(4,nelemabs))
 
   do i = 1, nelemabs
-     read(994,*) abs_surface(1,i), abs_surface(2,i), abs_surface(3,i), abs_surface(4,i)
-! The input format is currently limited: only two nodes per element can be listed.
-! If one of your elements has more than one edge along a given absorbing contour (e.g., if that contour has a corner)
-! then list it twice, putting the first edge on the first line and the second edge on the second line.
-! if one of your elements has a single point along the absording contour rather than a full edge, do NOT list it
-! (it would have no weight in the contour integral anyway because it would consist of a single point).
-! If you are using 9-node elements, list only the first and last point of the edge and not the intermediate point
-! located around the middle of the edge; the right 9-node curvature will be restored automatically by the code.
-     if (abs_surface(2,i) /= 2) stop 'only two nodes per element should be listed for absorbing edges'
+    read(994,*) abs_surface(1,i), abs_surface(2,i), abs_surface(3,i), abs_surface(4,i)
+    if (abs_surface(2,i) /= 2) then
+      print *,'The input format is currently limited: only two nodes per element can be listed.'
+      print *,'If one of your elements has more than one edge along a given absorbing contour (e.g., if that contour has a corner)'
+      print *,'then list it twice, putting the first edge on the first line and the second edge on the second line.'
+      print *,'if one of your elements has a single point along the absording contour rather than a full edge, do NOT list it'
+      print *,'(it would have no weight in the contour integral anyway because it would consist of a single point).'
+      print *,'If you are using 9-node elements, list only the first and last point of the edge and not the intermediate point'
+      print *,'located around the middle of the edge; the right 9-node curvature will be restored automatically by the code.'
+      stop 'only two nodes per element should be listed for absorbing edges'
+    endif
   enddo
 
   close(994)