PRO ts_mean,obj,hydro=hydro,dissipation=dissipation
;
; Calculates mean value of saturated urms and brms values (t1<tt<t2).
; If mean_val_lin='true' the kinematic urms value is calculated 
; (t1_lin<tt<t2_lin).
;
lmhd=0
lhydro=0
ldiss=0
if (keyword_set(hydro)) then lhydro=-1
if (not keyword_set(hydro)) then lmhd=-1
if (keyword_set(dissipation)) then ldiss=-1
;print,'lhydro=',lhydro
;print,'lmhd=',lmhd

@param
default,t1,0.
default,t2,1000.
default,nsmooth,201
default,mean_val_lin,'false'
;
tt=obj.t
;
; Find urms or dissipation
;
if (ldiss) then begin
    kinetic=obj.epsK
    kinout='epsK='
endif else begin
    kinetic=obj.urms
    kinout='mean_urms='
end
;
; MHD or only hydro
;
if (lmhd) then begin
    ;
    ; Find brms or dissipation
    ;
    if (ldiss) then begin
        magnetic=obj.epsM
        magout='epsM='
    endif else begin
        magnetic=obj.brms
        magout='mean_brms='
    end
end
;
;  Kinematic or non linear phase
;
if mean_val_lin eq 'true' then begin
	good=where(tt ge t1_lin and tt le t2_lin)
endif else begin
	good=where(tt ge t1 and tt le t2)
endelse
;
;  Begin calculation and plotting
;
if mean_val_lin eq 'true' then begin
	mean_kinetic=mean(kinetic(good))
	print,kinout,mean_kinetic
	!p.multi=[0,1,1]
	plot,tt,kinetic
	ttt=[t1_lin,t2_lin] & oplot,ttt,ttt-ttt+mean_kinetic,col=122
endif else begin
    if (lmhd) then begin
        mean_kinetic=mean(kinetic(good))
	print,kinout,mean_kinetic
	!p.multi=[0,1,2]
	plot,tt,kinetic
	ttt=[t1,t2] & oplot,ttt,ttt-ttt+mean_kinetic,col=122
	;
	mean_magnetic=mean(magnetic(good))
	print,magout,mean_magnetic
	plot_io,tt,magnetic
	ttt=[t1,t2] & oplot,ttt,ttt-ttt+mean_magnetic,col=122
    endif else begin
        mean_kinetic=mean(kinetic(good))
	print,kinout,mean_kinetic
	!p.multi=[0,1,1]
	plot,tt,kinetic
        ttt=[t1,t2] & oplot,ttt,ttt-ttt+mean_kinetic,col=122 
    end
endelse
;
!p.multi=[0,1,1]
;
END