;$Id: power_frac.pro,v 1.2 2020/08/17 19:04:21 brandenb Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=14,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; !p.charsize=1.7 !x.margin=[7.8,.5] !y.margin=[3.2,.5] !p.multi=[0,1,2] @parameters ; power,'_mag','hel_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot,/lks yr_power_frac=[3e-9,3e-3] & yr_frac=[-.3,1.05] ; it1=-1 it3=10 nt=n_elements(t) it2=nt-1 col=grange(55,200,nt) ; !x.title='!6' !y.title='!8E!6!dM!n(!8k!6)' plot_oo,k[1:*],spec1[1:*,it2],yr=yr_power_frac,/nodata loadct,6 for it=it1,it2,it3 do oplot,k[1:*],spec1[1:*,it > 0],col=fix(col(it > 0)) loadct,5 ; !x.title='!8k!6' !y.title='!8kH!6!dM!n(!8k!6)/2!8E!6!dM!n(!8k!6)' plot_oi,k[1:*],k[1:*]*0,yr=yr_frac loadct,6 for it=it1,it2,it3 do oplot,k[1:*],.5*k[1:*]*spec2[1:*,it > 0]/spec1[1:*,it > 0],col=col(it > 0) for it=it1,it2,it3 do print,'t=',t[it > 0] loadct,5 ; cwd,run print,"$mv idl.ps ~/tex/notes/GW_decay_compens/fig/power_frac_"+run+".ps" END