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