Merge pull request #13341 from mcgratta/master

FDS Source: Issue #13340. Remove unused variables

Source/mesh.f90

@@ -54,7 +54,6 @@ MODULE MESH_VARIABLES
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: MIX_TIME    !< Mixing-controlled combustion reaction time (s)
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: STRAIN_RATE !< Strain rate \f$ |S|_{ijk} \f$ (1/s)
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: D_Z_MAX     !< \f$ \max D_\alpha \f$
-   REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: Q_DOT_PPP_S !< Heat release rate per unit volume in 3D pyrolysis model
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: PP_RESIDUAL !< Pressure Poisson residual (debug)
    REAL(EB), ALLOCATABLE, DIMENSION(:,:,:) :: LES_FILTER_WIDTH !< Characteristic cell dimension (m)
 
@@ -93,7 +92,6 @@ MODULE MESH_VARIABLES
    REAL(EB), ALLOCATABLE, DIMENSION(:) :: D_PBAR_DT   !< \f$ (\partial \overline{p}_m/\partial t)^n \f$
    REAL(EB), ALLOCATABLE, DIMENSION(:) :: D_PBAR_DT_S !< \f$ (\partial \overline{p}_m/\partial t)^* \f$
    REAL(EB), ALLOCATABLE, DIMENSION(:) :: U_LEAK
-   REAL(EB), ALLOCATABLE, DIMENSION(:) :: U_DUCT
    REAL(EB), ALLOCATABLE, DIMENSION(:) :: U_WIND      !< Component of wind in x direction
    REAL(EB), ALLOCATABLE, DIMENSION(:) :: V_WIND      !< Component of wind in y direction
    REAL(EB), ALLOCATABLE, DIMENSION(:) :: W_WIND      !< Component of wind in z direction
@@ -343,7 +341,7 @@ MODULE MESH_POINTERS
 REAL(EB), POINTER, DIMENSION(:,:,:) :: &
    U,V,W,US,VS,WS,DDDT,D,DS,H,HS,KRES,FVX,FVY,FVZ,FVX_B,FVY_B,FVZ_B,FVX_D,FVY_D,FVZ_D,RHO,RHOS, &
    MU,MU_DNS,TMP,Q,KAPPA_GAS,CHI_R,QR,QR_W,UII,RSUM,D_SOURCE, &
-   CSD2,MTR,MSR,WEM,MIX_TIME,CHEM_SUBIT,STRAIN_RATE,D_Z_MAX,Q_DOT_PPP_S,PP_RESIDUAL,LES_FILTER_WIDTH
+   CSD2,MTR,MSR,WEM,MIX_TIME,CHEM_SUBIT,STRAIN_RATE,D_Z_MAX,PP_RESIDUAL,LES_FILTER_WIDTH
 REAL(EB), POINTER, DIMENSION(:,:,:,:) :: ZZ,ZZS,REAC_SOURCE_TERM,DEL_RHO_D_DEL_Z,FX,FY,FZ, &
                                          SCALAR_WORK1,SCALAR_WORK2,SCALAR_WORK3,SCALAR_WORK4, &
                                          Q_REAC,AVG_DROP_DEN,AVG_DROP_TMP,AVG_DROP_RAD,AVG_DROP_AREA,M_DOT_PPP, &
@@ -353,7 +351,7 @@ MODULE MESH_POINTERS
 REAL(EB), POINTER, DIMENSION(:,:,:) :: PRHS
 REAL(EB), POINTER, DIMENSION(:,:) :: BXS,BXF,BYS,BYF,BZS,BZF, BXST,BXFT,BYST,BYFT,BZST,BZFT
 INTEGER, POINTER :: LSAVE,LWORK,LBC,MBC,NBC,LBC2,MBC2,NBC2,ITRN,JTRN,KTRN,IPS
-REAL(EB), POINTER, DIMENSION(:) :: P_0,RHO_0,TMP_0,D_PBAR_DT,D_PBAR_DT_S,U_LEAK,U_DUCT,U_WIND,V_WIND,W_WIND
+REAL(EB), POINTER, DIMENSION(:) :: P_0,RHO_0,TMP_0,D_PBAR_DT,D_PBAR_DT_S,U_LEAK,U_WIND,V_WIND,W_WIND
 REAL(EB), POINTER, DIMENSION(:,:) :: PBAR,PBAR_S,R_PBAR
 INTEGER, POINTER, DIMENSION(:,:,:) :: PRESSURE_ZONE
 TYPE (ZONE_MESH_TYPE), POINTER, DIMENSION(:) :: ZONE_MESH
@@ -508,7 +506,6 @@ SUBROUTINE POINT_TO_MESH(NM)
 MIX_TIME=>M%MIX_TIME
 Q=>M%Q
 Q_REAC=>M%Q_REAC
-Q_DOT_PPP_S=>M%Q_DOT_PPP_S
 CHI_R => M%CHI_R
 QR=>M%QR
 QR_W=>M%QR_W
@@ -574,7 +571,6 @@ SUBROUTINE POINT_TO_MESH(NM)
 KTRN=>M%KTRN
 IPS=>M%IPS
 U_LEAK=>M%U_LEAK
-U_DUCT=>M%U_DUCT
 U_WIND=>M%U_WIND
 V_WIND=>M%V_WIND
 W_WIND=>M%W_WIND