Merge pull request firemodels#11987 from ericvmueller/master

FDS Source: change undefined value for TIME OF ARRIVAL quantity to T_…
mcgratta · Jul 21, 2023 · 9310117 · 9310117
2 parents bfe3263 + def4b04
commit 9310117
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/Source/main.f90 b/Source/main.f90
@@ -1999,7 +1999,7 @@ SUBROUTINE INITIALIZE_MESH_EXCHANGE_1(NM)
    IF (LEVEL_SET_MODE>0 .OR. TERRAIN_CASE) THEN
       ALLOCATE(OM%PHI_LS(IMIN:IMAX,JMIN:JMAX))  ; OM%PHI_LS   = -1._EB
       ALLOCATE(OM%PHI1_LS(IMIN:IMAX,JMIN:JMAX)) ; OM%PHI1_LS  = -1._EB
-      ALLOCATE(OM%TOA(IMIN:IMAX,JMIN:JMAX))     ; OM%TOA      = -1._EB
+      ALLOCATE(OM%TOA(IMIN:IMAX,JMIN:JMAX))     ; OM%TOA      = T_END + 1._EB
       ALLOCATE(OM%U_LS(IMIN:IMAX,JMIN:JMAX))    ; OM%U_LS     =  0._EB
       ALLOCATE(OM%V_LS(IMIN:IMAX,JMIN:JMAX))    ; OM%V_LS     =  0._EB
       ALLOCATE(OM%Z_LS(IMIN:IMAX,JMIN:JMAX))    ; OM%Z_LS     =  0._EB

diff --git a/Source/vege.f90 b/Source/vege.f90
@@ -64,7 +64,7 @@ SUBROUTINE INITIALIZE_LEVEL_SET_FIRESPREAD_1(NM)
 ALLOCATE(M%PHI_LS(0:IBP1,0:JBP1)) ; CALL ChkMemErr('VEGE:LEVEL SET','PHI_LS',IZERO)   ; PHI_LS  => M%PHI_LS  ; PHI_LS = PHI_LS_MIN
 ALLOCATE(M%PHI1_LS(0:IBP1,0:JBP1)); CALL ChkMemErr('VEGE:LEVEL SET','PHI1_LS',IZERO)  ; PHI1_LS => M%PHI1_LS ; PHI1_LS = PHI_LS_MIN
 
-ALLOCATE(M%TOA(0:IBP1,0:JBP1)) ; CALL ChkMemErr('VEGE:TOA','TOA',IZERO)   ; TOA  => M%TOA  ; TOA = -1._EB
+ALLOCATE(M%TOA(0:IBP1,0:JBP1)) ; CALL ChkMemErr('VEGE:TOA','TOA',IZERO)   ; TOA  => M%TOA  ; TOA = T_END + 1._EB
 
 ! Wind speed components in the center of the first gas phsae cell above the ground.
 
@@ -512,7 +512,7 @@ SUBROUTINE LEVEL_SET_FIRESPREAD(T,DT,NM)
                   CFA => CFACE(CUT_FACE(ICF)%CFACE_INDEX(IW))
                   B1  => BOUNDARY_PROP1(CFA%B1_INDEX)
                   IF (PHI_LS(IIG,JJG)>=0._EB .AND. B1%T_IGN>1.E5_EB) CALL IGNITE_GRID_CELL
-                  IF (PHI_LS(IIG,JJG)>=0._EB .AND. TOA(IIG,JJG)<0._EB) TOA(IIG,JJG) = T
+                  IF (PHI_LS(IIG,JJG)>=0._EB .AND. TOA(IIG,JJG)>T_END) TOA(IIG,JJG) = T
                   B2 => BOUNDARY_PROP2(CFA%B2_INDEX)
                   B2%PHI_LS = PHI_LS(IIG,JJG)
                ENDDO
@@ -530,7 +530,7 @@ SUBROUTINE LEVEL_SET_FIRESPREAD(T,DT,NM)
             WC => WALL(IW)
             B1 => BOUNDARY_PROP1(WC%B1_INDEX)
             IF (PHI_LS(IIG,JJG)>=0._EB .AND. B1%T_IGN>1.E5_EB) CALL IGNITE_GRID_CELL
-            IF (PHI_LS(IIG,JJG)>=0._EB .AND. TOA(IIG,JJG)<0._EB) TOA(IIG,JJG) = T
+            IF (PHI_LS(IIG,JJG)>=0._EB .AND. TOA(IIG,JJG)>T_END) TOA(IIG,JJG) = T
             B2 => BOUNDARY_PROP2(WC%B2_INDEX)
             B2%PHI_LS = PHI_LS(IIG,JJG)
          ENDDO