F_colloidal is kept for h < 0 as the value at h = 0.
Many other modifications.