function ff = makerr (R);
smooth_level =90; %!!!!!
for i=1:288
   p5=polyfit([1:72],R(i,:),5);
   R(i,:)=polyval(p5,[1:72]);
end;
 re= min(R');
 for i=1:288
  ff(i) = find(R(i,:)==re(i));  
 end; %i
 ff= smooth([ff ff ff],smooth_level,'lowess');
 %plot(ff);
 ff = ff(289:288*2);
figure;
contourf(R,10);
hold on;
plot(ff,[1:288],'--r','linewidth',2.5);
