Replies: 2 comments 2 replies
-
Hi @rsx23ser , I appreciate that's frustrating. We'd need much more information though — what's a simpler example that works? What's happens if you run over smaller data? Does it work if the arrays are numpy arrays? See the links in #5404 for more |
Beta Was this translation helpful? Give feedback.
0 replies
-
@max-sixty Thanks a lot for your response! I will try to create a simpler example and repost. |
Beta Was this translation helpful? Give feedback.
2 replies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
-
I am trying to reproduce the map_blocks example given here (link), which finds the mixed-layer depth. I am using output files from the POP ocean model but cut up into smaller lon/lat boxes (using ncks), as opposed to reading in the global fields. The only change I have made to the subroutine in that link is that I am directly feeding in the density to the mld_dsigma routine instead of computing the density from salinity and temperature.
`def mld_dsigma(RHO, dsigma=0.03, rho_chunks='auto'):
The other routine interp_mld_dsigma is unchanged from what is in that link. When I execute the above routine, for some reason the calculation just stops after a point. I have uploaded a screenshot of the dask dashboard (link) which shows the point beyond which the code does not appear to proceed. Eventually I just kill the job. Has anybody experienced a similar issue before? Thanks for any help.
Beta Was this translation helpful? Give feedback.
All reactions