;$Id: pcomp_spec_christer_Re.pro,v 1.10 2021/05/01 13:27:02 brandenb Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=12,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; !p.charsize=1.7 !x.margin=[9.4,.5] !y.margin=[3.2,.5] ; ; sample to compare time series form different directories ; run='hel' yr=[1e-11,1e-1] yr=[1e-4,2e-2] ; V2 isothermal_helical_256_f4.0_dtf1.0_5adjagr_tsc_nu0.01_noshocknu ; V3 isothermal_helical_256_f4.0_dtf1.0_5adjagr_tsc_nu0.02_noshocknu ; V4 isothermal_helical_256_f4.0_dtf1.0_5adjagr_tsc_nu0.05_noshocknu christer_dir1='isothermal_helical_512_f4.0_dtf1.0_5adjagr_tsc_nu0.01_noshocknu_512proc' ;(used in other Fig) christer_dir2='isothermal_helical_256_f4.0_dtf1.0_5adjagr_tsc_nu0.02_noshocknu' christer_dir3='isothermal_helical_512_f4.0_dtf1.0_5adjagr_tsc_nu0.005_noshocknu_512proc' ;(ran until t=128) ; xr=[1.,256.] fact=.2 plot_oo,xr,yr,xtit='!8k!6/!8k!6!d1!n',ytit='!8P!dn!n!6(!8k!6)',/nodata ; xx=1.2 & xxx=.8 & siz=1.1 & ldy=.17 legend,siz=siz,xx*1.4,0.,10.^(-3.9+5*ldy),0,'!8St!6!dint!n=' legend,siz=siz,xx,xxx,10.^(-3.9),0,col=55,'!631' loadct,6 legend,siz=siz,xx,xxx,10.^(-3.9+1*ldy),0,col=122,'!63.1' loadct,5 legend,siz=siz,xx,xxx,10.^(-3.9+2*ldy),0,'!63.1!9X!610!u-1!n',col=122 legend,siz=siz,xx,xxx,10.^(-3.9+3*ldy),0,'!63.1!9X!610!u-2!n',col=155 legend,siz=siz,xx,xxx,10.^(-3.9+4*ldy),0,'!63.1!9X!610!u-3!n' ; xx=45. & xxx=30. & siz=1.1 legend,siz=siz,xx,xxx,10.^(-3.9+3*ldy),1,'!8Re!6=20, V2a' legend,siz=siz,xx,xxx,10.^(-3.9+2*ldy),2,'!8Re!6=45, V2' legend,siz=siz,xx,xxx,10.^(-3.9+1*ldy),0,'!8Re!6=95, V2b' ; fact2=5. fact2=5. fact2=5. fact2=5. fact2=4.*64. restore,'../../'+christer_dir1+'/mean_spec.sav' oplot,k[1:*],fact2*mspec1[1:*],li=2 oplot,k[1:*],fact2*mspec2[1:*],li=2,col=155 oplot,k[1:*],fact2*mspec3[1:*],li=2,col=122 loadct,6 oplot,k[1:*],fact2*mspec4[1:*],li=2,col=122 loadct,5 oplot,k[1:*],fact2*mspec5[1:*],li=2,col=55 ; fact2=4. restore,'../../'+christer_dir2+'/mean_spec.sav' oplot,k[1:*],fact2*mspec1[1:*],li=1 oplot,k[1:*],fact2*mspec2[1:*],col=155,li=1 oplot,k[1:*],fact2*mspec3[1:*],col=122,li=1 loadct,6 oplot,k[1:*],fact2*mspec4[1:*],col=122,li=1 loadct,5 oplot,k[1:*],fact2*mspec5[1:*],col=55,li=1 ; thick1=3. fact2=64.*4. restore,'../../'+christer_dir3+'/mean_spec.sav' oplot,k[1:*],fact2*mspec1[1:*],li=0 oplot,k[1:*],fact2*mspec2[1:*],col=155,li=0 oplot,k[1:*],fact2*mspec3[1:*],col=122,li=0 loadct,6 oplot,k[1:*],fact2*mspec4[1:*],col=122,li=0 loadct,5 oplot,k[1:*],fact2*mspec5[1:*],col=55,li=0 ; xyouts,180,.0140,'!6(!8a!6)',siz=siz ; print,"$mv idl.ps ~/tex/christer/isothermal_expwave_git/fig/pcomp_spec_christer_Re_"+run+".eps" ;print,"$mv idl.ps ~/tex/christer/isothermal_expwave_extra/fig/pcomp_spec_christer_Re_"+run+".eps" END