% Plot T in triangles
for ii = 1:Nt,
  for nn = 1:3,
    xtri(nn,ii)  = xy(1,tri2nod(nn,ii));
    ytri(nn,ii)  = xy(2,tri2nod(nn,ii));
    mTtri(nn,ii) = meanT(tri2nod(nn,ii));
    sTtri(nn,ii) = sdevT(tri2nod(nn,ii));
  end
end
subplot(211);
HT = patch(xtri,ytri,mTtri);
axis('equal');
set(HT,'LineStyle','none');
title('mean Temperature');
caxis([900,1300]);
HC = colorbar;
hold on; bladeplotgrid; hold off;

subplot(212);
HT = patch(xtri,ytri,sTtri);
axis('equal');
set(HT,'LineStyle','none');
title('standard deviation of Temperature');
HC = colorbar;
hold on; bladeplotgrid; hold off;