/*tababun interpolate on table of points to do 'element table' command, */
double tababun(double r0, 
  double depth, 
  long int iel);