clear;
switchtolayout;
setnamed("HEAT","enabled",1);
mesh("HEAT");
setnamed("HEAT","mesh lock",1);

runsweep("Pin");

T=getsweepresult("Pin","T");
Pin=T.scale/3.87147e-07;
delta_T = matrix(length(Pin));
Pin_array=size(T.T);
Pin_ind=find(Pin_array,length(Pin));

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

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

plot(Pin,norm_factor*delta_T,"Pin (mW)","delta_T (K)");