FDS Source: Initialize AIT zone temp Issue #15175

Author: drjfloyd

Source/read.f90

@@ -4667,7 +4667,7 @@ SUBROUTINE READ_REAC
 ALLOCATE(REACTION(N_REACTIONS+N_REVERSE),STAT=IZERO)
 ALLOCATE(SIMPLE_FUEL_DEFINED(N_REACTIONS),STAT=IZERO)
 SIMPLE_FUEL_DEFINED = .FALSE.
-
+WRITE(LU_ERR,*) 'RNA:',REACTION(1)%AIT_EXCLUSION_ZONE(1)%AIT
 ! Read and store the reaction parameters
 N_REVERSE = 0
 REAC_READ_LOOP: DO NR=1,N_REACTIONS
@@ -4790,6 +4790,7 @@ SUBROUTINE READ_REAC
          RN%AIT_EXCLUSION_ZONE(IZ)%CTRL_ID = AIT_EXCLUSION_ZONE_CTRL_ID(IZ)
       ENDIF
    ENDDO
+WRITE(LU_ERR,*) 'RNB:',REACTION(1)%AIT_EXCLUSION_ZONE(1)%AIT,AIT_EXCLUSION_ZONE_TEMPERATURE(1)
    RN%E                         = E*1000._EB
    RN%E_IN                      = E
    RN%EQUATION                  = EQUATION
@@ -5034,6 +5035,7 @@ SUBROUTINE READ_REAC
       CALL SHUTDOWN(MESSAGE) ; RETURN
    ENDIF
 ENDDO
+WRITE(LU_ERR,*) 'RNC:',REACTION(1)%AIT_EXCLUSION_ZONE(1)%AIT
 
 DEALLOCATE (REAC_FUEL_2)
 
@@ -5046,6 +5048,7 @@ SUBROUTINE SET_REAC_DEFAULTS
 A                           = -1._EB
 AUTO_IGNITION_TEMPERATURE   = -TMPM
 AIT_EXCLUSION_ZONE          = -1.E6_EB
+AIT_EXCLUSION_ZONE_TEMPERATURE  = -TMPM
 AIT_EXCLUSION_ZONE_CTRL_ID  = 'null'
 AIT_EXCLUSION_ZONE_DEVC_ID  = 'null'
 C                           = 0._EB
@@ -5542,13 +5545,13 @@ SUBROUTINE PROC_REAC_1
       ELSEIF(TRIM(SPECIES_MIXTURE(NS)%ID) == "O2" .OR. TRIM(SPECIES_MIXTURE(NS)%ID) == "OXYGEN") THEN
          I_O2 = NS
       ELSEIF(TRIM(SPECIES_MIXTURE(NS)%ID) == "CO2" .OR. TRIM(SPECIES_MIXTURE(NS)%ID) == "CARBON DIOXIDE") THEN
-         I_CO2 = NS   
+         I_CO2 = NS
       ELSEIF(TRIM(SPECIES_MIXTURE(NS)%ID) == "H2O" .OR. TRIM(SPECIES_MIXTURE(NS)%ID) == "WATER VAPOR") THEN
-         I_H2O = NS        
-      ENDIF   
+         I_H2O = NS
+      ENDIF
    ENDDO
 
-ENDIF  
+ENDIF
 
 END SUBROUTINE PROC_REAC_1
 
@@ -6605,7 +6608,7 @@ SUBROUTINE READ_PROP
          PY%SPECIFIC_HEAT(J) = 451.5_EB+0.001_EB*(REAL(MAX(293,MIN(1493,J)),EB)-293._EB)*(601.0_EB-451.5_EB)
       ENDDO
    ENDIF
-   PY%SPECIFIC_HEAT = PY%SPECIFIC_HEAT/TIME_SHRINK_FACTOR  
+   PY%SPECIFIC_HEAT = PY%SPECIFIC_HEAT/TIME_SHRINK_FACTOR
    IF (SPECIFIC_HEAT_RAMP /= 'null') CALL GET_RAMP_INDEX(SPECIFIC_HEAT_RAMP,'TEMPERATURE',PY%SPECIFIC_HEAT_RAMP_INDEX)
    PY%C_FACTOR                 = C_FACTOR
    PY%CHARACTERISTIC_VELOCITY  = CHARACTERISTIC_VELOCITY
@@ -6852,7 +6855,7 @@ SUBROUTINE SET_PROP_DEFAULTS
 DIAMETER                 = 0.001       ! m
 EMISSIVITY               = 0.85_EB
 HEAT_TRANSFER_COEFFICIENT= -1._EB      ! W/m2/K
-SPECIFIC_HEAT            = -1._EB      ! kJ/kg/K 
+SPECIFIC_HEAT            = -1._EB      ! kJ/kg/K
 SPECIFIC_HEAT_RAMP       = 'null'
 C_FACTOR                 = 0.0_EB
 CHARACTERISTIC_VELOCITY  = 1.0_EB      ! m/s
@@ -7019,7 +7022,7 @@ SUBROUTINE PROC_PROP
          PY%V_FACTOR = PY%PARTICLE_VELOCITY/SQRT(PY%OPERATING_PRESSURE)
       ENDIF
    ENDIF
-   
+
    ! Check units of specific heat
 
    IF (PY%SPECIFIC_HEAT_RAMP_INDEX > 0) THEN
@@ -8857,8 +8860,8 @@ SUBROUTINE READ_SURF(QUICK_READ)
       IF (DELAMINATION_TMP(NL) > 0._EB) SF%DELAMINATION_TMP(NL) = DELAMINATION_TMP(NL) + TMPM
       IF (DELAMINATION_DENSITY(NL) > 0._EB) SF%DELAMINATION_DENSITY(NL) = DELAMINATION_DENSITY(NL)
    ENDDO
-   
-      
+
+
    ! Adjust min and max temperature and convert C to K
 
    IF (TMP_FRONT >= -TMPM) TMPMIN = MIN(TMPMIN,TMP_FRONT+TMPM)
@@ -13700,7 +13703,7 @@ SUBROUTINE READ_DEVC
    IF (PROP_ID/='null') THEN
       CALL GET_PROPERTY_INDEX(PROP_INDEX,'DEVC',PROP_ID)
       IF (PROPERTY(PROP_INDEX)%IGNITION_ZONE) XB_XYZ= .FALSE.
-   ENDIF   
+   ENDIF
    IF (POINTS==1 .AND. SPATIAL_STATISTIC=='null' .AND. ALL(XB>-1E6_EB) .AND. ALL(XYZ<=-1E6_EB) .AND. XB_XYZ) THEN
       XYZ(1) = 0.5_EB*(XB(1)+XB(2))
       XYZ(2) = 0.5_EB*(XB(3)+XB(4))
@@ -14929,7 +14932,7 @@ SUBROUTINE PROC_DEVC
       WRITE(MESSAGE,'(3A)')  'ERROR(899): DEVC ',TRIM(DV%ID),' should not contain a + sign.'
       CALL SHUTDOWN(MESSAGE) ; RETURN
    ENDIF
-   
+
    ! Check for HVAC outputs with no HVAC inputs
 
    IF ((DV%DUCT_ID/='null' .OR. DV%NODE_ID(1)/='null') .AND. .NOT. HVAC_SOLVE) THEN
@@ -15006,7 +15009,7 @@ SUBROUTINE PROC_DEVC
          ENDIF
          IGNITION_ZONES(N_IGNITION_ZONES)%DEVC_INDEX = N
          IGNITION_ZONES(N_IGNITION_ZONES)%DEVC_ID = DV%ID
-      ENDIF   
+      ENDIF
    ENDIF
 
    IF(TRIM(ODE_SOLVER)=='CVODE' .AND. N_IGNITION_ZONES >= 1) THEN
@@ -15526,7 +15529,7 @@ SUBROUTINE READ_PROF
          CALL SHUTDOWN(MESSAGE) ; RETURN
       ENDIF
    ENDIF
-   
+
    IF (INIT_ID/='null') THEN
       INIT_DO: DO I=1,N_INIT
          IF (INIT_ID==INITIALIZATION(I)%ID) THEN