Adding hurs (realtive humidity) derivation script #2397
Conversation
malininae
added
enhancement
| tdps_cube = cubes.extract_cube(NameConstraint(var_name='tdps')) | ||
| tas_cube = cubes.extract_cube(NameConstraint(var_name='tas')) | ||
|
|
||
| cubes_difference = iris.analysis.maths.subtract(tas_cube, tdps_cube) |
There was a problem hiding this comment.
| cubes_difference = iris.analysis.maths.subtract(tas_cube, tdps_cube) | |
| cubes_difference = tas_cube - tdps_cube |
Would this also work? Did you read the section on cube maths in the iris documentation?
There was a problem hiding this comment.
Thanks, so when I was trying it with simple tas_cube-tdps_cube the issue was that I was running out of memory, while it wasn't a case with subtract. I will test it again with the time diagnostic and will let you know.
There was a problem hiding this comment.
That would be surprising because it should do the same thing, i.e. the minus operator on an Iris cube is implemented using the exact same function:
https://github.com/SciTools/iris/blob/c5e0353da11c263656c6c77e5e56eafefcede077/lib/iris/cube.py#L3968
This could probably be done with code like this: import dask.array as da
hurs_cb.data = da.ma.where(hurs_cb.core_data() > 100., hurs_cb.core_data(), 100.) |
|
If you would like a quicker response in the future, it usually helps if you tag people so they get a notification, e.g. @bouweandela |
Hi folks,
Following the April 2024 monthly meeting, I'm adding the relative humidity derivation (
hurs). Before being able to open a real PR, I would need to add filtering in case thehurs_cb.data>100%. Also need to write test. Will do it before the workshop.This PR is related to #2254 .