; This is to read files created by the GSE for MAY the SEP attenuator life test

pro read_attenuator_test_file,file =  file,  $
                              durationa =  durationa,  durationb =  durationb
  ;file ='C:\Users\rlillis\Work\MAVEN\SEP_instrument\GSE\actTest_test_2.dat'
  if not keyword_set (file) then stop
  
string1 = 'gh'
String2 = string1
openr,unit,file,/get_lun
Date = 'na'
Time = 'na'
Pulse_duration =  sqrt (-9.3)
side='na'; actuating side A (B)means pulling the attenuator out (in)
sense_A =0b; 0 means switch A is closed, i.e. attenuator is in, 1 means switch A is open, or attenuator out
sense_B =0b; 0 means switch B is closed, i.e. attenuator is out, 1 means switch B is open, or attenuator in
Count = 0L
While (eof(unit) eq 0) do begin
readf, unit, string1 

tmp1 = strsplit (string1, ' ',/extract)
print, tmp1[1]
date = [date, tmp1 [0]]
Tmp2 = strsplit (tmp1 [1],',',/extract)
Time = [time, tmp2 [0]]
sense_a = [sense_a,tmp2 [1]]
sense_b = [sense_b, tmp2 [2]]
readf, unit, string1
tmp1 = strsplit (string1, ' ',/extract)
print, tmp1[1]
tmp2 = strsplit (tmp1 [1], ',',/extract)

pulse_duration = [pulse_duration, float (tmp2[2])]
Side = [side,'A']

if eof(unit) then break
readf, unit, string1 
tmp1 = strsplit (string1, ' ',/extract)
print, tmp1[1]

date = [date, tmp1 [0]]
Tmp2 = strsplit (tmp1 [1],',',/extract)
Time = [time, tmp2 [0]]
sense_a = [sense_a,tmp2 [1]]
sense_b = [sense_b, tmp2 [2]]

if eof(unit) then break
readf, unit, string1 & $
tmp1 = strsplit (string1, ' ',/extract )
print, tmp1[1]
tmp2 = strsplit ( tmp1 [1], ',',/extract)
pulse_duration = [pulse_duration, float (tmp2[2])]
Side = [side,'B']
Count = count +1
endwhile
aside = where (side eq 'A')
bside =  where (side eq 'B')
durationa=pulse_duration (aside)
Durationb = pulse_duration (bside)
print, 'Actuator A:',mean (durationa), stddev (durationa)
Print, 'actuator B:', Mean (durationb), stddev (durationb)
time_string =  date [1:*] + '/' + time[1:*]
time_double =  time_double (time_string)
minutes =  (time_double - time_double [0])/60.0


Plot, minutes[aside], Durationa,  $
  xtit =  'minutes',  ytit =  'stroke duration, ms',  $
  yr =  minmax ([durationa,  durationb]),  $
  title =  'Black: Open, Blue: Close'
Oplot, minutes[bside], durationb, Color = 54
xyouts,  100,  350,  'Open'
xyouts,  100,  305,  'Close',  Color =  54
Plot, Durationa,  $
  xtit =  'actuations',  ytit =  'stroke duration, ms',  $
  yr =  minmax ([durationa,  durationb])
Oplot, durationb, Color = 54
;stop
end