This is a helper function to get a vector of 3-character province codes.
getCHMDprovinces()
a character vector of 3 character province codes.
## Not run: (provs <- getCHMDprovinces())