%---------------

function [Zo1,Z11] = Z(Pr1,Tr1)
global Pr Tr Zo Z1
Zo1 = interp2(Pr,Tr,Zo,Pr1,Tr1);
Z11 = interp2(Pr,Tr,Z1,Pr1,Tr1);

%--------------