Dealing with HDF format in Matlab (SeaWiFS application)

Ocean color（SeaWiFS）のデータは全てHDFで保存されているので、一昨日論文提出ギリギリでなんとかClimatologyなクロロフィルデータのきれいな図を作りたかったので以下のフォーラム＆MatlabのHDF4の説明を参考に以下のスクリプトを作った。
http://oceancolor.gsfc.nasa.gov/forum/oceancolor/topic_show.pl?tid=3285

clear all;

file={'win' 'spr' 'sum' 'aut'};

abc='abcdefg';

figure;

for it=1:4

    subplot(2,2,it);

    data=hdfread(file{it},'l3m_data' );

    k=find(eq(data,-32767));

    data(k)=NaN;

    data=flipud(data);

    [h1 h2]=contourf(data(1640:1850,1450:1630),80);

    set(h2,'edgecolor','none');

    %caxis([0 1]);

    colorbar('SouthOutside');

    set(gca,'XTick',[]);

    set(gca,'YTick',[]);

    title(abc(it));

end