[cig-commits] r22402 - seismo/3D/SPECFEM3D/trunk/src/specfem3D

[danielpeter at geodynamics.org]

Author: danielpeter
Date: 2013-06-24 04:51:50 -0700 (Mon, 24 Jun 2013)
New Revision: 22402

Modified:
   seismo/3D/SPECFEM3D/trunk/src/specfem3D/gravity_perturbation.f90
Log:
updates user output in gravity_init()

Modified: seismo/3D/SPECFEM3D/trunk/src/specfem3D/gravity_perturbation.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/specfem3D/gravity_perturbation.f90	2013-06-24 11:25:18 UTC (rev 22401)
+++ seismo/3D/SPECFEM3D/trunk/src/specfem3D/gravity_perturbation.f90	2013-06-24 11:51:50 UTC (rev 22402)
@@ -31,7 +31,7 @@
        xstore, ystore, zstore, &
        xigll, yigll, zigll, &
        wxgll, wygll, wzgll, &
-       NGNOD, ibool, myrank
+       NGNOD, ibool, myrank, IMAIN
   use specfem_par_elastic, only : rho_vs
   implicit none
 
@@ -48,13 +48,28 @@
 
   open(unit=IIN_G,file='../DATA/gravity_stations',status='old',iostat=ier)
   if( ier /= 0 ) then
-     write(6,*) 'No Gravity stations to compute:  assume not a gravity simulation'
-     return
+    ! user output
+    if( myrank == 0 ) then
+      write(IMAIN,*)
+      write(IMAIN,*) 'no gravity simulation'
+      write(IMAIN,*)
+    endif
+
+    return
   endif
 
   GRAVITY_SIMULATION = .true.
 
   read(IIN_G,*) nstat,ntimgap
+
+  ! user output
+  if( myrank == 0 ) then
+    write(IMAIN,*)
+    write(IMAIN,*) 'incorporating gravity simulation'
+    write(IMAIN,*) '    gravity stations: ',nstat
+    write(IMAIN,*)
+  endif
+
   allocate(xstat(nstat))
   allocate(ystat(nstat))
   allocate(zstat(nstat))