diff --git a/Manuals/FDS_User_Guide/FDS_User_Guide.tex b/Manuals/FDS_User_Guide/FDS_User_Guide.tex index cbc6940d2da..d425a9e4284 100644 --- a/Manuals/FDS_User_Guide/FDS_User_Guide.tex +++ b/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} \\ diff --git a/Source/read.f90 b/Source/read.f90 index dc1f7260033..b8142c5dd2a 100644 --- a/Source/read.f90 +++ b/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