code to read and plot the data; need to extract subset  for us west