[cig-commits] [commit] devel: fixed reading adjoint traces for SH simulations (0eadab6)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem2d

On branch  : devel
Link       : https://github.com/geodynamics/specfem2d/compare/7a7ca08afb0578095ac833a886583b5fdee17ee3...caf697fd40a47a7cef1cbc80e13c600c3ec41167

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

commit 0eadab603028bb45faa02407fb005326f2e8d31d
Author: rmodrak <rmodrak at princeton.edu>
Date:   Thu Nov 13 16:51:58 2014 -0500

    fixed reading adjoint traces for SH simulations


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

0eadab603028bb45faa02407fb005326f2e8d31d
 src/specfem2D/specfem2D.F90 | 6 ++++++
 1 file changed, 6 insertions(+)

diff --git a/src/specfem2D/specfem2D.F90 b/src/specfem2D/specfem2D.F90
index 0081d13..f416d37 100644
--- a/src/specfem2D/specfem2D.F90
+++ b/src/specfem2D/specfem2D.F90
@@ -2377,6 +2377,9 @@
        write(filename, "('./SEM/Ux_file_single.bin.adj')")
        open(111,file=trim(filename),access='direct',recl=240+4*NSTEP,iostat = ios)
                if (ios /= 0) call exit_MPI(' file '//trim(filename)//'does not exist')
+       write(filename, "('./SEM/Uy_file_single.bin.adj')")
+       open(112,file=trim(filename),access='direct',recl=240+4*NSTEP,iostat = ios)
+               if (ios /= 0) call exit_MPI(' file '//trim(filename)//'does not exist')
        write(filename, "('./SEM/Uz_file_single.bin.adj')")
        open(113,file=trim(filename),access='direct',recl=240+4*NSTEP,iostat = ios)
                if (ios /= 0) call exit_MPI(' file '//trim(filename)//'does not exist')
@@ -2389,6 +2392,8 @@
           adj_sourcearray(:,:,:,:) = 0.0
           read(111,rec=irec,iostat=ios) r4head, adj_src_s(:,1)
                if (ios /= 0) call exit_MPI(' file '//trim(filename)//' read error')
+          read(112,rec=irec,iostat=ios) r4head, adj_src_s(:,2)
+               if (ios /= 0) call exit_MPI(' file '//trim(filename)//' read error')
           read(113,rec=irec,iostat=ios) r4head, adj_src_s(:,3)
                if (ios /= 0) call exit_MPI(' file '//trim(filename)//' read error')
           header2=int(r4head(29), kind=2)
@@ -2404,6 +2409,7 @@
          endif !  if(myrank == which_proc_receiver(irec))
        enddo ! irec
        close(111)
+       close(112)
        close(113)
        deallocate(adj_src_s)
     endif