add read harvest function, complete map plotting for yearly output