;$Id: pcomp_1d_initial.pro,v 1.2 2020/08/14 07:55:30 brandenb Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=10,yoffset=3 !p.charthick=2.2 & !p.thick=2.2 & !x.thick=2.2 & !y.thick=2.2 end ; siz2=1. !p.charsize=1. !x.margin=[7.8,.5] !y.margin=[3.2,.5] ; dirA='~/pencil-code/christer/isoth_expwave/lagrangian/1d/256a1/data' dirB='~/pencil-code/christer/isoth_expwave/eulerian/1d/256e/data' !p.multi=[0,2,2] ; da=1e-6 a1=1e-4 & a2=3e-4 & a3=1e-3 & a4=3e-3 & a5=1e-2 & a6=3e-2 & a7=1e-1 xcrit=1.5 xout=-1.3 yout1=.50 yout2=3.4 ; ;----------------------------------------------------------------------------- ivar=1 pc_read_var,obj=varB,/trimall,ivar=ivar,datadir=dirB xr=minmax(varB.x) ; yr=[-1,1]*.7 !x.title='!6' !y.title='!8u!6' plot,xr,yr,ps=3,/nodata oplot,varB.x,varB.uu[*,0] xyouts,xout,yout1,siz=siz2,'!6(!8a!6)' ; yr=[.95,3.8] !y.title='!7q!6' plot,xr,yr,ps=3,/nodata oplot,varB.x,exp(varB.lnrho) xyouts,xout,yout2,siz=siz2,'!6(!8b!6)' ; ;----------------------------------------------------------------------------- ivar=5 pc_read_var,obj=varB,/trimall,ivar=ivar,datadir=dirB xr=minmax(varB.x) ; yr=[-1,1]*.7 !x.title='!8x!6' !y.title='!8u!6' plot,xr,yr,ps=3,/nodata oplot,varB.x,varB.uu[*,0] xyouts,xout,yout1,siz=siz2,'!6(!8c!6)' ; yr=[.95,2.1] yr=[.95,3.8] !y.title='!7q!6' plot,xr,yr,ps=3,/nodata oplot,varB.x,exp(varB.lnrho) xyouts,xout,yout2,siz=siz2,'!6(!8d!6)' ; print,"$mv idl.ps ~/tex/christer/isothermal_expwave_git/fig/pcomp_1d_initial.eps" END