Add the function for calculating the radius of pore throats