pro st_thruster_load,trange=trange,probe=probe,tshift=tshift


if not keyword_set(source_options) then begin
    stereo_init
    source_options = !stereo
endif
mystereo = source_options

vb =mystereo.verbose

if not keyword_set(probe) then probe =  'a'

fileres = 3600l*24     ; one day resolution in the files

tr = timerange(trange)
n = ceil((tr[1]-tr[0])/fileres)  > 1
dates = dindgen(n)*fileres + tr[0]

pn = (byte(probe)-byte('a'))[0]

path = '.thruster/SW_'+strupcase(probe)+'_Thrusts_YYYYMMDD.txt'

relpathnames= time_string(dates,tformat= path)
files = file_retrieve(relpathnames,_extra = mystereo)

format = {dt:0d,  t1:0d,  t2:0d,  ur8:0d, bits:[0b,0b,0b,0b]}
data = read_asc(files[0],format=format,nheader=5)

t = (data.ur8 + 365*12 + 3 )* 24 * 3600
data.bits +=  [7,5,3,1] # replicate(1,n_elements(t))

;printdat,data

if not keyword_set(tshift) then tshift =0d
store_data,'st'+probe+'_thruster',data = {x:t+tshift, y:transpose(data.bits) },dlimit = {colors:'rbgm', labels:['A1','A2','A3','A4']}


end