clear;
mesh;
setnamed("HEAT","mesh lock",1);
runsweep;

T=getsweepresult("Pin","temperature");
Pin=T.Pin;
delta_T = matrix(length(Pin));

for (i=1:length(Pin)) { delta_T(i) = max(pinch(T.T,3,i))-300;}

norm_factor = pi*(19.5e-6/2)^2/85e-6^2;

plot(Pin,norm_factor*delta_T);