;$Id: ppLM_log_panels.pro,v 1.6 2020/11/22 13:41:41 brandenb Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=8.8,yoffset=3 ;!p.charthick=1.1 & !p.thick=1.1 & !x.thick=1.1 & !y.thick=1.1 !p.charthick=0.9 & !p.thick=0.9 & !x.thick=0.9 & !y.thick=0.9 col0=0 endif else begin col0=255 endelse ; ; Fig.8 of the paper (Nov 2020) ; siz=.8 !p.charsize=1.4 !x.margin=[6.3,.2] !y.margin=[3.2,.4] !x.range=[0,1.] fo="(i4)" dot='!s!9!A .!n!6!r' ; xtit='!8r!6/!8r!6!d*!n' ytit=dot+'!8E!6!dK!n and '+dot+'!8E!6!dM!n' ; icount=0 !p.multi=[0,4,3] ivar1=406 & ivar2=428 & ivar3=2 for ivar=ivar1,ivar2,ivar3 do begin ; file='VAR'+str(ivar)+'.sav' restore,'sav/'+file ; iy=0 s=2e1 s=1. yr=[-.1,1]*1.7 yr=[1e-5,2] ru2=dot2(uuu)*rrr b2=dot2(bbb) ny=n_elements(yy) ; dny=5 ny1=ny/3 ny2=ny/3*2 ; if icount ge 8 then !x.title=xtit else !x.title='!6' if icount mod 4 eq 0 then !y.title=ytit else !y.title='!6' ; plot_io,xx,xx,/nodata,yr=yr ;for iy=0,ny-1,dny do oplot,xx,li=0 ,col=55 ,4*!pi*xx^2*uuu[*,iy,0]*ru2[*,iy]/2. ;for iy=0,ny-1,dny do oplot,xx,li=0 ,col=122,4*!pi*xx^2*uuu[*,iy,0]* b2[*,iy]/2.*s ; for iy=0,ny1,dny do oplot,xx,li=1 ,col=55,4*!pi*xx^2*uuu[*,iy,0]*ru2[*,iy]/2. for iy=ny1,ny2,dny do oplot,xx,li=0 ,col=55,4*!pi*xx^2*uuu[*,iy,0]*ru2[*,iy]/2. for iy=ny2,ny-1,dny do oplot,xx,li=2,col=55,4*!pi*xx^2*uuu[*,iy,0]*ru2[*,iy]/2. ; for iy=0,ny1,dny do oplot,xx,li=1 ,col=122,4*!pi*xx^2*uuu[*,iy,0]* b2[*,iy]/2.*s for iy=ny1,ny2,dny do oplot,xx,li=0 ,col=122,4*!pi*xx^2*uuu[*,iy,0]* b2[*,iy]/2.*s for iy=ny2,ny-1,dny do oplot,xx,li=2,col=122,4*!pi*xx^2*uuu[*,iy,0]* b2[*,iy]/2.*s ; if icount eq 0 then begin x0=.35 & dx=.13 & s=.7 legend,x0,dx,10^(-3.4),1,siz=siz*s,col=col0,'!7h!6 < 30!uo!n' legend,x0,dx,10^(-4.0),0,siz=siz*s,col=col0,'!630!uo!n 60!uo!n' oplot,[1,1]*.4,[1e-3,yr[1]] endif else begin oplot,[1,1]*.4,yr endelse ; xyouts,.04,.3,siz=siz,'!8t!6='+str(t,fo=fo) icount=icount+1 endfor ; !p.multi=0 ; cwd,run print,"$mv idl.ps ~/tex/jakab/wind/fig/ppLM_log_panels_"+run+".ps" END