real->double precision to match caller

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