function Qw3=wpcn3(Pw3,well3)
x=well3(:,1);
Y=well3(:,2:4);
Qw3= interp1(x,Y,Pw3)';