SINGLE_PRECISION bug

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
jcwarner
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

SINGLE_PRECISION bug

#1 Unread post by jcwarner »

small but tricky bug when #define SINGLE_PRECISION
1) in get_idata we have
real(r8) :: time_save = 0.0_r8
but it needs to be
real(dp) :: time_save = 0.0_dp

2) And then there are multiple uses of:

time_save=time(ng)
time(ng)=8640000.0_r8
.... call get2d..
time(ng)=time_save

I think this should be :

time_save=time(ng)
time(ng)=8640000.0_dp
.... call get2d..
time(ng)=time_save

but maybe the r8 is ok.
-j

Post Reply