@@ -64,8 +64,6 @@ SUBROUTINE INITIALIZE_LEVEL_SET_FIRESPREAD_1(NM)
6464ALLOCATE (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
6565ALLOCATE (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
6666
67- ALLOCATE (M% TOA(0 :IBP1,0 :JBP1)) ; CALL ChkMemErr(' VEGE:TOA' ,' TOA' ,IZERO) ; TOA = > M% TOA ; TOA = T_END + 1._EB
68-
6967! Wind speed components in the center of the first gas phsae cell above the ground.
7068
7169ALLOCATE (M% U_LS(0 :IBP1,0 :JBP1)) ; CALL ChkMemErr(' VEGE:LEVEL SET' ,' U_LS' ,IZERO) ; U_LS = > M% U_LS ; U_LS = 0._EB
@@ -516,7 +514,6 @@ SUBROUTINE LEVEL_SET_FIRESPREAD(T,DT,NM)
516514 CFA = > CFACE(CUT_FACE(ICF)% CFACE_INDEX(IW))
517515 B1 = > BOUNDARY_PROP1(CFA% B1_INDEX)
518516 IF (PHI_LS(IIG,JJG)>= 0._EB .AND. B1% T_IGN> 1.E5_EB ) CALL IGNITE_GRID_CELL
519- IF (PHI_LS(IIG,JJG)>= 0._EB .AND. TOA(IIG,JJG)>T_END) TOA(IIG,JJG) = T
520517 B2 = > BOUNDARY_PROP2(CFA% B2_INDEX)
521518 B2% PHI_LS = PHI_LS(IIG,JJG)
522519 ENDDO
@@ -534,7 +531,6 @@ SUBROUTINE LEVEL_SET_FIRESPREAD(T,DT,NM)
534531 WC = > WALL(IW)
535532 B1 = > BOUNDARY_PROP1(WC% B1_INDEX)
536533 IF (PHI_LS(IIG,JJG)>= 0._EB .AND. B1% T_IGN> 1.E5_EB ) CALL IGNITE_GRID_CELL
537- IF (PHI_LS(IIG,JJG)>= 0._EB .AND. TOA(IIG,JJG)>T_END) TOA(IIG,JJG) = T
538534 B2 = > BOUNDARY_PROP2(WC% B2_INDEX)
539535 B2% PHI_LS = PHI_LS(IIG,JJG)
540536 ENDDO
@@ -619,7 +615,7 @@ SUBROUTINE FILL_BOUNDARY_VALUES
619615
620616USE COMPLEX_GEOMETRY, ONLY : CC_CGSC,CC_SOLID,CC_CUTCFE
621617INTEGER :: IW,IIO,JJO,N_INT_CELLS,NOM,IC
622- REAL (EB) :: PHI_LS_OTHER,U_LS_OTHER,V_LS_OTHER,Z_LS_OTHER,TOA_LS_OTHER
618+ REAL (EB) :: PHI_LS_OTHER,U_LS_OTHER,V_LS_OTHER,Z_LS_OTHER
623619TYPE (EXTERNAL_WALL_TYPE), POINTER :: EWC
624620LOGICAL :: SOLID_CELL
625621
@@ -647,7 +643,6 @@ SUBROUTINE FILL_BOUNDARY_VALUES
647643U_LS_OTHER = 0._EB
648644V_LS_OTHER = 0._EB
649645Z_LS_OTHER = 0._EB
650- TOA_LS_OTHER = 0._EB
651646DO JJO= EWC% JJO_MIN,EWC% JJO_MAX
652647 DO IIO= EWC% IIO_MIN,EWC% IIO_MAX
653648 IF (PREDICTOR) THEN
@@ -658,7 +653,6 @@ SUBROUTINE FILL_BOUNDARY_VALUES
658653 U_LS_OTHER = U_LS_OTHER + OMESH(NOM)% U_LS(IIO,JJO)
659654 V_LS_OTHER = V_LS_OTHER + OMESH(NOM)% V_LS(IIO,JJO)
660655 Z_LS_OTHER = Z_LS_OTHER + OMESH(NOM)% Z_LS(IIO,JJO)
661- TOA_LS_OTHER = TOA_LS_OTHER + OMESH(NOM)% TOA(IIO,JJO)
662656 ENDDO
663657ENDDO
664658N_INT_CELLS = (EWC% IIO_MAX- EWC% IIO_MIN+1 ) * (EWC% JJO_MAX- EWC% JJO_MIN+1 )
@@ -669,7 +663,6 @@ SUBROUTINE FILL_BOUNDARY_VALUES
669663 U_LS(II,JJ) = U_LS_OTHER/ REAL (N_INT_CELLS,EB)
670664 V_LS(II,JJ) = V_LS_OTHER/ REAL (N_INT_CELLS,EB)
671665 Z_LS(II,JJ) = Z_LS_OTHER/ REAL (N_INT_CELLS,EB)
672- TOA(II,JJ) = TOA_LS_OTHER/ REAL (N_INT_CELLS,EB)
673666 CASE (3 ) ! only grab a PHI_LS value from the other mesh if the (II,JJ) cell of the current mesh has no terrain surface
674667 SOLID_CELL = .FALSE.
675668 IF (CC_IBM) THEN
@@ -682,7 +675,6 @@ SUBROUTINE FILL_BOUNDARY_VALUES
682675 U_LS(II,JJ) = U_LS_OTHER/ REAL (N_INT_CELLS,EB)
683676 V_LS(II,JJ) = V_LS_OTHER/ REAL (N_INT_CELLS,EB)
684677 Z_LS(II,JJ) = Z_LS_OTHER/ REAL (N_INT_CELLS,EB)
685- TOA(II,JJ) = TOA_LS_OTHER/ REAL (N_INT_CELLS,EB)
686678 ENDIF
687679 CASE (- 3 ) ! only grab a PHI_LS value from the other mesh if the (II,JJ) cell of the current mesh has no terrain surface
688680 SOLID_CELL = .FALSE.
@@ -696,7 +688,6 @@ SUBROUTINE FILL_BOUNDARY_VALUES
696688 U_LS(II,JJ) = U_LS_OTHER/ REAL (N_INT_CELLS,EB)
697689 V_LS(II,JJ) = V_LS_OTHER/ REAL (N_INT_CELLS,EB)
698690 Z_LS(II,JJ) = Z_LS_OTHER/ REAL (N_INT_CELLS,EB)
699- TOA(II,JJ) = TOA_LS_OTHER/ REAL (N_INT_CELLS,EB)
700691 ENDIF
701692END SELECT
702693
0 commit comments