[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
cig_noreply at geodynamics.org
Tue Aug 26 06:52:21 PDT 2014
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
More information about the CIG-COMMITS
mailing list