interpolation code for 3d data !!



Hi,
   I am looking for a code(in fortran) for interpolation of a large set
 of 3 dimensional data (on an eqaully spaced grid). I tried the polynomial
 routine in numerical recipe -but it is very poor for my data. Probably
 spline or hermite polynomial based interpolation scheme would be the best.
 I would appreciate any help.
 thanks.
          Pradipta