flname = '/share/norp/xdr/2011/02/norp20110215_0154.xdr'
restore, flname
mint=10
norp_mkint,mint,mvd,tim,fi,fv,mvdav,timav,fiav,fvav
for m=0,6 do fiav(m,*)=fiav(m,*)-fiavg(m)

; cut out a portion of mint * nspan sec. starting nstrt sec. from the begining of the daily data
if (1 eq 2) then begin
  nspan=400
  nstrt=300
  fiavtmp=fiav 
  fiav =fltarr(7,nspan)
  mvdavtmp=mvdav 
  mvdav=bytarr(7,nspan)
  for m=0,6 do begin
    fiav(m,*) =fiavtmp (m,nstrt:nstrt+nspan-1)
    mvdav(m,*)=mvdavtmp(m,nstrt:nstrt+nspan-1)
  endfor
endif

if (1 eq 1) then begin
  set_plot,'ps'
  device, filename='freq3_094_170_340.ps'
  norp_plot,3,mvdav,timav,fiav,title='3 freq. 9.4, 17.0 and 34.0'
  norp_plot,4,mvdav,timav,fiav,/over
  norp_plot,5,mvdav,timav,fiav,/over
;  plot,  fiav(3,*), title='3 freq. 9.4, 17.0 and 34.0'
;  oplot, fiav(4,*)
;  oplot, fiav(5,*)
  device,/close
  set_plot,'x'
endif
if (1 eq 1) then begin
  norp_plot,3,mvdav,timav,fiav,title='3 freq. 9.4, 17.0 and 34.0'
  norp_plot,4,mvdav,timav,fiav,/over
  norp_plot,5,mvdav,timav,fiav,/over
endif

if (1 eq 1) then begin
  numdat=n_elements(fiav(0,*))
  freq2=freq(3:5)
  freq3=alog(freq2)
  lncoef=fltarr(2,numdat)
  nsweep=1
  lncoefsweep=fltarr(2,numdat,(nsweep*2+1)*(nsweep*2+1)*(nsweep*2+1))
  for i3 = -nsweep, nsweep do begin
  for i4 = -nsweep, nsweep do begin
  for i5 = -nsweep, nsweep do begin
    fiav2=fiav(3:5,*)
    fiav2(0,*)=fiav2(0,*)+float(i3)*fiavg(3)*0.025
    fiav2(1,*)=fiav2(1,*)+float(i4)*fiavg(4)*0.025
    fiav2(2,*)=fiav2(2,*)+float(i5)*fiavg(5)*0.025
    fiav3=alog(fiav2)
    for n = 0L, numdat - 1 do begin
      sgmt=fiav3(*,n)
      adrs=where(finite(sgmt,/nan))
      if (adrs[0] lt 0) then begin
        results=linfit(freq3,sgmt)
      endif else begin
        results=[0.0,0.0]
      endelse 
      if ((i3 eq 0) and (i4 eq 0) and (i5 eq 0)) then lncoef(*,n)=results
      isweep=(i3+nsweep)*(nsweep*2+1)*(nsweep*2+1)+(i4+nsweep)*(nsweep*2+1)+(i5+nsweep)
      lncoefsweep(*,n,isweep)=results
    endfor
  endfor
  endfor
  endfor

  lncoefmax=lncoef
  lncoefmin=lncoef
  for n = 0L, numdat - 1 do begin
    lncoefmax(0,n)=max(lncoefsweep(0,n,*))
    lncoefmax(1,n)=max(lncoefsweep(1,n,*))
    lncoefmin(0,n)=min(lncoefsweep(0,n,*))
    lncoefmin(1,n)=min(lncoefsweep(1,n,*))
  endfor

  set_plot,'ps'

  device, filename='linfitcoef0.ps'
  norp_plot,0,mvdav,timav,lncoef,yrange=[-5,15],title='linear fit coef.0'
  if (1 eq 2) then begin
    for i3 = -nsweep, nsweep do begin
    for i4 = -nsweep, nsweep do begin
    for i5 = -nsweep, nsweep do begin
      isweep=(i3+nsweep)*(nsweep*2+1)*(nsweep*2+1)+(i4+nsweep)*(nsweep*2+1)+(i5+nsweep)
      lncoeftmp=lncoefsweep(*,*,isweep)
      norp_plot,0,mvdav,timav,lncoeftmp,psym=3,/over
    endfor
    endfor
    endfor
  endif else begin
    norp_plot,0,mvdav,timav,lncoefmax,line=3,/over
    norp_plot,0,mvdav,timav,lncoefmin,line=3,/over
  endelse
  device,/close

  device, filename='linfitcoef1.ps'
  norp_plot,1,mvdav,timav,lncoef,yrange=[-4,2],title='linear fit coef.1'
  if (1 eq 2) then begin
    for i3 = -nsweep, nsweep do begin
    for i4 = -nsweep, nsweep do begin
    for i5 = -nsweep, nsweep do begin
      isweep=(i3+nsweep)*(nsweep*2+1)*(nsweep*2+1)+(i4+nsweep)*(nsweep*2+1)+(i5+nsweep)
      lncoeftmp=lncoefsweep(*,*,isweep)
      norp_plot,1,mvdav,timav,lncoeftmp,psym=3,/over
    endfor
    endfor
    endfor
  endif else begin
    norp_plot,1,mvdav,timav,lncoefmax,line=3,/over
    norp_plot,1,mvdav,timav,lncoefmin,line=3,/over
  endelse
  device,/close

  set_plot,'x'

  norp_plot,0,mvdav,timav,lncoef,yrange=[-5,15],title='linear fit coefs'
  norp_plot,1,mvdav,timav,lncoef,yrange=[-4,2], /over
endif else begin
  mvfreq=[4,5,6]
  norp_alpha,freq(mvfreq),fiav(mvfreq,*),mvdav(mvfreq,*),mvdfit,alpha_tk,alpha_tn,freqpk,fluxpk
  norp_alpha2delta,alpha_tn,delta
  print, max(delta), min(delta)
endelse

freq2=freq(3:5)
freq3=alog(freq2)
fiav2=fiav(3:5,*)
fiav3=alog(fiav2)
save,timav,lncoef,freq3,fiav3,filename='cdaw14_norp3_2011feb15_all.sav'
save,timav,filename='timav.sav'
save,lncoef,filename='lncoef.sav'
save,freq3,filename='freq3.sav'
save,fiav3,filename='fiav3.sav'

end