add the function for calculating pore throat radius