pc_read_ts,obj=ts
pc_read_param,obj=param,/param2
;
spawn,'touch parameters.pro'
@parameters.pro
;
default,t1,0.
default,t2,9e29
default,t1b,0.
default,t2b,9e29
;
;  Note: put a minus sign here, because of the minus sign
;  in the definition of kappa.
;
tt=ts.t
ucm=ts.kap33
lam=param.lam_testscalar
om=param.om_testscalar
c=cos(om*tt)
s=sin(om*tt)
;
good=where(tt gt t1 and tt le t2)
tgood=tt(good)
good2=where(tgood gt t1b and tgood le t2b)
;
!p.multi=[0,1,2]
;
ucmc=ucm*c
;ucmcm=accum(ucmc(good))
ucmcm=mean(ucmc(good))
kappa_re=mean(ucmcm(good2))
plot,tt,ucmc,yst=3
;oplot,tgood,ucmcm,col=122
;oplot,tgood(good2),accum(ucmcm(good2)),col=188
oplot,tgood(good2),tgood(good2)*0+ucmcm,col=188
;
ucms=ucm*s
;ucmsm=accum(ucms(good))
ucmsm=mean(ucms(good))
kappa_im=mean(ucmsm(good2))
plot,tt,ucms
;oplot,tgood,ucmsm,col=122
;oplot,tgood(good2),accum(ucmsm(good2)),col=188
oplot,tgood(good2),tgood(good2)*0+ucmsm,col=188
;
cwd,run
default,file,'to16.dat'
openw,1,file
fo='(2f7.3,2e10.2,2x,a)'
print,lam,om,kappa_re,kappa_im,run,fo=fo
printf,1,lam,om,kappa_re,kappa_im,run,fo=fo
close,1
;
spawn,'cat '+file+'>>../idl/'+file
END