Calculates the value of rho for every surface in the given rho file. Used later in rho calculation, to interpolate between surfaces
dommaschk_t