[cig-commits] [commit] master: bugfix in wavefield output (be819aa)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/9380934afd0d9a3e12a892767413896fb2fa7e60...be819aae57cc99f055989132fc778afa4778d3ec

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

commit be819aae57cc99f055989132fc778afa4778d3ec
Author: martinvandriel <martin at vandriel.de>
Date:   Tue Sep 16 18:48:32 2014 +0200

    bugfix in wavefield output


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

be819aae57cc99f055989132fc778afa4778d3ec
 SOLVER/wavefields_io.f90 | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/SOLVER/wavefields_io.f90 b/SOLVER/wavefields_io.f90
index 5c34287..93f01bc 100644
--- a/SOLVER/wavefields_io.f90
+++ b/SOLVER/wavefields_io.f90
@@ -923,6 +923,7 @@ subroutine dump_velo_global(v, dchi, istrain)
       ! tranform from +/- coordinates to s/phi
       f(:,:,:,1) = v(:,:,:,1) + v(:,:,:,2)
       f(:,:,:,2) = v(:,:,:,1) - v(:,:,:,2)
+      f(:,:,:,3) = v(:,:,:,3)
    else
       f = v
    endif
@@ -1009,6 +1010,7 @@ subroutine dump_disp_global(u, chi, istrain)
       ! tranform from +/- coordinates to s/phi
       f(:,:,:,1) = u(:,:,:,1) + u(:,:,:,2)
       f(:,:,:,2) = u(:,:,:,1) - u(:,:,:,2)
+      f(:,:,:,3) = u(:,:,:,3)
    else
       f = u
    endif