
dfname = ['fort.89'];
fid=fopen(dfname);
im1 = fread(fid,[116 300],'float'); status = fclose(fid);
pcolor(im1)
