pn='/homes/dmenemen/Arctic36km/MITgcm/exe/'; nx=210; ny=192; nme='SIheff'; for ts=36:36:262980 fn=[pn nme '.' myint2str(ts,10) '.data']; if exist(fn) fld=readbin(fn,[nx ny]); clf mypcolor(fld'); caxis([0 5]) thincolorbar title(['nme ' ts2dte(ts,2400)]) pause(.01) else break end end