[cig-commits] [commit] QA: added the station name to output_list_stations.txt. also removed an unused variable from src/specfem3D/gravity_perturbation.f90 (ce90249)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/specfem3d

On branch  : QA
Link       : https://github.com/geodynamics/specfem3d/compare/232818c4f9cbfeabaa93541c3898060edc5516a0...1fe92c91d0ab35e48fb6846c196d60390e9df02f

>---------------------------------------------------------------

commit ce90249acd62cb6a5902ab24836bb8e7e378bf14
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Mon Jan 6 23:56:47 2014 +0100

    added the station name to output_list_stations.txt.
    also removed an unused variable from src/specfem3D/gravity_perturbation.f90


>---------------------------------------------------------------

ce90249acd62cb6a5902ab24836bb8e7e378bf14
 src/specfem3D/gravity_perturbation.f90 | 2 +-
 src/specfem3D/locate_receivers.f90     | 4 ++--
 src/specfem3D/write_output_SU.f90      | 8 ++++----
 3 files changed, 7 insertions(+), 7 deletions(-)

diff --git a/src/specfem3D/gravity_perturbation.f90 b/src/specfem3D/gravity_perturbation.f90
index d1c113d..0af3c5f 100644
--- a/src/specfem3D/gravity_perturbation.f90
+++ b/src/specfem3D/gravity_perturbation.f90
@@ -37,7 +37,7 @@ subroutine gravity_init()
 
   integer, parameter :: IIN_G = 367
   integer ier
-  real(kind=CUSTOM_REAL), dimension(NGLLX,NGLLY,NGLLZ):: rho_elem,vs_elem
+  real(kind=CUSTOM_REAL), dimension(NGLLX,NGLLY,NGLLZ):: rho_elem
   integer :: i,j,k,iglob,ispec,istat
   double precision :: Jac3D
   ! coordinates of the control points
diff --git a/src/specfem3D/locate_receivers.f90 b/src/specfem3D/locate_receivers.f90
index 7785f0b..a7d59c2 100644
--- a/src/specfem3D/locate_receivers.f90
+++ b/src/specfem3D/locate_receivers.f90
@@ -226,7 +226,7 @@
           call exit_mpi(myrank,'error opening file '//trim(OUTPUT_FILES)//'/output_list_stations.txt')
 
         do irec=1,nrec
-          write(IOUT_SU,*) x_found(irec),y_found(irec),z_found(irec)
+          write(IOUT_SU,*) station_name(irec),network_name(irec),x_found(irec),y_found(irec),z_found(irec)
         enddo
 
         close(IOUT_SU)
@@ -922,7 +922,7 @@
       call exit_mpi(myrank,'error opening file '//trim(OUTPUT_FILES)//'/output_list_stations.txt')
 
     do irec=1,nrec
-      write(IOUT_SU,*) x_found(irec),y_found(irec),z_found(irec)
+      write(IOUT_SU,*) station_name(irec),network_name(irec),x_found(irec),y_found(irec),z_found(irec)
     enddo
 
     close(IOUT_SU)
diff --git a/src/specfem3D/write_output_SU.f90 b/src/specfem3D/write_output_SU.f90
index 74960c0..9c17a39 100644
--- a/src/specfem3D/write_output_SU.f90
+++ b/src/specfem3D/write_output_SU.f90
@@ -46,19 +46,19 @@
   integer :: i,ier
 
   allocate(x_found(nrec),y_found(nrec),z_found(nrec),stat=ier)
-  if( ier /= 0 ) stop 'error allocating array x_found y_found z_found'
+  if( ier /= 0 ) stop 'error allocating arrays x_found y_found z_found'
 
   ! reads in station locations from output_list file
-  open(unit=IIN_SU1,file=trim(OUTPUT_FILES)//'/output_list_stations.txt',status='unknown',iostat=ier)
+  open(unit=IIN_SU1,file=trim(OUTPUT_FILES)//'/output_list_stations.txt',status='old',iostat=ier)
   if( ier /= 0 ) stop 'error opening output_list_stations.txt file'
 
   do irec=1,nrec
-   read(IIN_SU1,*) x_found(irec),y_found(irec),z_found(irec)
+   read(IIN_SU1,*) station_name(irec),network_name(irec),x_found(irec),y_found(irec),z_found(irec)
   enddo
   close(IIN_SU1)
 
   ! reads in source locations from output_list file
-  open(unit=IIN_SU1,file=trim(OUTPUT_FILES)//'/output_list_sources.txt',status='unknown',iostat=ier)
+  open(unit=IIN_SU1,file=trim(OUTPUT_FILES)//'/output_list_sources.txt',status='old',iostat=ier)
   if( ier /= 0 ) stop 'error opening output_list_sources.txt file'
 
   read(IIN_SU1,*) x_found_source,y_found_source,z_found_source