Weekly_report/weekly_telem.pro at master · chandra-mta/Weekly_report · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
PRO WEEKLY_TELEM, tstart,tstop

; compile telemetry limit violations table for weekly report.
; tstart and tstop are of the form 20060626, yyyymmdd

outfile="telem.txt"
readcol,"telem.lst",msids,format='a'
maxarr=strarr(9,n_elements(msids))
maxarr(8,*)=0 ; mark all as no violations found yet
maxarr(0,*)=msids
minarr=maxarr

date=long(tstart)
day=1
while (date le tstop) do begin
  date_str=strcompress(string(date),/remove_all)
  ;print,"x"+date_str+"x" ; debug
  ;files=file_search('/data/mta/www/mp_reports/'+date_str+'20060620/*/data/*_summ.fits')
  files=findfile('/data/mta/www/mp_reports/'+date_str+'/*/data/*_summ.fits*')
  ;files=findfile('/data/mta/www/ap_report/'+date_str+'/*/data/*_summ.fits*')
  ;files=find_file('/data/mta/www/mp_reports/20060618/*/data/*_summ.fits')
  if (n_elements(files) gt 1) then begin
    for ifiles=0, n_elements(files)-1 do begin
      ;print,files(ifiles) ; debug
      telem=mrdfits(files(ifiles),1)
      for itelem=0,n_elements(telem)-1 do begin
        if (telem(itelem).yellow gt 0.1 or telem(itelem).red gt 0.1) then begin

          name=strcompress(telem(itelem).name,/remove_all)
          b=where(maxarr(0,*) eq name,bnum)
          ;b=where(maxarr(0,*) eq "HKP27V",bnum)
          ;print, telem(itelem).name,telem(itelem).max,bnum
          if (bnum eq 1) then begin
            viol=telem(itelem).max
            lim = limits(name,describe=des,unit=units)
            lim_db=lim(n_elements(lim)-1)
            if (viol gt lim_db.yel_hi) then begin
              maxarr(day,b(0))=string(viol,format='(F7.2)')
              maxarr(8,b(0))=1
            endif
            viol=telem(itelem).min
if (name eq 'AWD3TQI') then print,files(ifiles),name,viol,lim_db.yel_lo
            if (viol lt lim_db.yel_lo) then begin
              minarr(day,b(0))=string(viol,format='(F7.2)')
              minarr(8,b(0))=1
            endif
          endif
        endif
      endfor ;for itelem=0,n_elements(telem)-1 do begin
    endfor ; for (i=0, n_elements(files)-1) do begin
  endif ; if (n_elements(files) gt 1) then begin
  date=date+1
  day=day+1
print,date," ***********************************"
  if (date eq 20110532) then date=20110601
endwhile ; while (date le tstop) do begin

;print,maxarr
;print,"x"+maxarr(0,0)+"x"
;print,maxarr(0,9)
openw,OUNIT,outfile,/get_lun
for iout=0,n_elements(maxarr(0,*))-1 do begin
  if (maxarr(8,iout) eq 1 or minarr(8,iout) eq 1) then begin
    b=where(maxarr(*,iout) eq "", bnum)
    if (bnum ge 1) then maxarr(b,iout)="x"
    b=where(minarr(*,iout) eq "", bnum)
    if (bnum ge 1) then minarr(b,iout)="x"
    viols=strjoin(minarr(1:7,iout)," | ")
    b=where(minarr(1:7,iout) ne "x",bnum)
    if (bnum gt 0) then printf, ounit, minarr(0,iout), " | ",viols
    viols=strjoin(maxarr(1:7,iout)," | ")
    b=where(maxarr(1:7,iout) ne "x",bnum)
    if (bnum gt 0) then printf, ounit, maxarr(0,iout), " | ",viols
  endif
endfor ; for iout=0,n_elements(maxarr(0,*))-1 do begin
free_lun,OUNIT
help
end