From 1e40e25d183f8ede1c2bd8f4604e41606be76506 Mon Sep 17 00:00:00 2001 From: mcgratta Date: Tue, 13 Aug 2024 11:44:57 -0400 Subject: [PATCH] FDS Source: Issue #13235. Allow printout of Y+ and U_tau --- Source/cons.f90 | 1 + Source/read.f90 | 2 ++ Source/wall.f90 | 3 ++- 3 files changed, 5 insertions(+), 1 deletion(-) diff --git a/Source/cons.f90 b/Source/cons.f90 index 23359f5b413..726a8739294 100644 --- a/Source/cons.f90 +++ b/Source/cons.f90 @@ -268,6 +268,7 @@ MODULE GLOBAL_CONSTANTS LOGICAL :: CHECK_BOUNDARY_ONE_D_ARRAYS=.FALSE. !< Flag that indicates that ONE_D array dimensions need to be checked LOGICAL :: TENSOR_DIFFUSIVITY=.FALSE. !< If true, use experimental tensor diffusivity model for spec and tmp LOGICAL :: TEST_CHAR_MASS_TRANSFER_MODEL=.FALSE. !< Experimental flag to test mass transfer resistence in char model +LOGICAL :: OUTPUT_WALL_QUANTITIES=.FALSE. !< Flag to force call to WALL_MODEL INTEGER, ALLOCATABLE, DIMENSION(:) :: CHANGE_TIME_STEP_INDEX !< Flag to indicate if a mesh needs to change time step INTEGER, ALLOCATABLE, DIMENSION(:) :: SETUP_PRESSURE_ZONES_INDEX !< Flag to indicate if a mesh needs to keep searching for ZONEs diff --git a/Source/read.f90 b/Source/read.f90 index e5acec74978..31710d0413a 100644 --- a/Source/read.f90 +++ b/Source/read.f90 @@ -15835,6 +15835,8 @@ SUBROUTINE GET_QUANTITY_INDEX(SMOKEVIEW_LABEL,SMOKEVIEW_BAR_LABEL,OUTPUT_INDEX,O ENDIF ENDIF +IF (QUANTITY=='VISCOUS WALL UNITS' .OR. QUANTITY=='FRICTION VELOCITY') OUTPUT_WALL_QUANTITIES = .TRUE. + ! Assigne MATL_INDEX when MATL_ID is specified IF (MATL_ID/='null') THEN diff --git a/Source/wall.f90 b/Source/wall.f90 index 2aca081ae22..e54a85035d9 100644 --- a/Source/wall.f90 +++ b/Source/wall.f90 @@ -111,7 +111,8 @@ SUBROUTINE WALL_BC(T,DT,NM) CALL CALCULATE_RHO_D_F(B1,BC,WALL_INDEX=IW) ENDIF - IF (WC%BOUNDARY_TYPE==SOLID_BOUNDARY .AND. (ANY(SPECIES_MIXTURE%CONDENSATION_SMIX_INDEX>0).OR.DEPOSITION)) THEN + IF (WC%BOUNDARY_TYPE==SOLID_BOUNDARY .AND. & + (ANY(SPECIES_MIXTURE%CONDENSATION_SMIX_INDEX>0) .OR. DEPOSITION .OR. OUTPUT_WALL_QUANTITIES)) THEN CALL WALL_MODEL(SLIP_COEF,B2%U_TAU,B2%Y_PLUS,MU_DNS(BC%IIG,BC%JJG,BC%KKG)/RHO(BC%IIG,BC%JJG,BC%KKG),SF%ROUGHNESS,& 0.5_EB/B1%RDN,B1%U_TANG) ENDIF