1st index always refers to the region so mpk_P(1,:) is the Power spectrum in the first active region
this will be called once for each redshift bin