function dP=Pdrop2(Q)

load welldata2

X=Q(3);
if X > Xo(end);
    X=Xo(end);
end

Y=Q(2);
if Y > Yv(end);
    Y=Yv(end);
end

Z=Q(1);
if Z > Zg(end);
    Z=Zg(end);
end

dP = interp3(Xo,Yv,Zg,Pddat,X,Y,Z);         % Linear interpolation between datapoints
