Sourish Basu

geos 3.9.1, which got installed by macports as a dependency when I installed basemap.

Do I have to go back all the way to 3.3? I can install an older version, but there seems to have been several versions between 3.3 and 3.9.

Downgrading to geos 3.8.1 solved it! Any idea why an upgraded geos causes this? I'd be happy to file a geos bug report if I knew what exactly was causing...

I recently updated my netCDF4 module to 1.5.1.2 and was surprised by this default behavior of always returning masked arrays, since it broke my code. Googling led me here, and...

What does set_auto_mask do to variables in netcdf files that actually have masked elements? If a variable in the file has *some* masked elements (i.e., filled with _FillValue), then I'd...

### Update This problem also occurs with the `EckertIV` [projection](https://drive.google.com/file/d/1frrUh9YGFk72T47tVHABvw-Qh17aTfbm/view?usp=sharing) but **not** with the `EckertIII` [projection](https://drive.google.com/file/d/1o1eTAxZBW3sYM4Ckw5hDoBODtREvBZ6a/view?usp=sharing). Maybe that's a clue?

@greglucas Here's a MWE: ```python import cartopy.crs as ccrs import cartopy.feature as cfeature from matplotlib import pyplot as plt import numpy as np # get the data lats, lons, nobs...

@mathause Indeed, `ax.set_global()` did the trick. However, this is not default matplotlib behavior for `scatter`. E.g., the following code produces [this plot](https://drive.google.com/file/d/1iylZXaGJNZhqzz4MJX7fDYV8ZDjqrLq3/view?usp=sharing): ```python fig = plt.figure() ax = plt.axes() ax.set_xlim(-3.,3.)...

@greglucas Indeed, @mathause has the right solution. However, I'm wondering if automatically resizing the axes is the correct thing to do when plotting data on a map.

@greglucas OK, understood. However, in the case of a scatter plot on a map, it actually **contracted** the limits. Also, turning off autoscaling does not work for cartopy axes. If...