Rooki calls climate data operations on the rook processing service.
WARNING: used cmip5 demo data instead of c3s-cmip5.
from rooki import operators as ops
wf = ops.Diff(
ops.Subset(
ops.Input(
'inm', ['cmip5.output1.INM.inmcm4.rcp45.mon.ocean.Omon.r1i1p1.latest.zostoga']
),
time="2085-01-01/2120-12-30",
),
ops.Subset(
ops.Input(
'mpi', ['cmip5.output1.MPI-M.MPI-ESM-LR.rcp45.mon.ocean.Omon.r1i1p1.latest.zostoga']
),
time="2085-01-01/2120-12-30",
),
)
resp = wf.orchestrate()
resp
metalink output
print(resp.xml)
URLs in metalink document ...
resp.download_urls()
... and open with xarray
ds = resp.datasets()[0]
ds