Currently calculate_ps_lc supports having chunks of different sizes (in configuration space) with the chunk_size argument. It could be convenient if the chunks could automatically be set to match the frequency bands of a given experiment. I imagine having something like an argument called freq_bands, and given this input the chunks boundaries as well as the redshift that is associated with each chunk could be determined.