;
; WRITTEN:  Xuepu Zhao May 10, 2007
; MODIFIED: Xuepu Zhao May 23, 2009
;
pro updmsc_crnclog20090810,crn,clog,itl,dirp,updmsc,updscrnphd,SINL=sinl,hwd

    if n_params( ) lt 1 then begin
       print,'pro updmsc_crnclog20090810,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
    zgrid,xsz,ysz,phd,thd,lad,cth,sth
    sphd=strarr(xsz)
    for i=0,xsz-1 do begin
        CASE 1 OF
         phd(i) lt 10: sphd(i)='00'+strtrim(phd(i),2)
         phd(i) ge 10 and phd(i) lt 100: sphd(i)='0'+strtrim(phd(i),2)
         else: sphd(i)=strtrim(phd(i),2)
        ENDCASE
    endfor
    sphd=STRMID(sphd,0,7)
    scrnphd0=strarr(xsz)
    for i=0,xsz-1 do scrnphd0(i)=scrn0+':'+sphd(i)
    scrnphdl=strarr(xsz)
    scrnphdr=strarr(xsz)
    updscrnphd=strarr(xsz)
    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),*)
       for i=0,xsz-1 do scrnphdl(i)=scrnl+':'+sphd(i)
       updscrnphd(0:(-dlogitl-1))=scrnphdl((dlogitl+xsz):(xsz-1))
       updscrnphd(-dlogitl:(xsz-1))=scrnphd0(0:(xsz-1+dlogitl))
                  end
     dlogitl eq 0: begin
       updmsc=mscbr0
       updscrnphd=scrnphd0
                    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),*)
       for i=0,xsz-1 do scrnphdr(i)=scrnr+':'+sphd(i)
       updscrnphd(0:(xsz-1-dlogitl))=scrnphd0(dlogitl:(xsz-1))
       updscrnphd((xsz-dlogitl):(xsz-1))=scrnphdr(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