[cig-commits] [commit] master: Bugfix: Number of seismogram dumps was too small (46eb20d)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/4fdc3be0b9163702cd293f978f30a53c8313fa17...7b7289c7dbd6c6f7e20177c657193146d6e0b70e

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

commit 46eb20d1c57e947fdb3fff42c0ae04bc4fd03144
Author: Simon Stähler <staehler at geophysik.uni-muenchen.de>
Date:   Mon Aug 25 16:48:07 2014 +0200

    Bugfix: Number of seismogram dumps was too small
    
     - Led to subtle boundary violation in NetCDF output, where first sample of seismograms was overwritten by value of the last from another station


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

46eb20d1c57e947fdb3fff42c0ae04bc4fd03144
 SOLVER/parameters.F90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/SOLVER/parameters.F90 b/SOLVER/parameters.F90
index df451f2..3928cff 100644
--- a/SOLVER/parameters.F90
+++ b/SOLVER/parameters.F90
@@ -887,7 +887,7 @@ subroutine compute_numerical_parameters
   endif
   deltat_coarse = seis_dt
 
-  nseismo = ceiling(real(niter) / real(seis_it)) 
+  nseismo = floor(real(niter) / real(seis_it)) + 1
 
   ! Frequency of checkpointing. Hardcoded to every 5% of runtime
   check_it = niter / 20