diff --git a/sources/interpolations.F90 b/sources/interpolations.F90
index 7036f43..1cf144c 100644
--- a/sources/interpolations.F90
+++ b/sources/interpolations.F90
@@ -4501,8 +4501,8 @@ module interpolations
     real(kind=8), dimension(size(fc)) :: r
     real(kind=8), dimension(size(fc)) :: u
 
-    real(kind=8), parameter :: a1 = 5.0163016d-01
-    real(kind=8), parameter :: a2 = 2.5394716d-01
+    real(kind=8), parameter :: a1 = 5.0160152430504783780454277648572d-01
+    real(kind=8), parameter :: a2 = 2.5395800542237252167985944322018d-01
     real(kind=8), dimension(3), parameter ::                                   &
             di5 = [ a1, 1.0d+00, a2 ]
     real(kind=8), dimension(5), parameter ::                                   &
@@ -4628,8 +4628,8 @@ module interpolations
     real(kind=8), dimension(size(fc)) :: r
     real(kind=8), dimension(size(fc)) :: u
 
-    real(kind=8), parameter :: a1 = 6.6850691831375709863684029643567d-01
-    real(kind=8), parameter :: a2 = 3.3644225201902153852210572056440d-01
+    real(kind=8), parameter :: a1 = 6.7080731007136045597206649140123d-01
+    real(kind=8), parameter :: a2 = 3.4031334443871704940876966367699d-01
     real(kind=8), dimension(5), parameter ::                                   &
             di7 = [ (3.0d+00 * a1 + 2.0d+00 * a2 - 2.0d+00) / 8.0d+00,         &
                     a1, 1.0d+00, a2,                                           &
@@ -4769,8 +4769,8 @@ module interpolations
     real(kind=8), dimension(size(fc)) :: r
     real(kind=8), dimension(size(fc)) :: u
 
-    real(kind=8), parameter :: a1 = 6.6775815668016370843881215557557d-01
-    real(kind=8), parameter :: a2 = 4.0271201156292879402181110913917d-01
+    real(kind=8), parameter :: a1 = 6.7049430606778809680492458975703d-01
+    real(kind=8), parameter :: a2 = 4.0910274527161221589360057231066d-01
     real(kind=8), dimension(5), parameter ::                                   &
             di9 = [ (9.0d+00 * a1 + 5.0d+00 * a2 - 6.0d+00) / 2.0d+01,         &
                     a1, 1.0d+00, a2,                                           &