function to subset meshblocks takes ages so use spatialData.Rdata
fixed typo, meshblocks now subset correctly
plot DCBA