/*tabden interpolate on table of points for density with dlaw table command, by K Volk */
double tabden(double r0, 
  double depth);