Merge pull request firemodels#12679 from ericvmueller/master

FDS Source: prevent floating point error in value of M%CELLSI/J/K at …

Author: ericvmueller

Source/read.f90

@@ -1417,6 +1417,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)
       M%CELLSI(I) = MAX(M%CELLSI(I),-0.9_EB)
       M%CELLSI(I) = MIN(M%CELLSI(I),REAL(M%IBAR)+0.9_EB)
    ENDDO
@@ -1425,6 +1426,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)
       M%CELLSJ(J) = MAX(M%CELLSJ(J),-0.9_EB)
       M%CELLSJ(J) = MIN(M%CELLSJ(J),REAL(M%JBAR)+0.9_EB)
    ENDDO
@@ -1433,6 +1435,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)
       M%CELLSK(K) = MAX(M%CELLSK(K),-0.9_EB)
       M%CELLSK(K) = MIN(M%CELLSK(K),REAL(M%KBAR)+0.9_EB)
    ENDDO
@@ -6264,7 +6267,7 @@ SUBROUTINE READ_PROP
    IF (PY%VIEW_ANGLE<TWO_EPSILON_EB .OR. PY%VIEW_ANGLE>180._EB) THEN
       WRITE(MESSAGE,'(A,A,A)') 'ERROR(xxx): PROP ',TRIM(PY%ID),' VIEW_ANGLE must be between 0 and 180.'
       CALL SHUTDOWN(MESSAGE) ; RETURN
-   ENDIF  
+   ENDIF
 
 ENDDO READ_PROP_LOOP
 
@@ -13942,10 +13945,10 @@ SUBROUTINE PROC_DEVC
       ENDIF
 
       IF (DV%QUANTITY(1)=='null' .AND. PROPERTY(DV%PROP_INDEX)%QUANTITY/='null') DV%QUANTITY(1) = PROPERTY(DV%PROP_INDEX)%QUANTITY
-      
+
       IF (DV%PROP_INDEX > 0) THEN
          IF (PROPERTY(DV%PROP_INDEX)%VIEW_ANGLE < 180._EB) &
-            ORIENTATION_VIEW_ANGLE(DV%ORIENTATION_INDEX) = COS(PROPERTY(DV%PROP_INDEX)%VIEW_ANGLE/360._EB * PI) 
+            ORIENTATION_VIEW_ANGLE(DV%ORIENTATION_INDEX) = COS(PROPERTY(DV%PROP_INDEX)%VIEW_ANGLE/360._EB * PI)
       ENDIF
       ! Create an auto-ignition exclusion zone (AIT) in the cell containing a SPARK