;path to the model
ptk='./'

;model index
nmod='088200'

close,1
openr,1,ptk+'ltestokes.'+nmod+'.0000.0000',/f77_unformatted

;variable definitions
alpha=0.0
phi=0.0
xi=0.0
nsx=0
nsz=0
nwl=0

;viewing angles, set to 0, since solar disk centre
readu,1,alpha
readu,1,phi
readu,1,xi

;horizontal box sizes (288x288)
readu,1,nsx
readu,1,nsz

;number of wavelength points
readu,1,nwl

;array definitions

;wavelength axis
lambda=fltarr(nwl)

;continuum intensity
icont=fltarr(nsx,nsz)

;stokes profiles I, V, Q, U
stokes=fltarr(4,nwl,nsx,nsz)


readu,1,lambda
readu,1,icont
readu,1,stokes

close,1



end