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

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

im3 = im1 - im2;

subplot(2,1,1)

pcolor(im1)

subplot(2,1,2)

pcolor(im2)

subplot(2,1,2)

pcolor(im3)
