cosmo fx hdr

function f_ds = cosmo_fx(ds, f, split_by, dim, check)