It's also related to ```iris.coord_systems.GeogCS()``` not having a default value. Perhaps we should consider defaulting to either spherical with the PP value or to ellipsoidal with the wgs84 values (6378137.0 m with 1/f of 298.257). I'll create a module constant in the meantime.