xsz=2034 ysz=2034 if(n_elements(NEW) eq 0 ) then begin ;{ OLD=fltarr(xsz,ysz) & NEW=fltarr(xsz,ysz) endif ;} ofilenm='OPALS/ORIG/opal_L' & nfilenm='OPALS/opal_L' ;; ofilenm='OPALS/ORIG/opal_M' & nfilenm='OPALS/opal_M' ;; ofilenm='OPALS/ORIG/opal_S' & nfilenm='OPALS/opal_S' close,1 & openr,1,ofilenm & readu,1,OLD & close,1 close,1 & openr,1,nfilenm & readu,1,NEW & close,1 OLD=OLD*1.e12 NEW=NEW*1.e12 diffi=(OLD-NEW) plot,[100,1000],[-30,150],/nodata,ystyle=1 oplot,OLD(100:1000,400) oplot,NEW(100:1000,400) oplot,diffi(100:1000,400) end