power,'_kin','hel_kin',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot ; ; $Id: ppower_all.pro,v 1.2 2019/02/25 13:53:42 brandenb Exp $ ; ; 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' @parameters xr=[1,max(k)] istride=1 default,w,.05 default,t1,.0 default,iover,0 default,istride,1 default,istride,1 default,helplot,1 default,yr,[3e-12,3e-2] ampl15=3e-10 & ampl4=1e-13 ampl15=3e-6 & ampl4=1e-8 ; good=where(t ge t1) i1=min(good) i2=n-2 ; circ_sym,0.8,1 for i=i1,i2,istride do begin if iover eq 0 or i eq 0 then begin plot_oo,k,spec1(*,i),xr=xr,yr=yr,tit='t='+str(t(i)),/nodata endif oplot,k,abs(spec1(*,i)) if helplot then begin oplot,k,abs(spec2(*,i)),li=1 oplot,k,+spec2(*,i),ps=8,col=122 oplot,k,-spec2(*,i),ps=8,col=55 endif kk=[1,20] & oplot,kk,.01/kk^1.6667;,col=188 ;kk=[1,20] & oplot,kk,.001*ampl15*kk^2;,col=188 ;kk=[1,5] & oplot,kk,ampl4*kk^4,col=55 print,i,t(i) wait,w end ; END