added a few comments and made use of plot labelling functions consistent
extract all variable values before subsetting - runs much faster this way!