Skip to content
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

Problem with Cartopy vectors on various projections #2267

Closed
valeriupredoi opened this issue Aug 19, 2021 · 5 comments
Closed

Problem with Cartopy vectors on various projections #2267

valeriupredoi opened this issue Aug 19, 2021 · 5 comments
Labels

Comments

@valeriupredoi
Copy link
Contributor

Not sure if this affects us or not, but better safe than sorry - thanks to @ajheaps who flagged this to me SciTools/cartopy#1179 - can someone who knows their diagnostic plots vectors and uses various projections with Cartopy have a look and see if their results are wonky, please? 🍺

@zklaus
Copy link

zklaus commented Sep 1, 2021

The only use of quiver in our code base that I could find is in esmvaltool/diag_scripts/emergent_constraints/lif1f2.py. There, it seems like this is not an issue because plot and data use the same projection (plate carrée). However, maybe this should be double checked? @katjaweigel and @axel-lauer have worked on that program. Could you two have a quick look?

@katjaweigel
Copy link
Contributor

katjaweigel commented Sep 13, 2021

Yes, it uses the same projection (and at least the current map is rather close to the Equator than to the Pole), so this shouldn't be an issue in esmvaltool/diag_scripts/emergent_constraints/lif1f2.py

@zklaus
Copy link

zklaus commented Sep 13, 2021

Thanks @katjaweigel, closing. @valeriupredoi, please reopen if you still have doubts.

@zklaus zklaus closed this as completed Sep 13, 2021
@valeriupredoi
Copy link
Contributor Author

all good, was just about to do the same, cheers @katjaweigel 🍺

@lguez
Copy link

lguez commented Apr 21, 2023

I think it is an issue for your code in lif1f2.py, even if the projection for axx is PlateCarree. quiver expects u /cos(phi) as an argument and you feed it with u. So I think you get wrong direction and wrong magnitude.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

4 participants