In the extract step of my ESPO-G workflow, I add a few steps to handle RDRS.
I want to cut my region on rlat and rlon. From what I understand, clisops only cuts on lat and lon.
I want to add nans from my mask.
Potential Solution
Add a method sel to region that would not be a call to clisops, just
arg_sel = {dim:slice(*map(float, bounds)) for dim, bounds in region['sel'].items()}
ds= ds.sel(**arg_sel )
Add a mask arg to add NaNs. Ideally, I would want to be able to ask for the mask in search_data_catalogs with {mask:fx} in variables_and_freq. Then, if mask is True. Apply mask on the other variables. mask could also be a dataarray if the mask is not in a catalogue.
ds= ds.where(mask == True)
Additional context
Let me know you thoughts on the best way to implement the mask.
Right now, this is all in my ESPO-G workflow, but it is not very clean... It would be better in extract.
Contribution
[X] I would be willing/able to open a Pull Request to contribute this feature.
Addressing a Problem?
In the extract step of my ESPO-G workflow, I add a few steps to handle RDRS.
Potential Solution
selto region that would not be a call to clisops, justmaskarg to add NaNs. Ideally, I would want to be able to ask for the mask insearch_data_catalogswith{mask:fx}invariables_and_freq. Then, ifmaskis True. Apply mask on the other variables.maskcould also be a dataarray if the mask is not in a catalogue.Additional context
Let me know you thoughts on the best way to implement the mask.
Right now, this is all in my ESPO-G workflow, but it is not very clean... It would be better in extract.
Contribution