;$Id: pspecm_hel_delk_k1_kf100_combined.pro,v 1.2 2021/04/01 16:50:31 brandenb Exp $ if !d.name eq 'PS' then begin ;device,xsize=18,ysize=8,yoffset=3 device,xsize=18,ysize=6,yoffset=3 !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6 end ; xfac=1./6. yfac=1./2.8 xfac2=1./5.0 yfac2=1.8 ; dir1='M1152e_exp6k4_k1_kf100_delk03' dir2='M1152e_exp6k4_k1_kf100_delk1' dir3='M1152e_exp6k4_k1_kf100_delk3' dir4='M1152e_exp6k4_k1_kf100_delk10' ; dir1b='M1152e_exp6k4' dir2b='M1152e_exp6k4_delk10' dir3b='M1152e_exp6k4_delk50' dir4b='M1152e_exp6k4_delk200' ; xs=.34 ys1=.06 ys2=.6 ; xr=[2e-9,3.8e-5] yr_pspecm_scl3=[1d-18,1.3d-13] yr_pspechm_scl3=[2d-23,8.2d-18] ; thick4=5 siz=1.0 si2=1.0 !p.charsize=1.2 !x.margin=[7.8,0.5] !y.margin=[3.2,0.2] !p.multi=[0,2,1] ;fac_h=1.263e-18 fac2=2. half='!s!u 1!n!r!s-!r!d 2!n' specm='/specm.sav' wspecm='/specm.txt' label='' ; ; determine normalization ; default,maglabel,'mag' default,green_lines,0 default,plot_polarization,1 ytickf='logticks_exp' yr=yr_pspecm_scl3 ; ; new units ; T_in100GeV=1. & gstar_in100=1. T_in100GeV=1.5d-3 & gstar_in100=.15 ;gS_in100=1. ;gS_in100=.0391 tend=1. Hstar=2.066d10*T_in100GeV^2*gstar_in100^.5 H0=3.241d-18 ;a0=1.254d15*T_in100GeV*gS_in100^(1./3.) a0=1d12 fact=(Hstar/H0)^2*(tend/a0)^4 print,'fact=',fact ; ;Hstar=2.066d10 ;a0=1.2545d15 freq=Hstar/(2.*!pi*a0) H0=3.24d-18 ;changed from H0 = 100 km/s/Mpc to 70 km/s/Mpc ;H0=2.268545d-18 rhocrit_fac=3./(8.*!pi) ;fac_h=1.263e-18 ;changed H0 from 100 to 70 km/s/Mpc, and computed h_c ;from formula fac_h=sqrt(3./2.)*H0/!pi fac=(Hstar/H0)^2/rhocrit_fac*(1./a0)^4 facGW=1./(16.*!pi) circ_sym,0.7,1 mixed=0. ; ; 1st panel ; !p.title='!6';+run !x.title='!8f!6 [Hz]' !y.title='!6!8h!6!s!d0!n!r!u2!n!7X!6!dGW!n(!8f!6)' ; ;----------------------------------------------------------------------------- restore,'../'+dir1+specm grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot plot_oo,xr,yr,ytickf=ytickf,/nodata oplot,freq*k,OmegaGW ; restore,'../'+dir2+specm grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,freq*k,OmegaGW,col=55 ; loadct,6 restore,'../'+dir3+specm grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,freq*k,OmegaGW,col=122 loadct,5 ; restore,'../'+dir4+specm grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,freq*k,OmegaGW,col=155 ; ;----------------------------------------------------------------------------- restore,'../'+dir1b+specm grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,xfac*freq*k,yfac*OmegaGW ; restore,'../'+dir2b+specm grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,xfac*freq*k,yfac*OmegaGW,col=55 ; loadct,6 restore,'../'+dir3b+specm grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,xfac*freq*k,yfac*OmegaGW,col=122 loadct,5 ; restore,'../'+dir4b+specm grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,xfac*freq*k,yfac*OmegaGW,col=155 ; ; fits ; xx=grange(1e-9,3e-7,41,/log) oplot,xx,xx*1e-7 oplot,thick=3,xx,7e-8*(xx)^3/(2.3e-8^2+xx^2),li=2 oplot,thick=3,xx,9e-8*(xx)^3/(6.0e-9^2+xx^2),li=2 ; xout=1.4e-5 xyouts,xout,5.0e-14,'!6(a)',siz=si2 xx=[1.2e-6,1.5e-5] & oplot,xx,8.0d-30/xx^(8./3.) xyouts,3e-6,7e-15,'!9A!8k!6!u-8/3!n',siz=siz ; ; 2nd panel ; !y.title='!8h!6!dc!n(!8f!6)' default,kin,0 default,hdone,1 ytickf='logticks_exp' yr=yr_pspechm_scl3 circ_sym,0.7,1 ; ;----------------------------------------------------------------------------- restore,'../'+dir1+specm fach=1.263e-18/(freq*k[1:*]) grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot plot_oo,freq*k,sqrt(k*grah1m)/a0,xr=xr,yr=yr,ytickf=ytickf,/nodata oplot,freq*k,sqrt(k*grah1m)/a0 ; ; 2nd line ; restore,'../'+dir2+specm fach=1.263e-18/(freq*k[1:*]) grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,freq*k,sqrt(k*grah1m)/a0,col=55 ; ; 3rd line ; loadct,6 restore,'../'+dir3+specm fach=1.263e-18/(freq*k[1:*]) grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,freq*k,sqrt(k*grah1m)/a0,col=122 loadct,5 ; ; 4th line ; restore,'../'+dir4+specm fach=1.263e-18/(freq*k[1:*]) grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,freq*k,sqrt(k*grah1m)/a0,col=155 ; ;----------------------------------------------------------------------------- restore,'../'+dir1b+specm fach=1.263e-18/(freq*k[1:*]) grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,xfac2*freq*k,yfac2*sqrt(k*grah1m)/a0 ; ; 2nd line ; restore,'../'+dir2b+specm fach=1.263e-18/(freq*k[1:*]) grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,xfac2*freq*k,yfac2*sqrt(k*grah1m)/a0,col=55 ; ; 3rd line ; loadct,6 restore,'../'+dir3b+specm fach=1.263e-18/(freq*k[1:*]) grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,xfac2*freq*k,yfac2*sqrt(k*grah1m)/a0,col=122 loadct,5 ; ; 4th line ; restore,'../'+dir4b+specm fach=1.263e-18/(freq*k[1:*]) grav_tot=grav1m+grah1m/tm[0]^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,xfac2*freq*k,yfac2*sqrt(k*grah1m)/a0,col=155 ; xyouts,xout,2.6e-18,'!6(b)',siz=si2 xx=[1e-6,1.8e-5] & oplot,xx,5d-33/xx^(7./3.) xyouts,3e-6,6e-20,'!9A!8k!6!u-7/3!n',siz=siz ; print,"$mv idl.ps ~/GitHub/Yutong/MGW-NANOGrav/Figures/pspecm_hel_delk_k1_kf100_combined.eps" END