Fix for projections using a spherical earth #336
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Some projections use a spherical earth. This is indicated in the proj4-string by either
the
+R=...
identifier orby setting the
+a=
and+b=
equal.In both cases the
parse_proj4_string
returns a dict with theR
-key, but this key is not checked byproj4_to_cartopy
resulting in the default WGS84 globe for the cartopy CRS, giving small differences in the locations of the gridpoints.This is now fixed by explicitly checking for the
R
-key inproj4_to_cartopy
and defining the correct globe keywords.I needed this to take advantage of the plot_precip_fields function for plotting precipitation coming from the ALARO NWP model.