Merge pull request #13226 from mcgratta/master

FDS Source: Issue #13224. Fix problem with cell detection
firemodels · Jul 26, 2024 · efe28a1 · efe28a1
2 parents 9e3cdc3 + 45764b0
commit efe28a1
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/Source/read.f90 b/Source/read.f90
@@ -1418,7 +1418,7 @@ SUBROUTINE READ_TRAN
    M%CELLSI_HI=NIPX+NIPXF
    DO I=M%CELLSI_LO,M%CELLSI_HI
       M%CELLSI(I) = GINV(REAL(I,EB)/M%RDXINT,1,NM)*M%RDXI
-      IF(MOD(I,500)==0) M%CELLSI(I) = REAL(NINT(M%CELLSI(I)),EB)
+      IF (ABS(M%CELLSI(I)-REAL(NINT(M%CELLSI(I)),EB)) < 1.E-14_EB) M%CELLSI(I) = REAL(NINT(M%CELLSI(I)),EB)
       M%CELLSI(I) = MAX(M%CELLSI(I),-0.9_EB)
       M%CELLSI(I) = MIN(M%CELLSI(I),REAL(M%IBAR)+0.9_EB)
    ENDDO
@@ -1427,7 +1427,7 @@ SUBROUTINE READ_TRAN
    M%CELLSJ_HI=NIPY+NIPYF
    DO J=M%CELLSJ_LO,M%CELLSJ_HI
       M%CELLSJ(J) = GINV(REAL(J,EB)/M%RDYINT,2,NM)*M%RDETA
-      IF(MOD(J,500)==0) M%CELLSJ(J) = REAL(NINT(M%CELLSJ(J)),EB)
+      IF (ABS(M%CELLSJ(J)-REAL(NINT(M%CELLSJ(J)),EB)) < 1.E-14_EB) M%CELLSJ(J) = REAL(NINT(M%CELLSJ(J)),EB)
       M%CELLSJ(J) = MAX(M%CELLSJ(J),-0.9_EB)
       M%CELLSJ(J) = MIN(M%CELLSJ(J),REAL(M%JBAR)+0.9_EB)
    ENDDO
@@ -1436,7 +1436,7 @@ SUBROUTINE READ_TRAN
    M%CELLSK_HI=NIPZ+NIPZF
    DO K=M%CELLSK_LO,M%CELLSK_HI
       M%CELLSK(K) = GINV(REAL(K,EB)/M%RDZINT,3,NM)*M%RDZETA
-      IF(MOD(K,500)==0) M%CELLSK(K) = REAL(NINT(M%CELLSK(K)),EB)
+      IF (ABS(M%CELLSK(K)-REAL(NINT(M%CELLSK(K)),EB)) < 1.E-14_EB) M%CELLSK(K) = REAL(NINT(M%CELLSK(K)),EB)
       M%CELLSK(K) = MAX(M%CELLSK(K),-0.9_EB)
       M%CELLSK(K) = MIN(M%CELLSK(K),REAL(M%KBAR)+0.9_EB)
    ENDDO