function write_gauss, intervall, rday, intervall_a4000d1, intervall_a4000d2, intervall_a4000c2, intervall_a4000s2, intervall_p4000d1, intervall_p4000d2, intervall_p4000c2, intervall_p4000s2
;--------------nimmt die über das intervall angesammelten werte und
;berechnet daraus den gaussfit

intervall_gauss='/data/etnas1/ulysses/ket/rate/shifttab/gauss-'+strtrim(string(intervall),1)+'.dat'
openw, 5, intervall_gauss, /append, width=400

if (n_elements(intervall_a4000d1) gt 1 and n_elements(intervall_a4000d2) gt 1 and n_elements(intervall_a4000c2) gt 1  and n_elements(intervall_a4000s2) gt 1 and n_elements(intervall_p4000d1) gt 1 and n_elements(intervall_p4000d2) gt 1 and n_elements(intervall_p4000c2) gt 1  and n_elements(intervall_p4000s2) gt 1) then begin

;---------------baue histogramm

printf, 5, rday, mean(intervall_a4000d1[1:*]), $
stddev(intervall_a4000d1[1:*]),  $
mean(intervall_a4000d2[1:*]), $
stddev(intervall_a4000d2[1:*]), $
mean(intervall_a4000c2[1:*]), $
stddev(intervall_a4000c2[1:*]),  $
mean(intervall_a4000s2[1:*]), $
stddev(intervall_a4000s2[1:*]), $
mean(intervall_p4000d1[1:*]), stddev(intervall_p4000d1[1:*]), $
mean(intervall_p4000d2[1:*]), stddev(intervall_p4000d2[1:*]), $ 
mean(intervall_p4000c2[1:*]), stddev(intervall_p4000c2[1:*]), $
mean(intervall_p4000s2[1:*]), stddev(intervall_p4000s2[1:*])
endif
close, 5
free_lun, 5
end