pmax = np.minimum(pmax, -_TINY)  # Make sure -15 < logsf < 0