Skip to content

FDS Source: fix runtime issue in Issue #15890 #15929

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
drjfloyd merged 1 commit into from
Feb 6, 2026
+5 −4
Merged

FDS Source: fix runtime issue in Issue #15890 #15929

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
9 changes: 5 additions & 4 deletions Source/func.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2813,20 +2813,21 @@ END SUBROUTINE GET_Y_SURF

SUBROUTINE GET_TMP_REF(N_MATL,NR)
INTEGER, INTENT(IN) :: N_MATL,NR
REAL(EB) :: HEATING_RATE,DT=0.01_EB,DTDT,RR_MAX,REACTION_RATE,TMP,RHO_S
REAL(EB) :: HEATING_RATE,DTMP=0.01_EB,DT,DTMPDT,RR_MAX,REACTION_RATE,TMP,RHO_S
TYPE(MATERIAL_TYPE), POINTER :: ML

ML=> MATERIAL(N_MATL)

IF (ML%RATE_REF(NR) > 0._EB) THEN
HEATING_RATE = ML%RATE_REF(NR)
HEATING_RATE = ML%HEATING_RATE(NR)
ELSE
HEATING_RATE = TGA_HEATING_RATE
ENDIF

TMP = 0._EB
ML%TMP_REF(NR) = 0._EB
DTDT = HEATING_RATE/60._EB
DTMPDT = HEATING_RATE/60._EB
DT = DTMP/DTMPDT
RHO_S = ML%RHO_S
IF (ABS(ML%E(NR)) < TWENTY_EPSILON_EB) THEN
RR_MAX = ML%A(NR)*RHO_S**ML%N_S(NR)
Expand All @@ -2835,7 +2836,7 @@ SUBROUTINE GET_TMP_REF(N_MATL,NR)
ENDIF

DO WHILE (INT(TMP)<I_MAX_TEMP)
TMP = TMP + DTDT * DT
TMP = TMP + DTMP
REACTION_RATE = ML%A(NR)*RHO_S**ML%N_S(NR)*EXP(-ML%E(NR)/(R0*TMP))*TMP**ML%N_T(NR)*ML%X_O2_PYRO**ML%N_O2(NR)
IF (REACTION_RATE > RR_MAX) THEN
ML%TMP_REF(NR) = TMP
Expand Down
Loading