;$Id: ppower_all2.pro,v 1.2 2026/06/30 12:35:16 brandenb Exp $
col0=255
@parameters
print,"$sed.csh data/param2.nml"
pc_read_param,/param2,obj=param2
pc_read_param,obj=param
pc_read_ts,obj=ts
;
power,'o','_kin',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot,/lks
power,'_lr','_ro',k=k,spec1=spec3,spec2=spec4,i=n,tt=t,/noplot,/lks
dk=k[1]-k[0]
;
;  This routine is used to view all power spectra
;  that were written during the run (controlled by dspec)
;  The plot range yr is set to a default value.
;  The time interval of plotting is "w" (default value)
;
;   5-oct-02/axel: written
;
spawn,'if (! -f param.pro) touch parameters.pro'
nt=n_elements(t)
xr=[k[1],max(k)]
k1=k & k1[0]=1. & k1=1./k1 & k1[0]=0.
EMag=fltarr(nt)
xiMag=fltarr(nt)
Lu=fltarr(nt)
IA=fltarr(nt)
tA=fltarr(nt)
;
;  This routine also computes Cxi and CEE.
;  We take the value at t=1e5.
;
default,w,.01
default,i1,0
default,i2,nt-1
i2=nt-1
default,istride,1
;
close,1
openw,1,'kpm.txt'
printf,1,"   i     t        EMag        xiMag         CM          Lu          IA          Cxi         CEE"
fo="(i4,e10.2,7e12.4)"
for i=i1,i2,istride do begin
  plot_oo,k,spec1(*,i),xr=xr,yr=yr,tit='t='+str(t(i))
  oplot,k,spec2(*,i)*k^2,col=55
  oplot,k,spec3(*,i),col=155
  ;
  wait,w
endfor
close,1
;
spawn,'cvs add kpm.txt'
;save,file='grav1_last.sav',k,grav1_last,grav1_next
END