;
; WRITTEN:  Xuepu Zhao May 10, 2007
; MODIFIED: Xuepu Zhao May 23, 2009
;
pro updbccm_crnclog,crn,clog,itl,bccm0,bccmr,bccml,updbccm,updscrncld,SINL=sinl,hwd

    if n_params( ) lt 1 then begin
       print,'pro updbccm_crnclog,crn,clog,itl,updbccm,/sinl,hwd'
       print,'clog:Interested Carrington longitude in degrees'
       print,'itl=60_Left, 180_Center, itl and hwd in degrees'
       print,'bccm=B component Carrington Map, can be Bl or Br CM'
       return
    endif
  ; dirp='/scr20/xuepu/MDIBRSC/phof360_type1/' ;It is changeable
    scrn0=strtrim(crn,2)
    sz=SIZE(bccm0) 
    xsz=sz(1) & ysz=sz(2)                            ; in pixels
    updbccm=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)
    updscrncld=strarr(xsz)
    dlogitl=(clog-itl)*ppd                           ; in pixels
    CASE 1 OF
     dlogitl lt 0: begin
       scrnl=strtrim(crn+1,2)
       updbccm(0:(-dlogitl-1),*)=bccml((dlogitl+xsz):(xsz-1),*)
       updbccm(-dlogitl:(xsz-1),*)=bccm0(0:(xsz-1+dlogitl),*)
       for i=0,xsz-1 do scrnphdl(i)=scrnl+':'+sphd(i)
       updscrncld(0:(-dlogitl-1))=scrnphdl((dlogitl+xsz):(xsz-1))
       updscrncld(-dlogitl:(xsz-1))=scrnphd0(0:(xsz-1+dlogitl))
                  end
     dlogitl eq 0: begin
       updbccm=bccm0
       updscrncld=scrnphd0
                    end
     dlogitl gt 0: begin
       scrnr=strtrim(crn-1,2)
       updbccm(0:(xsz-1-dlogitl),*)=bccm0(dlogitl:(xsz-1),*)
       updbccm((xsz-dlogitl):(xsz-1),*)=bccmr(0:(dlogitl-1),*)
       for i=0,xsz-1 do scrnphdr(i)=scrnr+':'+sphd(i)
       updscrncld(0:(xsz-1-dlogitl))=scrnphd0(dlogitl:(xsz-1))
       updscrncld((xsz-dlogitl):(xsz-1))=scrnphdr(0:(dlogitl-1))
                     end
     else: 

    ENDCASE
    if keyword_set(sinl) then begin
       tek_color
       ofl='updbccm_'+STRTRIM(crn,2)+'_'+STRTRIM(clog,2)+'_'+STRTRIM(itl,2)+'.ps'
       psout,012,ofl
       xticsc_crnclog09,crn,clog,itl,xticscb,xticsct
       bupdbccm=BYTSCL(updbccm,max=60, min=-60)
       pscimg06, bupdbccm,' ',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