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