Added utility routines cbrt to compute cubic roots (used by some planetary physics packages) and getin_p (to read from def files in parallel). EM