The computation of non-Darcy properties is expensive, so we cache the results. Make sure that the calling function calls this function before using the properties.