[cig-commits] [commit] master: round to nearest integer instead of down. Do to floating point inaccuracy, this led to (42a3e78)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/axisem/compare/7b7289c7dbd6c6f7e20177c657193146d6e0b70e...42a3e78e1630c29b5752d3fcfd65c0ffbe45d070

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

commit 42a3e78e1630c29b5752d3fcfd65c0ffbe45d070
Author: martinvandriel <martin at vandriel.de>
Date:   Tue Aug 26 15:51:40 2014 +0200

    round to nearest integer instead of down. Do to floating point inaccuracy, this led to
    
    wrong numbers before


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

42a3e78e1630c29b5752d3fcfd65c0ffbe45d070
 SOLVER/parameters.F90 | 12 ++++++------
 1 file changed, 6 insertions(+), 6 deletions(-)

diff --git a/SOLVER/parameters.F90 b/SOLVER/parameters.F90
index 3928cff..ceb1a41 100644
--- a/SOLVER/parameters.F90
+++ b/SOLVER/parameters.F90
@@ -1392,9 +1392,9 @@ subroutine write_parameters
         write(55,22)ibeg,'  ibeg: beginning gll index for wavefield dumps'
         write(55,22)iend,'iend: end gll index for wavefield dumps'
         write(55,21)shift_fact,'source shift factor [s]'
-        write(55,22)int(shift_fact/deltat),'source shift factor for deltat'
-        write(55,22)int(shift_fact/seis_dt),'source shift factor for seis_dt'
-        write(55,22)int(shift_fact/deltat_coarse),'source shift factor for deltat_coarse'
+        write(55,22)nint(shift_fact/deltat),'source shift factor for deltat'
+        write(55,22)nint(shift_fact/seis_dt),'source shift factor for seis_dt'
+        write(55,22)nint(shift_fact/deltat_coarse),'source shift factor for deltat_coarse'
         write(55,23)trim(rec_file_type),'receiver file type'
         write(55,21)dtheta_rec,'receiver spacing (0 if not even)'
         write(55,24)use_netcdf,'use netcdf for wavefield output?'
@@ -1465,9 +1465,9 @@ subroutine write_parameters
         call nc_write_att_int(  ibeg,                  'ibeg')
         call nc_write_att_int(  iend,                  'iend')
         call nc_write_att_real( shift_fact,            'source shift factor in sec')
-        call nc_write_att_int(  int(shift_fact/deltat),  'source shift factor for deltat')
-        call nc_write_att_int(  int(shift_fact/seis_dt), 'source shift factor for seis_dt')
-        call nc_write_att_int(  int(shift_fact/deltat_coarse), 'source shift factor for deltat_coarse')
+        call nc_write_att_int(  nint(shift_fact/deltat),  'source shift factor for deltat')
+        call nc_write_att_int(  nint(shift_fact/seis_dt), 'source shift factor for seis_dt')
+        call nc_write_att_int(  nint(shift_fact/deltat_coarse), 'source shift factor for deltat_coarse')
         call nc_write_att_char( trim(rec_file_type),   'receiver file type')
         call nc_write_att_real( dtheta_rec,            'receiver spacing (0 if not even)')
         write(clogic,*) use_netcdf