Merge pull request #15503 from drjfloyd/master

FDS Source: Add error for nonsensical E value (Issue #15497).

Author: drjfloyd

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -14190,6 +14190,7 @@ \chapter{Error Codes}
 266  \> \ct{REAC ... FUEL ... is not a predefined or tracked species.} \> Section~\ref{info:simple_chemistry}  \\
 267  \> \ct{REAC ... has no consumed species.} \> Section~\ref{info:REAC_Diagnostics}  \\
 268  \> \ct{BACK tracked species ... not found.} \> Section~\ref{info:BACK}  \\
+269  \> \ct{MATL ... The specified pyrolysis parameters result in A >...} \> Section~\ref{info:kinetic_parameters}  \\
      \>                                                                               \>                                     \\
 299  \> \ct{SURF ... has no solid or particle for TGA.} \> Section~\ref{info:TGA_DSC_MCC} \\
 300  \> \ct{N_LAYER_CELLS_MAX should be at least ... for ...} \> Section~\ref{info:solid_phase_stability} \\

Source/read.f90

@@ -7347,6 +7347,12 @@ SUBROUTINE READ_MATL
             PEAK_REACTION_RATE = 2._EB*ML%HEATING_RATE(NR)/ML%PYROLYSIS_RANGE(NR)
          ENDIF
          ML%E(NR) = EXP(1._EB)*PEAK_REACTION_RATE*R0*ML%TMP_REF(NR)**2/ML%HEATING_RATE(NR)
+         ! 0.0001 HUGE_EB is used so the reaction rate is not overflowed when multiplied by RHO.
+         IF (LOG(0.0001_EB*HUGE_EB)*R0*ML%TMP_REF(NR) < ML%E(NR)) THEN
+            WRITE(MESSAGE,'(3A,E10.3)') 'ERROR(269): MATL ',TRIM(ML%ID),' The specified pyrolysis parameters result in A >',&
+               0.0001_EB*HUGE_EB
+            CALL SHUTDOWN(MESSAGE) ; RETURN
+         ENDIF
          ML%A(NR) = EXP(1._EB)*PEAK_REACTION_RATE*EXP(ML%E(NR)/(R0*ML%TMP_REF(NR)))
       ENDIF