pro norp_clb_sn $
  ,mclbx,nclb,npclb,nintclb,stat,vatt,okclb,integ $
  ,timej,timejclb_m,amp,sigma_m $
  ,iatt=iatt,freq=freq

;+
; NAME:
;       NORP_CLB_SN
;
; HISTORY:
;       2002-9-02 TY fix; wrong ATT status detection
;
;-

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

whr=[-1]
for mclb=0,mclbx-1 do begin
    whr=[whr,lindgen(nclb(1,mclb)-nclb(0,mclb))+nclb(0,mclb)]
endfor
whr=whr(1:*)

timejclb=timej(whr)
ampclb=amp(whr)
statclb=stat(whr)
integclb=integ(whr)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

npa=norp_gt_npa(statclb,iatt=iatt,freq=freq)

ficlb=vatt(npa)*float(ampclb)

sigma_m=fltarr(mclbx)
timejclb_m=replicate(timejclb(0),mclbx)
okclb=bytarr(mclbx)

nclb_st=0 & nclb_ed=-1
for mclb=0,mclbx-1 do begin
  nclb_st=nclb_ed+1
  nclb_ed=nclb_st+(nclb(1,mclb)-nclb(0,mclb))-1

  whr=where((npa(nclb_st:nclb_ed) ne npclb),count)
  if (count ge 1) then okclb(mclb)=0 else okclb(mclb)=1
  whr=where((integclb(nclb_st:nclb_ed) ne nintclb),count)
  if (count ge 1) then okclb(mclb)=0 

  timejclb_m(mclb) $
   =anytim((anytim(timejclb(nclb_st))+anytim(timejclb(nclb_ed)))/2.0d,/ints)

  ave=total(ficlb(nclb_st:nclb_ed))/total(integclb(nclb_st:nclb_ed))
  sigma_m(mclb)= $
    sqrt(total((ficlb(nclb_st:nclb_ed)/integclb(nclb_st:nclb_ed)-ave)^2) $
                 /(nclb_ed-nclb_st))*sqrt(float(nintclb))

endfor

return
end