diff --git a/src/phase_mechanical_eigen_thermalexpansion.f90 b/src/phase_mechanical_eigen_thermalexpansion.f90
index d03ea0eb2..3c7654dc9 100644
--- a/src/phase_mechanical_eigen_thermalexpansion.f90
+++ b/src/phase_mechanical_eigen_thermalexpansion.f90
@@ -98,14 +98,14 @@ module subroutine thermalexpansion_LiAndItsTangent(Li, dLi_dTstar, ph,me)
 
   associate(prm => param(kinematics_thermal_expansion_instance(ph)))
     Li = dot_T * ( &
-                  prm%A(1:3,1:3,1)*1.0_pReal &                                                      ! constant  coefficient
-                + prm%A(1:3,1:3,2)*(T - prm%T_ref)**1 &                                             ! linear    coefficient
-                + prm%A(1:3,1:3,3)*(T - prm%T_ref)**2 &                                             ! quadratic coefficient
+                  prm%A(1:3,1:3,1) &                                                                ! constant  coefficient
+                + prm%A(1:3,1:3,2)*(T - prm%T_ref)**1.0_pReal &                                     ! linear    coefficient
+                + prm%A(1:3,1:3,3)*(T - prm%T_ref)**2.0_pReal &                                     ! quadratic coefficient
                 ) / &
          (1.0_pReal &
-               + prm%A(1:3,1:3,1)*(T - prm%T_ref)**1 / 1. &
-               + prm%A(1:3,1:3,2)*(T - prm%T_ref)**2 / 2. &
-               + prm%A(1:3,1:3,3)*(T - prm%T_ref)**3 / 3. &
+               + prm%A(1:3,1:3,1)*(T - prm%T_ref)**1.0_pReal / 1.0_pReal &
+               + prm%A(1:3,1:3,2)*(T - prm%T_ref)**2.0_pReal / 2.0_pReal &
+               + prm%A(1:3,1:3,3)*(T - prm%T_ref)**3.0_pReal / 3.0_pReal &
          )
   end associate
   dLi_dTstar = 0.0_pReal