--------------------------------------------------------------------------
--------------------------------------------------------------------------
SATIRE-T_wu18_tsi.txt

----------
HEADER:

Solar Total Irradiance, 1642 ~ 31.05.2017

File structure:
1. JD: array [112349], julian day
  1642 ~ 1709, yearly cadence (given at July 2)
  1.1.1710 ~ 5.31.2017, daily cadence
2. TSI(JD): [112349], W/m2
Reference: Wu, C.-J., Krivova, N. A., Solanki, S. K., Usoskin, I. G., A&A, 620, A120, 12

----------
IDL code to read the file:

filename='PATH_TO_FILE/SATIRE-T_wu18_tsi.txt'
temp=read_ascii(filename,data_start=9,header=header,comment_sym=';')
temp=temp.(0)
date=temp[0,*]
tsi=temp[1,*]

;time in fraction of year
caldat,date,m,d,y
t2=double(julday(m,d,y)-julday(1,1,y))
t1=double(julday(12,31,y)-julday(1,1,y)+1)
time=t2/t1 + y
--------------------------------------------------------------------------
--------------------------------------------------------------------------

--------------------------------------------------------------------------
--------------------------------------------------------------------------
SATIRE-T_wu18_ssi.txt

----------
HEADER:

Solar Spectral Irradiance reconstruction using the SATIRE-T model, 1642 ~ 5.31.2017
 
File structure
1. wavelength: array[1070],nm
2. wavelength bin: array[1070],nm
3. date: array[112349], julian day
   1642 ~ 1709, yearly cadence (given at July 2)
   1.1.1710 ~ 5.31.2017, daily cadence
4. SSI(wavelength,date): array[1070,112349], W m-2 nm-1
Reference: Wu, C.-J., Krivova, N. A., Solanki, S. K., Usoskin, I. G., A&A, 620, A120, 12

----------
IDL code to read the file:

filename='PATH_TO_FILE/SATIRE-T_wu18_ssi.txt'
N=112349
header=strarr(11)
satire_wl=fltarr(1070)
satire_dwl=fltarr(1070)
date=lonarr(N)
ssi=dblarr(1070,N)

openr,1,filename
readf,1,header
readf,1,satire_wl
readf,1,satire_dwl
readf,1,date
readf,1,ssi
close,1

;Calculate TSI
N=112349
tsi=dblarr(N)
for i=0L,N-1 do tsi[i]=total(satire_dwl*ssi[*,i])

;time in fraction of year
caldat,date,m,d,y
t2=double(julday(m,d,y)-julday(1,1,y))
t1=double(julday(12,31,y)-julday(1,1,y)+1)
time=t2/t1 + y
--------------------------------------------------------------------------
--------------------------------------------------------------------------