[cig-commits] r20249 - seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D

[elliott.sales.de.andrade at geodynamics.org]

Author: elliott.sales.de.andrade
Date: 2012-05-31 12:08:09 -0700 (Thu, 31 May 2012)
New Revision: 20249

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/write_movie_volume.f90
Log:
Fix output of blank lines.

Apparently, empty strings are not allowed in Fortran, and this outputs
an error message on xlf. Because the default flags say to halt on any
message above a warning, this stops the compile. In a small test, not
specifying any variables is sufficient to output a blank line on ifort,
gfortran, and xlf.



Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90	2012-05-31 19:06:57 UTC (rev 20248)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/setup_sources_receivers.f90	2012-05-31 19:08:09 UTC (rev 20249)
@@ -339,7 +339,7 @@
   if(myrank == 0) then
 
     ! finishes vtk file
-    write(IOVTK,*) ""
+    write(IOVTK,*) ! blank line
     close(IOVTK)
 
     !  we should know NSOURCES+nrec at this point...

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/write_movie_volume.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/write_movie_volume.f90	2012-05-31 19:06:57 UTC (rev 20248)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/write_movie_volume.f90	2012-05-31 19:08:09 UTC (rev 20249)
@@ -925,7 +925,7 @@
     !write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
     write(IOVTK,'(3e18.6)') xval,yval,zval
   enddo
-  write(IOVTK,*) ""
+  write(IOVTK,*) ! blank line
 
   ! defines cell on coarse corner points
   ! note: indices for vtk start at 0
@@ -958,12 +958,12 @@
     endif
 
   enddo
-  write(IOVTK,*) ""
+  write(IOVTK,*) ! blank line
 
   ! type: hexahedrons
   write(IOVTK,'(a,i12)') "CELL_TYPES ",nspec
   write(IOVTK,*) (12,ispec=1,nspec)
-  write(IOVTK,*) ""
+  write(IOVTK,*) ! blank line
 
   ! x components
   write(IOVTK,'(a,i12)') "POINT_DATA ",nglob
@@ -992,7 +992,7 @@
                         + glob_data(2,i)*glob_data(2,i) &
                         + glob_data(3,i)*glob_data(3,i))
   enddo
-  write(IOVTK,*) ""
+  write(IOVTK,*) ! blank line
 
   close(IOVTK)
 
@@ -1090,7 +1090,7 @@
         write(IOVTK,'(3e18.6)') xval,yval,zval
       enddo
     enddo
-    write(IOVTK,*) ""
+    write(IOVTK,*) ! blank line
 
     ! defines cell on coarse corner points
     ! note: indices for vtk start at 0
@@ -1129,12 +1129,12 @@
 
       enddo
     enddo
-    write(IOVTK,*) ""
+    write(IOVTK,*) ! blank line
 
     ! type: hexahedrons
     write(IOVTK,'(a,i12)') "CELL_TYPES ",nspec*NPROCTOT_VAL
     write(IOVTK,*) (12,ispec=1,nspec*NPROCTOT_VAL)
-    write(IOVTK,*) ""
+    write(IOVTK,*) ! blank line
 
     ! x components
     write(IOVTK,'(a,i12)') "POINT_DATA ",nglob*NPROCTOT_VAL
@@ -1171,7 +1171,7 @@
                           + store_val_uz_all(i,iproc)**2 )
       enddo
     enddo
-    write(IOVTK,*) ""
+    write(IOVTK,*) ! blank line
 
     close(IOVTK)