updated ssm1.ft - floating-invalid fixed in SSM1OT

Fixed floating-invalid error that could arise when ssg(1,iGroup) > 2.1E+09 Replaced with suggested fix from Vivek Bedekar
wjqiu · Mar 21, 2019 · 01ff722 · 01ff722
1 parent 791d0ea
commit 01ff722
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/src/ssm1.f b/src/ssm1.f
@@ -1477,7 +1477,9 @@ SUBROUTINE SSM1OT(KPER,KSTP,NTRANS,TIME2)
         iGroup=ss(7,num)
         IF(iGroup.le.0) CYCLE
         ctmp=ssg(4,iGroup)
-        iFlag=int(ssg(1,iGroup))
+        !iFlag=int(ssg(1,iGroup))
+        iFlag = 0
+        IF(ABS((-999.)-ssg(1,iGroup)) .LE. 1.0E-03) iFlag=-999
         IF(iFlag.ne.-999) THEN
           ssg(1,iGroup)=-999
           WRITE(IU,1004) kper,kstp,ntrans,time2,iGroup,k,i,j,ctmp