Broadcast properly cell_measures when using extract_shape with decomposed: True

[sloosvel]

Description

This PR makes sure that time-dependent cell measures are broadcasted properly to a cube with a shape_id coordinate.
Closes #2347

Link to documentation:

---

Before you get started

• ☝ Create an issue to discuss what you are going to do

Checklist

It is the responsibility of the author to make sure the pull request is ready to review. The icons indicate whether the item will be subject to the 🛠 Technical or 🧪 Scientific review.

• 🧪 The new functionality is relevant and scientifically sound
• 🛠 This pull request has a descriptive title and labels
• 🛠 Code is written according to the code quality guidelines
• 🧪 and 🛠 Documentation is available
• 🛠 Unit tests have been added
• 🛠 Changes are backward compatible
• 🛠 Any changed dependencies have been added or removed correctly
• 🛠 The list of authors is up to date
• 🛠 All checks below this pull request were successful

---

To help with the number pull requests:

• 🙏 We kindly ask you to review two other open pull requests in this repository

[codecov]

Codecov Report

All modified and coverable lines are covered by tests ✅

Project coverage is 94.10%. Comparing base (3a0b49f) to head (66c7ee0).

Additional details and impacted files

@@           Coverage Diff           @@
##             main    #2348   +/-   ##
=======================================
  Coverage   94.10%   94.10%           
=======================================
  Files         246      246           
  Lines       13477    13482    +5     
=======================================
+ Hits        12682    12687    +5     
  Misses        795      795           

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

[schlunma]

Thanks Saskia, looks good to me, just have two minor comments.

[sloosvel]

Thanks for the comments @schlunma , I added a comment and changed a bit the condition to also make it work when only one region is extracted and the result cube has no shape_id dimension. Here is a recipe to test both cases:

  description: |
    This a recipe to monitor experiments by automatizing the creation of plots
    from raw preprocessor output.
  authors:
    - vegas-regidor_javier
  maintainer:
    - loosveldt-tomas_saskia


datasets:
- dataset: EC-Earth3-CC
  ensemble: r1i1p1f1
  grid: gn
  exp: historical
  mip: Omon
  project: CMIP6

preprocessors:
  extract_two_regions:
    extract_shape:
      shapefile: AR6
      method: contains
      crop: False
      decomposed: True
      ids: {Acronym: [EAO, WAF]}
    volume_statistics:
      operator: mean
  extract_one_region:
    extract_shape:
      shapefile: AR6
      method: contains
      crop: False
      decomposed: True
      ids: {Acronym: [EAO]}
    volume_statistics:
      operator: mean

diagnostics:
  timeseries:
    variables:
      thetao_1:
        short_name: thetao
        timerange: '1950/1952'
        preprocessor: extract_two_regions
      thetao_2:
        short_name: thetao
        timerange: '1950/1952'
        preprocessor: extract_one_region
    scripts: null

[schlunma]

Thanks Saskia for making the changes, everything looks good now! The output of the recipe looks as expected 🚀