Needs to: * Allocate planar recharge rates or constants to uppermost active cells. Blocked by: https://github.com/Deltares/imod-python/issues/949 We think boundary conditions always need to be assigned to a active cell, but this can be verified.