add in the first plotting routine for Fig.1.