[cig-commits] [commit] master: need these arrays in strain_only mode (bcc4f79)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/607f803cf074063627513d235f9ed0837fc1dd44...b6457db24acdde4a4e1c08935ae1b22adf87f5bf

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

commit bcc4f7977a6475740fef3812cd6cfe94dd9554bc
Author: martinvandriel <martin at vandriel.de>
Date:   Wed Oct 15 20:31:08 2014 +0200

    need these arrays in strain_only mode
    
    strange that this did not segfault with gfrotran, while intel now does


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

bcc4f7977a6475740fef3812cd6cfe94dd9554bc
 SOLVER/def_grid.f90 | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/SOLVER/def_grid.f90 b/SOLVER/def_grid.f90
index 85a1946..7aa538f 100644
--- a/SOLVER/def_grid.f90
+++ b/SOLVER/def_grid.f90
@@ -245,7 +245,8 @@ subroutine deallocate_preloop_arrays
   endif
 
   ! These terms are needed to compute the gradient!
-  if (.not. dump_wavefields .or. dump_type/='fullfields') then
+  if (.not. dump_wavefields .or. &
+        (dump_type /= 'fullfields' .and. dump_type /= 'strain_only')) then
      if (.not. anel_true .and. .not. dump_snaps_glob) then
         if (lpr) write(6,*)'  deallocating pointwise solid arrays...'
         deallocate(DsDeta_over_J_sol)