%delRny= [50.,  70., 100., 140., 190.,240., 290., 340., 390., 440.,490.,...
%540., 590., 640., 690.,];

delRny= [14.8461, 17.7881, 21.3130, 25.5365, 30.5969, 36.6601, 43.9247,...
   52.6290, 63.0581, 75.5539, 90.5259, 108.4648, 129.9585, 155.7114,...
  186.5677, 223.5385, 267.8356, 320.9107, 384.5033, 460.6977, 551.9909,...
  661.3751, 792.4352, 949.4666];

t_ref = fitTemp(delRny);