Fix an oddball result with i*2 numbers: abs(-32768)=-32768.

2020-09-12 09:00:39 -04:00 · 2020-09-12 09:00:39 -04:00 · 263675cac4
parent 14dad11f2e
commit 263675cac4
1 changed files with 2 additions and 0 deletions
--- a/lib/blanker.f90
+++ b/lib/blanker.f90
@ -8,6 +8,8 @@ subroutine blanker(iwave,nz,ndropmax,npct,c_bigfft)
  fblank=0.01*npct
  hist=0
  do i=1,nz
+! ### NB: if iwave(i)=-32768, abs(iwave(i))=-32768 ###
+     if(iwave(i).eq.-32768) iwave(i)=-32767
     n=abs(iwave(i))
     hist(n)=hist(n)+1
  enddo