function Qw1=wpcn1(Pw1,well1)
x=well1(:,1);
Y=well1(:,2:4);
Qw1= interp1(x,Y,Pw1)';
 
