script for plot1