athrz.pro

;   .r ../vis/idl/pltath
pro athrz, n
COMMON SHARE1,nx,ny,nz,nvar,nscalars
COMMON SHARE2,x,y,z
COMMON SHARE3,time,dt,gamm1,isocs
COMMON SHARE4,d,e,p,vx,vy,vz,bx,by,bz,s,phi
tag='CylNewtZ8'
tag='HKDisk'
strnum=string(n, format='(I04)')
pfact=5./3.
readvtk, tag+'.'+strnum+'.vtk', pfact

xmn=  4
xmx=  44
ymn= -0.785398
ymx=  0.785398
zmn= -15
zmx=  15
xr=xmx-xmn
yr=ymx-ymn
zr=zmx-zmn

 xx=findgen(nx)/nx*xr+xmn & yy=findgen(ny)/ny*yr+ymn  &   zz=findgen(nz)/nz*zr+zmn 
 x1=xx & x2=yy & x3=zz

dd =reform( (d[*,ny/2,*]))
dd=dd/max(dd)
vr =reform(vx[*,ny/2,*])
vp =reform(vy[*,ny/2,*])
vzi=reform(vz[*,ny/2,*])
br =reform(bx[*,ny/2,*])
bp =reform(by[*,ny/2,*])
bzi=reform(bz[*,ny/2,*])
en=reform(e[*,ny/2,*])

vp[0,0]=vp[0,0]+1e-6
pos=cglayout([4,2], xgap=7, ygap=12, oxmargin=[7,5], oymargin=[12,5])

dptr=ptrarr(12)
dptr[0]=ptr_new(dd)
dptr[1]=ptr_new(vr)
dptr[2]=ptr_new(vp)
dptr[3]=ptr_new(vzi)
dptr[4]=ptr_new(br)
dptr[5]=ptr_new(bp)
dptr[6]=ptr_new(bzi)
dptr[7]=ptr_new(en)

sarr=strarr(12)
sarr[0]="density"
sarr[1]="V!DR!N"
sarr[2]="V!D!7u!X!N"
sarr[3]="V!DZ!N"
sarr[4]="B!DR!N"
sarr[5]="B!D!7u!X!N"
sarr[6]="B!DZ!N"
sarr[7]="E!DT!N"

cgloadct,33
;tvlct,255,255,255,0
;tvlct,255,255,255,1
cgerase,0

;cgdisplay, xs=800,ys=800
for i=0,7 do begin
d=*dptr[i]
imin=min(d)
imax=max(d)
r=cgscalevector(d,1,254)
px=pos[*,i]
if ( i eq 0 ) then begin
cgimage, r, pos=px
endif else begin
cgimage, r, pos=px, /noerase
endelse 
cgcontour,r,xx,zz,  $ 
    /noerase, /nodata, $
    pos=px,  $
    title=sarr[i], $
    xtitle="R",$
    ytitle="Z"
    cgcolorbar, range=[imin,imax], pos=[px[0],px[1]-0.09, px[2], px[1]-0.08]

endfor

cgtext, 0.4,0.97, "athena GT4 turbulent accretion disk , t="+string(time/61.5)+" orbits", /normal

fname=tag+'-rz-'+strnum
print, fname
im=cgsnapshot(filename=fname, /jpeg, /nodialog)


 end