html_url,issue_url,id,node_id,user,created_at,updated_at,author_association,body,reactions,performed_via_github_app,issue https://github.com/pydata/xarray/issues/2499#issuecomment-431717076,https://api.github.com/repos/pydata/xarray/issues/2499,431717076,MDEyOklzc3VlQ29tbWVudDQzMTcxNzA3Ng==,1217238,2018-10-22T00:25:37Z,2018-10-22T00:25:37Z,MEMBER,"I think the problem is actually this line: `da_precip.groupby('point').apply(spi_gamma` The problem is that `spi_gamma` accesses/sets `data_array.values`, which mean it's loading data into NumPy arrays. You need to work entirely with `.data` / dask arrays in order for dask to provide useful speedups. I don't know what `indices.spi` is doing, but you can probably parallelize it either with `dask.array.map_blocks` or xarray's `apply_ufunc`.","{""total_count"": 1, ""+1"": 1, ""-1"": 0, ""laugh"": 0, ""hooray"": 0, ""confused"": 0, ""heart"": 0, ""rocket"": 0, ""eyes"": 0}",,372244156