Skip to content

Fix: read baselines from grib instead of zarr for meteograms#164

Merged
jonasbhend merged 4 commits into
unify_regions_stationsfrom
fix_for_meteogram
May 28, 2026
Merged

Fix: read baselines from grib instead of zarr for meteograms#164
jonasbhend merged 4 commits into
unify_regions_stationsfrom
fix_for_meteogram

Conversation

@jonasbhend
Copy link
Copy Markdown
Contributor

Summary of changes

  • pass on grib root instead of baseline zarr
  • load_forecasts already knows how to handle grib root

@jonasbhend jonasbhend requested a review from cosunae May 27, 2026 14:45
@jonasbhend
Copy link
Copy Markdown
Contributor Author

I am now trying this out to make sure it works. If it does, could you @cosunae please add it to your PR (or we merge in sequence, once unify... is merged).

@jonasbhend
Copy link
Copy Markdown
Contributor Author

ok, this is now ready. Works with the interpolators-ich1.yaml config when enabling meteograms as well.

@jonasbhend jonasbhend changed the title Read baselines from grib instead of zarr for meteograms Fix: read baselines from grib instead of zarr for meteograms May 28, 2026
@frazane
Copy link
Copy Markdown
Contributor

frazane commented May 28, 2026

Let's merge this! Thanks a lot @jonasbhend

@jonasbhend jonasbhend force-pushed the fix_for_meteogram branch from ac2a42c to c4843f9 Compare May 28, 2026 09:05
@jonasbhend jonasbhend force-pushed the fix_for_meteogram branch from c4843f9 to 0b610e1 Compare May 28, 2026 09:08
ds = ds.drop_vars(["U_10M", "V_10M"])
if param == "SP":
ds[param] = ekm_wind.speed(ds.U, ds.V)
ds[param] = (ds.U**2 + ds.V**2) ** 0.5
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

the use of ekm was implemented as part of a review, is it not working?

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

no it is not working because earthkit wants a numpy array, but we are giving it a xarray DataArray. we could probably just change this to ekm_wind.speed(ds.U.values, ds.V.values) and then it should work. Should I revert that?

@jonasbhend jonasbhend merged commit b87dcf6 into unify_regions_stations May 28, 2026
@cosunae cosunae self-requested a review May 28, 2026 09:24
Copy link
Copy Markdown
Contributor

@cosunae cosunae left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

lgtm

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants