one parameter epsilon added to take into account the rugosity of the grains, the parameter is now set to zero