; ; Plot graphs of polar and equatorial intensities ; xcen = 1. ycen = 1. poise_scale=1. pangle=1. close,1 & openr,1,'centers' & readf,1,xcen,ycen,pangle,poise_scale & close,1 print,xcen,ycen,pangle,poise_scale ;!p.title='POISE 98 pB0 vs Height' ;!y.title='log pB0' !p.title='POISE 98 B0 vs Height' !y.title='log B0' !x.title='Solar Radius' !p.charsize=2 !x.ticklen=1 !y.ticklen=1 !p.thick=2 !p.font=1 !y.range=[-11,-6] pangle = -20.9 img=fltarr(2034,2034) ;close,1 & openr,1,'pBmerge' & readu,1,img & close,1 close,1 & openr,1,'Imerge' & readu,1,img & close,1 img=rot(img,pangle,1.,xcen,ycen,/interp) fred=rebin(img,339,339) ;this will average data fred=alog10(fred*1.e-10) height=findgen(339)-339./2. height=height*3.1*6./960. plot,height(220:338),fred(220:338,170) oplot,height(220:338),fred(170,220:338),line=2 ;xyouts,2.4,-8.5,'Equator' ;xyouts,1.9,-9.8,'Pole' xyouts,2.4,-7.5,'Equator' xyouts,1.9,-8.8,'Pole' stop end