; This is to read files created by the GSE for MAY the SEP attenuator life test pro mav_sep_atten_test_load source = mav_file_source() pathname = 'maven/sep/test/atten/actTest_20110216.dat' file = file_retrieve(pathname,_extra= source) 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) 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) tmp2 = strsplit (tmp1 [1], ',',/extract) pulse_duration = [pulse_duration, float (tmp2[2])] Side = [side,'A'] readf, unit, string1 tmp1 = strsplit (string1, ' ',/extract) 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 ) tmp2 = strsplit ( tmp1 [1], ',',/extract) pulse_duration = [pulse_duration, float (tmp2[2])] Side = [side,'B'] Count = count +1 endwhile durationa=pulse_duration (where (side eq 'A')) Durationb = pulse_duration (where (side eq 'B')) print, 'Actuator A:',mean (durationa), stddev (durationa) Print, 'actuator B:', Mean (durationb), stddev (durationb) Plot, Durationa Oplot, durationb, Color = 54 stop end