Replies: 2 comments 2 replies
-
Hi @jamaa, You are not doing anything wrong, but sadly sdba is currently implemented with datetime grouping in mind. I thought it would still work with an integer coordinate like your Without modifying xclim, it would work if somehow your i.e. in a better world you should do able to do As a side note, the first error you gave is not the real error, that's just the consequence of an internal error on the |
Beta Was this translation helpful? Give feedback.
-
Hi @aulemahal , thanks for the quick response. I see. But if I understand you correctly, the following should work (my first attempt from above, similar to the example in the docs): group = sdba.Grouper("time.dayofyear", add_dims=['number'])
loci = sdba.LOCI.train(da_observed, da_hindcast, thresh="1 mm/d", group=group) But this also gives me an error. Here is the full stacktrace
|
Beta Was this translation helpful? Give feedback.
Uh oh!
There was an error while loading. Please reload this page.
Uh oh!
There was an error while loading. Please reload this page.
-
Setup Information
Context
I am trying to use the Local Intensity Scaling (LOCI) bias-adjustment method on an ensemble forecast of daily precipitation, by training it on hindcasts. I want the adjustment factors to be computed per lead time (no. of days -
step
), but across all realizations of the ensemble (number
).Steps To Reproduce
Result:
In the code above, "time.dayofyear" is a proxy for lead time, ideally I want to use the coordinate "step" instead, but I am not sure how to achieve that? I've also tried the following:
which results in
I've been fiddling around with this for days now and would appreciate if anybody has some ideas or hints, maybe I am doing it completely wrong?
Beta Was this translation helpful? Give feedback.
All reactions