diff --git a/src/fix_nh.cpp b/src/fix_nh.cpp
index b254636111..a78a109396 100644
--- a/src/fix_nh.cpp
+++ b/src/fix_nh.cpp
@@ -835,9 +835,17 @@ void FixNH::final_integrate()
 {
   nve_v();
 
+  // re-compute temp before nh_v_press()
+  // only needed for temperature computes with BIAS:
+  //   b/c some biases store per-atom values (e.g. temp/profile) 
+  //   per-atom values are invalid if reneigh/comm occurred
+  //     since temp->compute() in initial_integrate()
+
+  if (which == BIAS) t_current = temperature->compute_scalar();
+
   if (pstat_flag) nh_v_press();
 
-  // compute new T,P
+  // compute new T,P after velocities rescaled by nh_v_press()
   // compute appropriately coupled elements of mvv_current
 
   t_current = temperature->compute_scalar();