PRO STEREO_Read, sta, stb

; Program to read the STEREO A and B ASCII magnetic files downloaded directly from the
; UCLA site.

data_template = {ST_Ascii, version:1.0, datastart:0L, delimiter:0B, missingvalue:0.0, $
                  commentsymbol:'', fieldcount:0L, fieldtypes:fltarr(11), fieldnames:strarr(11), $
                 fieldlocations:lonarr(11), fieldgroups:lonarr(11)}


STA_DataFile = 'D:\Data\STEREO\Mag\STA_L1_MAG_20070810-20070910.dat'
STB_DataFile = 'D:\Data\STEREO\Mag\STB_L1_MAG_20070810-20070910.dat'

data_template.version = 1.0
data_template.datastart = 94
data_template.delimiter = 32
data_template.missingvalue = 1.e+32
data_template.commentsymbol = ''
data_template.fieldcount = 11
data_template.fieldtypes = [3,3,3,3,3,3,3,4,4,4,4]
data_template.fieldnames = ['Year','Month','Day','Hr','Mn','Sec','mSec','Br','Bt','Bn','B']
data_template.fieldlocations = [0,5,9,12,15,18,21,27,39,52,66]
data_template.fieldgroups = [0,1,2,3,4,5,6,7,8,9,10]

; Define common STEREO Mag data structure

mag_data = {Mag, Yr:0, Month:0, Day:0, DOY:0.0d0, MJD2000:0.0d0, Hr:0, Mn:0, Sec:0, mSec:0, $
            Br:0.0d0, Bt:0.0d0, Bn:0.0d0, B:0.0d0}

; Read STEREO A data

sta_data = read_ascii(STA_DataFile, template=data_template)

sta = replicate({Mag},n_elements(sta_data.year))
sta.yr = sta_data.year
sta.month = sta_data.month
sta.day = sta_data.day
sta.hr = sta_data.hr
sta.mn = sta_data.mn
sta.sec = sta_data.sec
sta.msec = sta_data.msec
sta.mjd2000 = julday(sta.month,sta.day,sta.yr,sta.hr,sta.mn,sta.sec+sta.msec/1000.d0)-2451544.5d0
doy_md, sta.yr, doy, sta.month, sta.day, /reverse
sta.doy = ((((sta.msec/1000.d0)+sta.sec)/60.d0+sta.mn)/60.d0+sta.hr)/24.0+doy
sta.br = sta_data.br
sta.bt = sta_data.bt
sta.bn = sta_data.bn
sta.b = sta_data.b

; Read STEREO B data

stb_data = read_ascii(STB_DataFile, template=data_template)

stb = replicate({Mag},n_elements(stb_data.year))
stb.yr = stb_data.year
stb.month = stb_data.month
stb.day = stb_data.day
stb.hr = stb_data.hr
stb.mn = stb_data.mn
stb.sec = stb_data.sec
stb.msec = stb_data.msec
stb.mjd2000 = julday(stb.month,stb.day,stb.yr,stb.hr,stb.mn,stb.sec+stb.msec/1000.d0)-2451544.5d0
doy_md, stb.yr, doy, stb.month, stb.day, /reverse
stb.doy = ((((stb.msec/1000.d0)+stb.sec)/60.d0+stb.mn)/60.d0+stb.hr)/24.0+doy
stb.br = stb_data.br
stb.bt = stb_data.bt
stb.bn = stb_data.bn
stb.b = stb_data.b

end