diff --git a/RESAMPLE/fpdisc.f b/RESAMPLE/fpdisc.f
index 51ad06a74f95ab795279df4a4d298696b419955d..5e1d047e3f92e14a002d01fc3385826aa6f352f8 100644
--- a/RESAMPLE/fpdisc.f
+++ b/RESAMPLE/fpdisc.f
@@ -4,12 +4,12 @@
 !c  ..scalar arguments..
       integer n,k2,nest
 !c  ..array arguments..
-      real t(n),b(nest,k2)
+      double precision t(n),b(nest,k2)
 !c  ..local scalars..
-      real an,fac,prod
+      double precision an,fac,prod
       integer i,ik,j,jk,k,k1,l,lj,lk,lmk,lp,nk1,nrint
 !c  ..local array..
-      real h(12)
+      double precision h(12)
 !c  ..
       k1 = k2-1
       k = k1-1