; ; WRITTEN: Xuepu Zhao May 10, 2007 ; MODIFIED: Xuepu Zhao May 23, 2009 ; pro updmsc_crnclog20090806,crn,clog,itl,updmsc,SINL=sinl,hwd if n_params( ) lt 1 then begin print,'pro updmsc_crnclog20090806,crn,clog,itl,updmsc,/sinl,hwd' print,'clog:Interested Carrington longitude in degrees' print,'itl=60_Left, 180_Center, itl and hwd in degrees' return endif dirp='/scr20/xuepu/MDIBRSC/phof360_type1/br360_' ;It is changeable scrn0=strtrim(crn,2) mscf0=dirp+scrn0+'.fits' mscbr0=readfits(mscf0,ht,/silent) ; No data gap sz=SIZE(mscbr0) xsz=sz(1) & ysz=sz(2) ; in pixels updmsc=fltarr(xsz,ysz) ppd=xsz/360. ; X-pixel number of per degree dlogitl=(clog-itl)*ppd ; in pixels CASE 1 OF dlogitl lt 0: begin scrnl=strtrim(crn+1,2) mscfl=dirp+scrnl+'.fits' mscbrl=readfits(mscfl,hd,/silent) updmsc(0:(-dlogitl-1),*)=mscbrl((dlogitl+xsz):(xsz-1),*) updmsc(-dlogitl:(xsz-1),*)=mscbr0(0:(xsz-1+dlogitl),*) end dlogitl eq 0: begin updmsc=mscbr0 end dlogitl gt 0: begin scrnr=strtrim(crn-1,2) mscfr=dirp+scrnr+'.fits' mscbrr=readfits(mscfr,hd,/silent) updmsc(0:(xsz-1-dlogitl),*)=mscbr0(dlogitl:(xsz-1),*) updmsc((xsz-dlogitl):(xsz-1),*)=mscbrr(0:(dlogitl-1),*) end else: ENDCASE if keyword_set(sinl) then begin tek_color ofl='updmsc_'+STRTRIM(crn,2)+'_'+STRTRIM(clog,2)+'_'+STRTRIM(itl,2)+'.ps' psout,012,ofl xticsc_crnclog09,crn,clog,itl,xticscb,xticsct bupdmsc=BYTSCL(updmsc,max=60, min=-60) ; pscimg06, bupdmsc,'UpdatedSC with CR'+scrn0+' Clog0='+strtrim(clog,2)+' ITL='+STRTRIM(itl,2),px,py,sy='sla' pscimg06, bupdmsc,' ',px,py,sy='sla' updsc_annot09,1,xticb=xticscb,xtict=xticsct,longi='carrington',/sinl oplot,[itl,itl],[-1,1],thick=4,color=2 oplot,[itl-hwd,itl-hwd]*ppd,[-1,1],thick=2,color=2 oplot,[itl+hwd,itl+hwd]*ppd,[-1,1],thick=2,color=2 psdone print,ofl endif end