Skip to content

FDS Source: Add LOGTEN = LOG(10._EB) to prec #16112

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
Apr 2, 2026
+5 −5
Merged

FDS Source: Add LOGTEN = LOG(10._EB) to prec #16112

Show file tree
Hide file tree
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
7 changes: 3 additions & 4 deletions Source/chem.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -710,13 +710,12 @@ END SUBROUTINE CALC_FALLOFF_DBIDC_AND_DBIDT
REAL(EB) FUNCTION DDC_TROE(P_RATIO, F, DPRDC, TMP, RNI)
REAL(EB), INTENT(IN) :: P_RATIO, F, DPRDC, TMP
INTEGER, INTENT(IN) :: RNI
REAL(EB) :: LOGPR, LOGTEN, LOGFCENT, C, N, DLOGPRDC, DPARENTDC
REAL(EB) :: LOGPR, LOGFCENT, C, N, DLOGPRDC, DPARENTDC
TYPE(REACTION_TYPE), POINTER :: RN
REAL(EB), PARAMETER :: D=0.14_EB

RN => REACTION(RNI)
LOGPR = LOG10(MAX(P_RATIO, TWO_EPSILON_EB))
LOGTEN = LOG(10.0)
IF (RN%T2_TROE <-1.E20_EB) THEN
LOGFCENT = LOG10(MAX((1 - RN%A_TROE)*EXP(-TMP*RN%RT3_TROE) + &
RN%A_TROE*EXP(-TMP*RN%RT1_TROE),TWO_EPSILON_EB))
Expand Down Expand Up @@ -747,13 +746,13 @@ END FUNCTION DDC_TROE
REAL(EB) FUNCTION DDTMP_TROE(P_RATIO, F, DPRDT, TMP, RNI)
REAL(EB), INTENT(IN) :: P_RATIO, F, DPRDT, TMP
INTEGER, INTENT(IN) :: RNI
REAL(EB) :: FCENT, LOGPR, LOGTEN, LOGFCENT, DFCENTDT, C, N, DCDT, DNDT, DPARENTDT, DLOGFCENTDT, DLOGPRDT
REAL(EB) :: FCENT, LOGPR, LOGFCENT, DFCENTDT, C, N, DCDT, DNDT, DPARENTDT, DLOGFCENTDT, DLOGPRDT
TYPE(REACTION_TYPE), POINTER :: RN
REAL(EB), PARAMETER :: D=0.14_EB

RN => REACTION(RNI)
LOGPR = LOG10(MAX(P_RATIO, TWO_EPSILON_EB));
LOGTEN = LOG(10.0);

IF (RN%T2_TROE <-1.E20_EB) THEN
FCENT = MAX((1 - RN%A_TROE)*EXP(-TMP*RN%RT3_TROE) + &
RN%A_TROE*EXP(-TMP*RN%RT1_TROE),TWO_EPSILON_EB)
Expand Down
2 changes: 1 addition & 1 deletion Source/dump.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -7825,7 +7825,7 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
ELSEIF (Y_INDEX>0) THEN
MEC = SPECIES(Y_INDEX)%MASS_EXTINCTION_COEFFICIENT
ENDIF
GAS_PHASE_OUTPUT_RES = Y_SPECIES*RHO(II,JJ,KK)*MEC/2.3_EB
GAS_PHASE_OUTPUT_RES = Y_SPECIES*RHO(II,JJ,KK)*MEC/LOGTEN

CASE(99) ! PRESSURE POISSON RESIDUAL
GAS_PHASE_OUTPUT_RES = PP_RESIDUAL(II,JJ,KK)
Expand Down
1 change: 1 addition & 0 deletions Source/prec.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -52,5 +52,6 @@ MODULE PRECISION_PARAMETERS
RFPI=1._EB/(4._EB*PI), FOTHPI = FOTH*PI, CR2=2._EB**(1._EB/3._EB)
INTEGER, PARAMETER :: INTEGER_ZERO=0,INTEGER_ONE=1,INTEGER_TWO=2,INTEGER_THREE=3
REAL(EB), PARAMETER :: DEG2RAD=PI/180.0_EB
REAL(EB), PARAMETER :: LOGTEN = LOG(10._EB)

END MODULE PRECISION_PARAMETERS
Loading