pro plot_hk_info,filename
;NAME
;	plot_hk_info
;PURPOSE
;	To plot the data produced by ave_hk_info.
;RESTRICTIONS
;	Ave_hk_info needs to have been run and a file
;	named craigdata or a user's filename needs to exist from that run,
;SAMPLE CALLING SEQUENCE
;	plot_hk_info
;HISTORY
;	written by Dan Goldberg   25-Aug-95

  if n_elements(filename) eq 0 then begin
     input,'Please enter the filename you used ',filename,'ave_info.txt'
  endif
   spawn, 'grep "number of mnemonics:" ' + filename,value
   nummne=strmid(value,20,4)
   nummne=strcompress(nummne)
   nummne=fix(nummne)
   spawn, 'grep "mnemonics used:" '+filename, return
   farside=16+nummne(0)*8
   mne=strmid(return,16,farside)
   mnes=str2arr(mne, delim=' ')
   namestr='dummy1,dummy2,dummy3,mean,min,max,strd dev'
   names=str2arr(namestr)
   for i=0,nummne(0)-1 do begin
      cmd1='grep "'+mnes(i)+' "  $MDI_CAL_INFO/hk_filter_info.txt'
      spawn, cmd1, tempdescr
      start=strpos(tempdescr(0),'descr=')
      descr=strmid(tempdescr(0),start(0)+6,75)
      descr=strtrim(descr,2)
      cmd='grep "'+mnes(i)+' " '+filename
      spawn, cmd,tempsarr
      sarr=strarr(n_elements(tempsarr)-1)
      sarr=tempsarr(1:n_elements(tempsarr)-1)
      sarr=strmid(sarr,0,90)
      farr=dblarr(7,n_elements(sarr))
      reads, sarr, farr
      !p.multi=[0,2,2]
      for j=3,6 do begin
	  utplot, farr(2,*), farr(j,*), '1-jan-91',/ynozero,  $
		ytitle=names(j),title=mnes(i)+': '+descr,charsize=0.7,psym=1
          average=total(farr(j,*))/n_elements(farr(j,*))
      tempmin=string(min(farr(j,*)),format='(F10.3)')
      tempmax=string(max(farr(j,*)),format='(F10.3)')
      tempstrdev=string(stdev(farr(j,*)),format='(F10.3)')
      tempaverage=string((total(farr(j,*))/n_elements(farr(j,*))),format='(F10.3)')
      minstr='min: '+strtrim(tempmin,2)
      maxstr='max: '+strtrim(tempmax,2)
      strdevstr='strd dev: '+strtrim(tempstrdev,2)
      meanstr='mean: '+strtrim(tempaverage,2)
      spawn, 'grep "'+mnes(i)+' " $MDI_CAL_INFO/hk_filter_info.txt',filtered_info
      minval=strmid(filtered_info(0),strpos(filtered_info(0),'min=')+4,7)
      maxval=strmid(filtered_info(0),strpos(filtered_info(0),'max=')+4,7)
      filteredstr='filter_min: '+strtrim(minval,2)+'  filter_max: '+strtrim(maxval,2)
      percent_filtered=(total(farr(1,*))*100)/total(farr(0,*))
      filteredstr2='percent filtered: '+strtrim(string(percent_filtered,format='(F10.3)'),2)
    ;  filteredstr='filtered: '+strtrim(string(total(farr(1,*)),format='(F10.3)'),2)
      xcoord=0.05
      ycoord=0.95
      plottimes,xcoord,ycoord,meanstr, 0.6
      plottimes,xcoord,ycoord-0.04,minstr,0.6
      plottimes,xcoord,ycoord-0.08,maxstr,0.6
      plottimes,xcoord,ycoord-0.12,strdevstr,0.6
      plottimes,xcoord,ycoord-0.16,filteredstr,0.6
      plottimes,xcoord,ycoord-0.20,filteredstr2,0.6
      endfor
      plottime
      xyouts, 0, .005, 'plot_hk_info ver 1.0', /normal, charsize=.75
      pprint
      pause
   endfor
end