New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Meshfill plot with Mercator projection issue #1671

Open
sankhesh opened this Issue Nov 10, 2015 · 2 comments

Comments

Projects
None yet
2 participants
@sankhesh
Contributor

sankhesh commented Nov 10, 2015

This use case currently fails. Test case provided in 9d6bed6

Run ctest -V -R vcs_test_animate_projected_meshfill_mercator

Error Reported:

File "<uvcdat>/install/lib/python2.7/site-packages/cdms2/coord.py", line 289, in subSlice
newbounds = bounds[specs]   # bounds can be a numarray or DatasetVariable
IndexError: only integers, slices (`:`), ellipsis (`...`), numpy.newaxis (`None`) and integer or boolean arrays are valid indices
@doutriaux1

This comment has been minimized.

Member

doutriaux1 commented Nov 10, 2015

@sankhesh I fixed @aashish24 code so that it now passes (branch: issue_1671_meshfill_mercoator). But there is a more fundamental issue here, cdms2 does not "shrink" curvilinear grids when subsetting but rather masks out cells that are not within that range, so @aashish24 code does not help in this case.

It maybe better to copy the graphic method and set the limits ourselves, if you do so it works
in the test I add simply do:

m.datawc_y1=-85
m.datawc_y2=85

Then it works and plots:
meshfill_mercator

@doutriaux1

This comment has been minimized.

Member

doutriaux1 commented Nov 10, 2015

@aashish24 @chaosphere2112 @williams13 not sure how to proceed from here, suggestions welcome.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment