diff --git a/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex b/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex
index 1fd02d2f204..657fe352a87 100644
--- a/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex
+++ b/Manuals/FDS_Verification_Guide/FDS_Verification_Guide.tex
@@ -2252,7 +2252,7 @@ \section{Monin-Obukhov Similarity Profiles (\texorpdfstring{\ct{MO_velocity_prof
 \label{MO_velocity_profile_stable}
 \label{MO_velocity_profile_unstable}
 
-Atmospheric turbulence is affected by the stability of the boundary layer.  A stable layer (cool, heavy air at ground level) will suppress turbulence, while an unstable layer (warm, light air at ground level) will enhance turbulent mixing as buoyant plumes rise.  The theory governing these flows, Monin-Obukhov similarity theory, is discussed at some length in the FDS User Guide \cite{FDS_Users_Guide}.  Here we examine velocity profiles from a stable boundary layer and an unstable boundary layer.  When mean forcing is used for driving the wind field the Monin-Obukhov parameters determine the shape of the mean streamwise velocity profile.  These cases use a very tight \ct{DT_MEAN_FORCING_2} of 0.1 s in order to drive the flow field directly to the specified profile, therefore comfirming the target profile is being computed correctly in FDS.
+Atmospheric turbulence is affected by the stability of the boundary layer.  A stable layer (cool, heavy air at ground level) will suppress turbulence, while an unstable layer (warm, light air at ground level) will enhance turbulent mixing as buoyant plumes rise.  The theory governing these flows, Monin-Obukhov similarity theory, is discussed at some length in the FDS User Guide \cite{FDS_Users_Guide}.  Here we examine velocity profiles from a stable boundary layer and an unstable boundary layer. The Monin-Obukhov parameters determine the shape of the mean streamwise velocity profile.
 
 \begin{figure}[ht]
    \begin{tabular*}{\textwidth}{l@{\extracolsep{\fill}}r}
@@ -5941,11 +5941,11 @@ \subsection{Heating a Metal Sphere via Radiation and Convection}
 \label{particle_heating_convection}
 \label{particle_heating_radiation}
 
-A small metal sphere with mass, $m_{\rm s}=0.005$~kg, is suspended in a 1~m cube filled with $m_{\rm g}=0.318$~kg nitrogen with a specified specific heat, $c_p=1$~kJ/(kg$\cdot$K), and initial temperature, $T_{\rm g,i}=1073.15$~K. The metal has a specified specific heat, $c_{\rm s}=1$~kJ/(kg$\cdot$K), and initial temperature, $T_{\rm s,i}=293.15$~K. The walls of the box are adiabatic. In the first case, \ct{particle\_heating\_convection}, the sphere is heated via convection only and there is no radiation heat transfer. In the second case, \ct{particle\_heating\_radiation}, the sphere is heated via radiation only. The initial heat flux in both cases is approximately 75~kW/m$^2$, and the final temperature in both cases, $T_{\rm f}$, is found from solving an equation that equates the internal energy gained by the solid with the internal energy lost by the gas:
+Ten small metal spheres, each with mass, $m_{\rm s}=0.005$~kg, are suspended in a 1~m cube filled with $m_{\rm g}=0.318$~kg nitrogen with a specified specific heat, $c_p=1$~kJ/(kg$\cdot$K), and initial temperature, $T_{\rm g,i}=1073.15$~K. The metal has a specified specific heat, $c_{\rm s}=1$~kJ/(kg$\cdot$K), and initial temperature, $T_{\rm s,i}=293.15$~K. The walls of the box are adiabatic. In the first case, \ct{particle\_heating\_convection}, the spheres are heated via convection only and there is no radiation heat transfer. In the second case, \ct{particle\_heating\_radiation}, the spheres are heated via radiation only. The initial heat flux in both cases is approximately 75~kW/m$^2$, and the final temperature in both cases, $T_{\rm f}$, is found from solving an equation that equates the internal energy gained by the solid with the internal energy lost by the gas:
 \be
-   m_{\rm s} \, c_{\rm s} \, (T_{\rm f}-T_{\rm s,i}) = m_{\rm g} \, c_{v} \, (T_{\rm g,i}-T_{\rm f})  \quad ; \quad  c_v = c_p - \frac{R}{W} = 1 - \frac{8.3145}{28} \approx 0.703 \; \hbox{kJ/(kg}\cdot\hbox{K})
+   10 \, m_{\rm s} \, c_{\rm s} \, (T_{\rm f}-T_{\rm s,i}) = m_{\rm g} \, c_{v} \, (T_{\rm g,i}-T_{\rm f})  \quad ; \quad  c_v = c_p - \frac{R}{W} = 1 - \frac{8.3145}{28} \approx 0.703 \; \hbox{kJ/(kg}\cdot\hbox{K})
 \ee
-The final temperature, $T_{\rm f}=1056.1$~K or 782.9~$^\circ$C, as shown in Fig.~\ref{particle_heating_figs}.
+The final temperature, $T_{\rm f}=931.4$~K or 658.3~$^\circ$C, as shown in Fig.~\ref{particle_heating_figs}.
 \begin{figure}[ht]
 \includegraphics[height=2.2in]{SCRIPT_FIGURES/particle_heating_convection}
 \includegraphics[height=2.2in]{SCRIPT_FIGURES/particle_heating_radiation}
diff --git a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv
index d4399cf7a8d..77a8604deb1 100644
--- a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv
+++ b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv
@@ -486,10 +486,10 @@ d,particle_drag_U100_N1600,Sprinklers_and_Sprays/particle_drag_U100_N1600_git.tx
 d,particle_drag_U150_N1600,Sprinklers_and_Sprays/particle_drag_U150_N1600_git.txt,Sprinklers_and_Sprays/particle_drag_U150_N1600.csv,1,2,T,U,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_drag_U150_N1600_devc.csv,2,3,Time,U-VEL,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Gas phase velocity (particle\_drag\_F),Time (s),Velocity (m/s),0,2,1,0,200,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_drag_F,Absolute Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
 d,particle_drag_U10_N16,Sprinklers_and_Sprays/particle_drag_U10_N16_git.txt,Sprinklers_and_Sprays/particle_drag_U10_N16.csv,1,2,T,F,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_drag_U10_N16_devc.csv,2,3,Time,drag force,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Total drag force (particle\_drag\_A),Time (s),Force (N),0,100,1,0,1,-1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_drag_sum_A,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
 d,particle_flux,Sprinklers_and_Sprays/particle_flux_git.txt,Sprinklers_and_Sprays/particle_flux.csv,1,2,Time,mass,Expected (mass),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_flux_devc.csv,2,3,Time,mass,FDS (mass),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Particle Mass (particle\_flux),Time (s),Mass (kg),0,20,1,0,0.8,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/particle_flux,Relative Error,end,0.01,Sprinklers and Sprays,bs,b,TeX
-d,particle_heating_convection,Sprinklers_and_Sprays/particle_heating_convection_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_convection_devc.csv,2,3,Time,T_gas,FDS (T\_gas),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_convection),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_convection,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
-f,particle_heating_convection,Sprinklers_and_Sprays/particle_heating_convection_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_convection_devc.csv,2,3,Time,T_ball,FDS (T\_ball),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_convection),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_convection,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
-d,particle_heating_radiation,Sprinklers_and_Sprays/particle_heating_radiation_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_radiation_devc.csv,2,3,Time,T_gas,FDS (T\_gas),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_radiation),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_radiation,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
-f,particle_heating_radiation,Sprinklers_and_Sprays/particle_heating_radiation_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_radiation_devc.csv,2,3,Time,T_ball,FDS (T\_ball),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_radiation),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_radiation,Relative Error,end,0.01,Sprinklers and Sprays,kd,k,TeX
+d,particle_heating_convection,Sprinklers_and_Sprays/particle_heating_convection_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_convection_devc.csv,2,3,Time,T_gas,FDS (T\_gas),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_convection),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_convection,Relative Error,end,0.02,Sprinklers and Sprays,kd,k,TeX
+f,particle_heating_convection,Sprinklers_and_Sprays/particle_heating_convection_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_convection_devc.csv,2,3,Time,T_ball,FDS (T\_ball),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_convection),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_convection,Relative Error,end,0.02,Sprinklers and Sprays,kd,k,TeX
+d,particle_heating_radiation,Sprinklers_and_Sprays/particle_heating_radiation_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,Analytical,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_radiation_devc.csv,2,3,Time,T_gas,FDS (T\_gas),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_radiation),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_radiation,Relative Error,end,0.02,Sprinklers and Sprays,kd,k,TeX
+f,particle_heating_radiation,Sprinklers_and_Sprays/particle_heating_radiation_git.txt,Sprinklers_and_Sprays/particle_heating.csv,1,2,Time,Temp,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_heating_radiation_devc.csv,2,3,Time,T_ball,FDS (T\_ball),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (particle\_heating\_radiation),Time (s),Temperature (°C),0,900,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_heating_radiation,Relative Error,end,0.02,Sprinklers and Sprays,kd,k,TeX
 d,particle_isotropic_radi,Sprinklers_and_Sprays/particle_isotropic_radi_git.txt,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Delta e_ball,E\_w ball,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Q_rad ball,Q\_r ball,k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Enthalpy (particle\_isotropic\_radiation),Time (s),Enthalpy (kJ),0,50,1,0,12,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_isotropic_radi,Relative Error,end,0.015,Sprinklers and Sprays,ko,k,TeX
 f,particle_isotropic_radi,Sprinklers_and_Sprays/particle_isotropic_radi_git.txt,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Delta e_cyl,E\_w cyl,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Q_rad cyl,Q\_r cyl,r--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Enthalpy (particle\_isotropic\_radiation),Time (s),Enthalpy (kJ),0,50,1,0,12,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_isotropic_radi,Relative Error,end,0.015,Sprinklers and Sprays,ko,k,TeX
 f,particle_isotropic_radi,Sprinklers_and_Sprays/particle_isotropic_radi_git.txt,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Delta e_plate,E\_w plate,b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/particle_isotropic_radi_devc.csv,2,3,Time,Q_rad plate,Q\_r plate,b--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Enthalpy (particle\_isotropic\_radiation),Time (s),Enthalpy (kJ),0,50,1,0,12,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/particle_isotropic_radi,Relative Error,end,0.015,Sprinklers and Sprays,ko,k,TeX
diff --git a/Utilities/Matlab/scripts/atmospheric_boundary_layer.m b/Utilities/Matlab/scripts/atmospheric_boundary_layer.m
index 79e43a7c5d1..274eedcd56f 100644
--- a/Utilities/Matlab/scripts/atmospheric_boundary_layer.m
+++ b/Utilities/Matlab/scripts/atmospheric_boundary_layer.m
@@ -67,6 +67,11 @@
 
 T = T + (theta_0-T(12));
 
+ERROR = abs(u(end)-M2.data(end,2));
+if ERROR>2.
+   display(['Matlab Warning: atmospheric_boundary_layer Case ',num2str(i),' velocity out of tolerance. ERROR = ',num2str(ERROR),' m/s'])
+end
+
 plot(u,z,'ko'); hold on
 plot(M2.data(:,2),M2.data(:,1),'k-'); hold off
 set(gca,'FontName',Font_Name)
@@ -99,6 +104,11 @@
 set(gca,'Units',Plot_Units)
 set(gca,'Position',[Plot_X Plot_Y Plot_Width Plot_Height])
 
+ERROR = abs(T(end)-273.15-M2.data(end,3));
+if ERROR>1.0
+   display(['Matlab Warning: atmospheric_boundary_layer Case ',num2str(i),' temperature out of tolerance. ERROR = ',num2str(ERROR),' K'])
+end
+
 plot(T-273.15,z,'ko'); hold on
 plot(M2.data(:,3),M2.data(:,1),'k-'); hold off
 set(gca,'FontName',Font_Name)
diff --git a/Verification/Atmospheric_Effects/MO_velocity_profile_stable.fds b/Verification/Atmospheric_Effects/MO_velocity_profile_stable.fds
index 894828f85b1..e0d2a3227d8 100644
--- a/Verification/Atmospheric_Effects/MO_velocity_profile_stable.fds
+++ b/Verification/Atmospheric_Effects/MO_velocity_profile_stable.fds
@@ -7,7 +7,7 @@
 
 &TIME T_END=10. /
 
-&MISC P_INF=90630, HUMIDITY=12.0, CHECK_HT=T /
+&MISC P_INF=90630, HUMIDITY=12.0 /
 &WIND L=69.4, U_STAR=0.369, Z_0=0.008, THETA_STAR=0.152, TMP_REF=30.8, Z_REF=1. /
 &RADI RADIATION=.FALSE. /
 
diff --git a/Verification/Complex_Geometry/geom_channel2.fds b/Verification/Complex_Geometry/geom_channel2.fds
index 9176ef8c635..1bcc8061a4a 100644
--- a/Verification/Complex_Geometry/geom_channel2.fds
+++ b/Verification/Complex_Geometry/geom_channel2.fds
@@ -15,7 +15,7 @@
 
 &RADI RADIATION=F /
 &PRES VELOCITY_TOLERANCE=1.E-6 /
-&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0., CHECK_HT=T /
+&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0. /
 &SURF ID='INLET', MASS_FLUX(1)=1., SPEC_ID(1)='TRACER', TAU_MF(1)=0./
 
 &VENT MB='XMIN', SURF_ID='INLET' /
diff --git a/Verification/Complex_Geometry/geom_channel_tmp.fds b/Verification/Complex_Geometry/geom_channel_tmp.fds
index 40a3b86fe70..0fcfe94901c 100644
--- a/Verification/Complex_Geometry/geom_channel_tmp.fds
+++ b/Verification/Complex_Geometry/geom_channel_tmp.fds
@@ -18,7 +18,7 @@
 
 &RADI RADIATION=F /
 &PRES SOLVER='ULMAT'/
-&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0., CHECK_HT=T /
+&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0. /
 &SURF ID='INLET', FREE_SLIP=T, MASS_FRACTION(1)=1., SPEC_ID(1)='TRACER', VEL=-1, TAU_V=0., TMP_FRONT=600., RAMP_T='t1'/
 
 &RAMP ID='t1', T=0.,F=1./
diff --git a/Verification/Complex_Geometry/geom_channel_tmp2.fds b/Verification/Complex_Geometry/geom_channel_tmp2.fds
index 0401178c03b..b06230f2996 100644
--- a/Verification/Complex_Geometry/geom_channel_tmp2.fds
+++ b/Verification/Complex_Geometry/geom_channel_tmp2.fds
@@ -14,7 +14,7 @@
 
 &RADI RADIATION=F /
 &PRES VELOCITY_TOLERANCE=1.E-6 /
-&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0., CHECK_HT=T /
+&MISC NOISE=F, STRATIFICATION=F, GVEC(3)=0. /
 &SURF ID='INLET', TMP_FRONT=600., RAMP_T='t1'/
 
 &RAMP ID='t1', T=0.,F=1./
diff --git a/Verification/Complex_Geometry/geom_obst.fds b/Verification/Complex_Geometry/geom_obst.fds
index 3c4fcacdd66..1e9ef12d69d 100644
--- a/Verification/Complex_Geometry/geom_obst.fds
+++ b/Verification/Complex_Geometry/geom_obst.fds
@@ -3,7 +3,6 @@
 &TIME T_END=1.0/
 &MESH IJK=15,15,15,XB=-1.0,1.0,-1.0,1.0,-1.0,1.0 /
 
-&MISC CHECK_HT=T /
 &REAC SOOT_YIELD=0.01,FUEL='PROPANE'/
 &SURF ID='BURNER',HRRPUA=600.0 / 
 
diff --git a/Verification/Heat_Transfer/back_wall_test_2.fds b/Verification/Heat_Transfer/back_wall_test_2.fds
index f93032ee1e9..e7eeecc3a75 100644
--- a/Verification/Heat_Transfer/back_wall_test_2.fds
+++ b/Verification/Heat_Transfer/back_wall_test_2.fds
@@ -1,11 +1,10 @@
-&HEAD CHID='back_wall_test_2', TITLE='Test 1-D heat transfer through rotated GEOM obstruction' /
+&HEAD CHID='back_wall_test_2_noht', TITLE='Test 1-D heat transfer through rotated GEOM obstruction' /
 
 &MESH XB=0,1,0,1,0,1, IJK=40,40,40 /
 &MESH XB=2,3,0,1,0,1, IJK=40,40,40 /
 &MESH XB=4,5,0,1,0,1, IJK=40,40,40 /
 
 &TIME T_END=20., DT=0.1 / 
-&MISC CFL_MAX=0.8, CHECK_HT=T / 
 &SURF ID='HOT', COLOR='RED', TMP_FRONT=1000., TAU_T=0., EMISSIVITY=1, HEAT_TRANSFER_COEFFICIENT=0 /
 &SURF ID='COLD',COLOR='BLUE',TMP_FRONT=20.,   TAU_T=0., EMISSIVITY=1, HEAT_TRANSFER_COEFFICIENT=0 /
 
diff --git a/Verification/Pressure_Solver/ulmat_2zone.fds b/Verification/Pressure_Solver/ulmat_2zone.fds
index 4d38c294544..06cc40a77ca 100644
--- a/Verification/Pressure_Solver/ulmat_2zone.fds
+++ b/Verification/Pressure_Solver/ulmat_2zone.fds
@@ -3,20 +3,11 @@
 &MESH IJK=10,10,5, XB=0,1,0,1,0,0.5/
 &MESH IJK=10,10,5, XB=0,1,0,1,0.5,1/
 
-#&MESH IJK=5,10,10, XB=0,0.5,0,1,0,1./
-#&MESH IJK=5,10,10, XB=0.5,1,0,1,0,1./
-
-#&MESH IJK=5,10,5, XB=0,0.5,0,1,0,0.5/
-#&MESH IJK=5,10,5, XB=0,0.5,0,1,0.5,1/
-#&MESH IJK=5,10,5, XB=0.5,1,0,1,0,0.5/
-#&MESH IJK=5,10,5, XB=0.5,1,0,1,0.5,1/
-
-
 &TIME T_END=10/
 
 &PRES CHECK_POISSON=T, SOLVER='ULMAT'/
 
-&MISC NOISE=F, CHECK_HT=T /
+&MISC NOISE=F /
 
 &RADI RADIATION=F/
 
diff --git a/Verification/Sprinklers_and_Sprays/particle_heating.csv b/Verification/Sprinklers_and_Sprays/particle_heating.csv
index e4617b5b708..e31317e3104 100644
--- a/Verification/Sprinklers_and_Sprays/particle_heating.csv
+++ b/Verification/Sprinklers_and_Sprays/particle_heating.csv
@@ -1,4 +1,4 @@
 Time,Temp
-700,782.9
-800,782.9
-900,782.9
+700,658.26
+800,658.26
+900,658.26
diff --git a/Verification/Sprinklers_and_Sprays/particle_heating_convection.fds b/Verification/Sprinklers_and_Sprays/particle_heating_convection.fds
index 44f5dcefecd..3d91a1c25e0 100644
--- a/Verification/Sprinklers_and_Sprays/particle_heating_convection.fds
+++ b/Verification/Sprinklers_and_Sprays/particle_heating_convection.fds
@@ -1,8 +1,8 @@
-&HEAD CHID='particle_heating_convection', TITLE='Test of particle energy conservation' /
+&HEAD CHID='particle_heating_convection_noht', TITLE='Test of particle energy conservation' /
 
-&TIME T_END=900, DT=0.1, WALL_INCREMENT=1 /
+&TIME T_END=900, DT=0.02, WALL_INCREMENT=1 /
 
-&MISC TMPA=800., CHECK_HT=T /
+&MISC TMPA=800. /
 
 &MESH IJK=10,10,10 XB=0.0,1.0,0.0,1.0,0.0,1.0 /
 
@@ -25,8 +25,8 @@
 
 &INIT ID                   = 'ball',
       PART_ID              = 'ball'
-      XYZ                  = 0.45,0.45,0.45
-      N_PARTICLES          = 1 /
+      XB                   = 0.0,1.0,0.0,1.0,0.0,1.0
+      N_PARTICLES          = 10 /
 
 &MATL ID                 = 'metal'
       DENSITY            = 5000.
diff --git a/Verification/Sprinklers_and_Sprays/particle_heating_radiation.fds b/Verification/Sprinklers_and_Sprays/particle_heating_radiation.fds
index 9cc0cd251e6..637df388dea 100644
--- a/Verification/Sprinklers_and_Sprays/particle_heating_radiation.fds
+++ b/Verification/Sprinklers_and_Sprays/particle_heating_radiation.fds
@@ -1,8 +1,8 @@
-&HEAD CHID='particle_heating_radiation', TITLE='Test of particle energy conservation' /
+&HEAD CHID='particle_heating_radiation_noht', TITLE='Test of particle energy conservation' /
 
-&TIME T_END=900, DT=0.1, WALL_INCREMENT=1 /
+&TIME T_END=900, DT=0.02, WALL_INCREMENT=1 /
 
-&MISC TMPA=800., CHECK_HT=T /
+&MISC TMPA=800. /
 
 &MESH IJK=10,10,10 XB=0.0,1.0,0.0,1.0,0.0,1.0 /
 
@@ -25,8 +25,8 @@
 
 &INIT ID                   = 'ball',
       PART_ID              = 'ball'
-      XYZ                  = 0.45,0.45,0.45
-      N_PARTICLES          = 1 /
+      XB                   = 0.0,1.0,0.0,1.0,0.0,1.0
+      N_PARTICLES          = 10 /
 
 &MATL ID                 = 'metal'
       DENSITY            = 5000.
@@ -42,4 +42,6 @@
 &DEVC XB=0,1,0,1,0,1, QUANTITY='CONVECTIVE HEAT FLUX', ID='wall con flux', SURF_ID='WALL', SPATIAL_STATISTIC='SURFACE INTEGRAL' /
 &DEVC XB=0,1,0,1,0,1, QUANTITY='RADIATIVE HEAT FLUX', ID='wall rad flux', SURF_ID='WALL', SPATIAL_STATISTIC='SURFACE INTEGRAL' /
 
+&BNDF QUANTITY='CONVECTIVE HEAT FLUX', CELL_CENTERED=T /
+
 &TAIL /