function h=plot_tf_comparisons_0


    mbw_4229_01=(rand(1,14)-0.9).*80;
    mbw_4229_02=(rand(1,14)-0.9).*80;
    mbw_0388_03=(rand(6,14)-0.9).*80;
    mbw_4252_04=(rand(8,14)-0.9).*80;
    mbw_4229_04=(rand(2,14)-0.9).*80;
    mbw_4229_05=(rand(5,14)-0.9).*80;

    clf;
    subplot(221);
    plot(mbw_4229_01(6),mbw_4229_01(13),"r");
    hold on;
    plot(mbw_4229_02(6),mbw_4229_02(13),"g");
    plot(mbw_4229_04(:,6),mbw_4229_04(:,13),"b");
    plot(mbw_4229_05(:,6),mbw_4229_05(:,13),"m");
    plot([-70,10],[-70,10],"color",[0.8,0.8,0.8]);
#    hold off;
    set(gca,"xlim",[-70,10],"ylim",[-70,10]);
    xlabel("MBW frost point (^oC)");
    ylabel(" frost point (^oC)");
    title("Comparison with s/n 4229")
    legend("1^{st} October (unheated?)","2^{nd} October (unheated?)",...
                             "4^{th} October, heated","5^{th} October, heated",...
                             "location","southeast");


    subplot(2,2,2);
    plot([-70,10],[-70,10],"color",[0.8,0.8,0.8]);
    hold on;
    plot(mbw_0388_03(:,6),mbw_0388_03(:,12),"b");
#    hold off;
    set(gca,"xlim",[-70,10],"ylim",[-70,10]);
    xlabel("MBW frost point (^oC)");
    ylabel(" frost point (^oC)");
    title("Comparison with s/n 388")

    h_a(3) = subplot(2,2,3);
    plot([-70,10],[-70,10],"color",[0.8,0.8,0.8]);
    hold on;
    plot(mbw_4252_04(:,6),mbw_4252_04(:,12),"b");
#    hold off;
    set(gca,"xlim",[-70,10],"ylim",[-70,10]);
    xlabel("MBW frost point (^oC)");
    ylabel(" frost point (^oC)");
    title("Comparison with s/n 4252")


endfunction;