Add trilinear interpolation (interp3d.f90)