pro plt1, array
sd=stdev(array,mean)
print, 'sd=', sd
print, 'mean=' ,mean
openw, 1, 'in.dat'
sz=size(array)
print, 'sz=', sz
count=sz(sz(0)+2)
print, 'count=', count
printf, 1, count
printf, 1, array
printf, 1, sd
printf, 1, mean
close, 1
openr, 1, 'in.dat'
arr_size=1
readf, 1, arr_size
new_array=fltarr(arr_size)
readf, 1, new_array
print, 'array read back'
print, new_array
readf, 1, a
readf, 1, b
print, 'mean=', a
print, 'sd=', b
end