add option to set the NHI where Lyman limit absorption is included. Default value is log10(NHI) > 19