pro cat, fn if n_elements(fn) le 0 then fn = '/Users/kw/Desktop/NORM_jelly_intensity_plot.dat' ; if n_elements(fn) le 0 then fn = '/Users/kw/Desktop/jelly_intensity_plot.dat' if !d.name eq 'X' then begin device, decomposed=0 endif if !d.name eq 'PS' then begin device, /color , bits=8 endif d=read_ascii(fn) d=d.(0) print, size(d) x=reform(d(0,*)) y=reform(d(1,*)) z=reform(d(2,*)) loadct, 0 plot, x, y, psym=3,/isotropic,/nodata $ ,xstyle=3,ystyle=3,xrange=[-22.5,22.5],yrange=[-22.5,22.5] xlhs =!x.crange(0) xrhs =!x.crange(1) ylhs =!y.crange(0) yrhs =!y.crange(1) xbox=[xlhs,xrhs,xrhs,xlhs] ybox=[ylhs,ylhs,yrhs,yrhs] ; polyfill, /data, xbox, ybox, color=0 ; use polyfill to give black background ; choose x/ystyle=3 to have axes, 5 without ; make region of polyfill smaller to not overlap y-axis print,'testing' loadct, 8 ; call program rainbow.pro for our colortable, setting grey=240 ; rainbow,247 const = 0. ; if tabulated z, const=0. If tabulated (1+z), const=1 zoff=z-const ; offset z, if need be b = bytscl(z^(1./4.)) for i=0,255 do begin w=where(b eq i, cnt) if(cnt gt 0) then begin oplot,x(w),y(w),psym=3 $ ,color = i endif endfor ; close program properly, all colours will show properly in .ps if !d.name eq 'PS' then device,/close end