Skip to content

Commit c6b23f2

Browse files
committed
FDS Source: Ensure continuity of SM3D TEMPERATURE
1 parent 1480d4e commit c6b23f2

File tree

3 files changed

+2
-10
lines changed

3 files changed

+2
-10
lines changed

Source/dump.f90

Lines changed: 0 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -4428,7 +4428,6 @@ SUBROUTINE DUMP_SMOKE3D(T,DT,NM)
44284428
REAL(FB) :: DXX,STIME
44294429
REAL(EB), POINTER, DIMENSION(:,:,:) :: FF
44304430
REAL(FB), ALLOCATABLE, DIMENSION(:) :: QQ_PACK
4431-
REAL(EB) :: FR_C
44324431
TYPE(SMOKE3D_TYPE), POINTER :: S3
44334432

44344433
! Miscellaneous settings
@@ -4453,13 +4452,6 @@ SUBROUTINE DUMP_SMOKE3D(T,DT,NM)
44534452
ENDDO
44544453
ENDDO
44554454

4456-
! Adjust the temperature as it is used in the expression for the radiation source term
4457-
4458-
IF (S3%DISPLAY_TYPE=='TEMPERATURE' .AND. RTE_SOURCE_CORRECTION) THEN
4459-
FR_C = RTE_SOURCE_CORRECTION_FACTOR**0.25_EB
4460-
WHERE (CHI_R*Q>QR_CLIP) FF = (FF+TMPM)*FR_C - TMPM
4461-
ENDIF
4462-
44634455
! Interpolate data to cell nodes
44644456

44654457
DO K=0,KBAR

Source/read.f90

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -16026,7 +16026,7 @@ SUBROUTINE READ_SM3D
1602616026
IF (N<=N_SMOKE3D_RESERVED) THEN
1602716027
IF (N==1) THEN ; QUANTITY='DENSITY' ; SPEC_ID='SOOT' ; ENDIF
1602816028
IF (N==2) THEN ; QUANTITY='HRRPUV' ; ENDIF
16029-
IF (N==3) THEN ; QUANTITY='TEMPERATURE' ; ENDIF
16029+
IF (N==3) THEN ; QUANTITY='EFFECTIVE FLAME TEMPERATURE' ; ENDIF
1603016030
ELSE
1603116031
CALL CHECKREAD('SM3D',LU_INPUT,IOS) ; IF (STOP_STATUS==SETUP_STOP) RETURN
1603216032
IF (IOS==1) EXIT READ_SM3D_LOOP

Source/smvv.f90

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -979,7 +979,7 @@ SUBROUTINE SMOKE3D_TO_FILE(NM,TIME,DX,SMOKE3D_INDEX,VALS,NX,NY,NZ)
979979
MAX_VAL = MAX(VAL_FDS,MAX_VAL)
980980
ENDDO
981981

982-
ELSEIF (S3%DISPLAY_TYPE=='TEMPERATURE') THEN
982+
ELSEIF (S3%DISPLAY_TYPE=='TEMPERATURE' .OR. S3%DISPLAY_TYPE=='EFFECTIVE FLAME TEMPERATURE') THEN
983983

984984
DO I=1,NVALS
985985
VAL_FDS = MIN(TEMP_MAX_SMV,MAX(TEMP_MIN_SMV,VALS(I)))

0 commit comments

Comments
 (0)