[cig-commits] [commit] master: strain_info.dat???? should not be written by default (f021186)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/201b2a8c3949bb5acd44c5b1f0eeaba0012815d3...d9174f4e9c01eb069c0079f9a67d155de1997361

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

commit f021186cd81d75a64539f35eda9fbf00ba87e6a9
Author: Simon Stähler <staehler at geophysik.uni-muenchen.de>
Date:   Wed Nov 12 13:29:28 2014 +0100

    strain_info.dat???? should not be written by default
    
     - This file is useful only for debugging and several thousand of it cause trouble on parallel file systems


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

f021186cd81d75a64539f35eda9fbf00ba87e6a9
 SOLVER/parameters.F90 | 16 ++++++++++------
 1 file changed, 10 insertions(+), 6 deletions(-)

diff --git a/SOLVER/parameters.F90 b/SOLVER/parameters.F90
index 9951704..e4af391 100644
--- a/SOLVER/parameters.F90
+++ b/SOLVER/parameters.F90
@@ -1061,12 +1061,16 @@ subroutine compute_numerical_parameters
   if (dump_wavefields) then
      nstrain = floor(real(niter)/real(strain_it)) + 1
 
-     open(unit=2900+mynum,file=datapath(1:lfdata)//'/strain_info.dat'//appmynum)
-     write(2900+mynum,*) nstrain 
-     do ielem = 1, nstrain 
-        write(2900+mynum,*)real(ielem)*t_0/real(strain_samp),ielem*strain_it
-     enddo
-     close(2900+mynum)
+     ! This causes problems for massively parallel jobs on parallel file systems 
+     ! GPFS allows only a few hundreds of files in one directory
+     if (diagfiles) then
+        open(unit=2900+mynum,file=datapath(1:lfdata)//'/strain_info.dat'//appmynum)
+        write(2900+mynum,*) nstrain 
+        do ielem = 1, nstrain 
+           write(2900+mynum,*)real(ielem)*t_0/real(strain_samp),ielem*strain_it
+        enddo
+        close(2900+mynum)
+     end if
 
      if (lpr) then
         write(6,*)