; ; check intensity ; mp1=fltarr(2034,2034) mp2=mp1 mp3=mp1 mp4=mp1 close,1 & openr,1,'mp1' & readu,1,mp1 close,1 & openr,1,'mp2' & readu,1,mp2 close,1 & openr,1,'mp3' & readu,1,mp3 close,1 & openr,1,'mp4' & readu,1,mp4 img = (mp1+mp2+mp3+mp4)/4. img = img - 20. limg = alog10(img) limg = limg-10 merge = mp1 close,1 & openr,1,'Imerge' & readu,1,merge merge = alog10(merge) merge = merge-10 !y.range=[-10,-5] plot,limg(*,1017) oplot,merge(*,1017) stop end