diff --git a/ProcessLib/HT/HTFEM.h b/ProcessLib/HT/HTFEM.h
index cbb55c4a4952bef43ce09d13201123fc5165fd2e..a237c8d51a025283d5ec63789de44132c0a1e035 100644
--- a/ProcessLib/HT/HTFEM.h
+++ b/ProcessLib/HT/HTFEM.h
@@ -306,6 +306,8 @@ public:
             auto const& N = ip_data.N;
             auto const& dNdx = ip_data.dNdx;
 
+            pos.setIntegrationPoint(ip);
+
             double T_int_pt = 0.0;
             double p_int_pt = 0.0;
             NumLib::shapeFunctionInterpolate(local_x, N, T_int_pt, p_int_pt);