pro mrp_rd_avg,arg0,arg1,arg2,arg3,arg4,freq=freq,avgdir=avgdir

;+
; NAME:
;       NORP_RD_AVG
;
; PURPOSE:
;       This procedure is for reading and daily averaged NoRP data.
;
; CALLING SEQUENCE:
;       norp_rd_avg,st_day,ed_day,day,fiavg,fvavg
;       norp_rd_avg,day_in,day,fiavg,fvavg
;
; INPUTS:
;
; OPTIONAL INPUT KEYWORDS:
; OUTPUTS:
; HISTORY:
;-


if n_elements(avgdir) eq 0 then avgdir=getenv('NORP_AVG')

case n_params() of
  5 : begin                         ; "st_time,ed_time"
      st_time=arg0
      ed_time=arg1
  end
  4 : begin                         ; "date"
      st_time=arg0
      ed_time=anytim(st_time,/ints)
      ed_time.day=ed_time.day+1
  end
  else : begin
  end
endcase

freq=[17]
mx=2
mxf=strtrim(mx,2)

st_daytim=anytim(st_time,/ints)
ed_daytim=anytim(ed_time,/ints)

ndayx=ed_daytim.day-st_daytim.day+1
ex=anytim(transpose([[lonarr(ndayx)],[st_daytim.day+lindgen(ndayx)]]),/ex)
yyyy=reform(ex(6,*))
mm=reform(ex(5,*))
dd=reform(ex(4,*))
yy=yyyy-yyyy/100*100
yymm=string(yy*100+mm)
if (n_elements(yymm) ge 2) then yymm=yymm(uniq(yymm))
yymmdd=$
  string(yymm,format="(i4.4)")
;yymmdd=$
;  string(transpose([[yy],[mm]]),format="(i2.2,i2.2)")
;year='19'+string(yy,format="(i2.2)")
  
if (n_elements(yymmdd) ge 2) then yymmdd=yymmdd(uniq(yymmdd))
files='mav'+yymmdd
dir=avgdir
infil=file_list(dir,files)
  if (n_elements(infil) eq 1) then begin
  if (infil(0) eq '') then begin
    print,'No files found'
    return
  endif
  endif

mfilx=n_elements(infil)

mfil=0
text0=rd_tfile(infil(mfil),11,nocomment='!')
dummy=0
if (text0(0) ne '') then text=text0 else begin
  text=strarr(11,2)
  dummy=2
endelse
if ((size(text))(0) eq 1) then begin
  text=strarr(11,2)
  text(*,1)=text0
  dummy=1
endif


if (mfilx ge 2) then begin
 for mfil=1,mfilx-1 do begin
  text0=rd_tfile(infil(mfil),11,nocomment='!')
 if (text0(0) ne '') then begin
  if ((size(text0))(0) eq 1) then text=[[text],[text0]] else $
    text=concat2d(text,text0)
 endif
 endfor
endif

if (dummy ne 0) then text=text(*,dummy:*)


sz=size(text)

date=reform(text(0,*))
date=anytim(date,/ints)

whr=where(((date.day ge st_daytim.day) $
       and (date.day lt ed_daytim.day)) ,cnts)

if (cnts eq 0) then begin
  print,'No observation in this duration'
  goto,Lend
endif

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;  sort & uniq

date=date(whr)
avg1=float(reform(text(1:*,whr)))

nday0=n_elements(date)

  date=date(0:nday0-1)
  avg1=avg1(*,0:nday0-1)
  fiavg1=avg1(2*lindgen(mx/2),0:nday0-1)
  fvavg1=avg1(2*lindgen(mx/2)+1,0:nday0-1)
  fiavg=fltarr(mx/2+1,nday0)
  fiavg(0:mx/2-1,*)=fiavg1
  fvavg=fltarr(mx/2+1,nday0)
  fvavg(0:mx/2-1,*)=fvavg1
if (nday0 eq 1) then begin
  date=date(0)
  avg1=reform(avg1)
  fiavg1=avg1(2*lindgen(mx/2))
  fvavg1=avg1(2*lindgen(mx/2)+1)
  fiavg=fltarr(mx/2+1)
  fiavg(0:mx/2-1)=fiavg1
  fvavg=fltarr(mx/2+1)
  fvavg(0:mx/2-1)=fvavg1
endif

case n_params() of
  5 : begin                         ; "st_time,ed_time"
      arg2=date
      arg3=fiavg
      arg4=fvavg
  end
  4 : begin                         ; "date"
      arg1=date
      arg2=fiavg
      arg3=fvavg
  end
  else : begin
  end
endcase
Lend:

return
end