if !d.name eq 'PS' then begin
  default,ysize,12
  device,xsize=18,ysize=ysize,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/mhd/decay/fig/pns.ps
;
@parameters
default,siz,1.8
default,urms0,vA0
!x.margin=[6.3,1.3]
!y.margin=[3.2,.4]
!p.charsize=siz
;
default,pq_it1,3
default,pq_it2,30
default,pq_it3,2
it1=pq_it1 & it2=pq_it2 & it3=pq_it3
;
;  first plot
;
q=[0,1]
q2=[0,.89]
default,yr_pq_loc,[0,2]
default,pq_label,'!6'
!x.title='!8q!di!n!6(!8t!6)'
!y.title='!8p!di!n!6(!8t!6)'
plot,q,2*(1-q),xr=xr_pq_loc,yr=yr_pq_loc
bet=0. & oplot,q,q*(1.+bet),li=1 ;(beta=0 line)
bet=1. & oplot,q2,q2*(1.+bet),li=1 ;(beta=1 line)
bet=2. & oplot,q,q*(1.+bet),li=1 ;(beta=2 line)
bet=3. & oplot,q,q*(1.+bet),li=1 ;(beta=3 line)
bet=4. & oplot,q,q*(1.+bet),li=1 ;(beta=4 line)
oplot,q,q*0+.57,li=2 ;(horizontal p=1/2 line)
xyouts,.26,1.83,'!7b!6=4',siz=siz
xyouts,.49,1.85,'!7b!6=3',siz=siz*.8
xyouts,.58,1.60,'!7b!6=2',siz=siz
xyouts,.72,1.30,'!7b!6=1',siz=siz
xyouts,.78,0.92,'!7b!6=0',siz=siz
xyouts,.88,1.80,pq_label,siz=siz
!x.title='!6'
;
print,'NOTE: pq.sav must have been written by pkft_new.pro'
restore,'pq.sav'
help,qK
;
default,w,0.
nt=n_elements(qK)
for it=it1,it2<nt,it3 do begin
  psiz=.4+1.2*float(it-it1)/float(it2-it1)
  circ_sym,psiz,0
  oplot,[1,1]*qK[it],[1,1]*pK[it],col=55,ps=8
  if ihydro eq 0 then begin
    circ_sym,psiz,1
    oplot,[1,1]*qM[it],[1,1]*pM[it],col=122,ps=8
  endif
  print,it,exp(lndt(it))*urms0*k0,exp(lndt(it))*vA0*k0,qM[it],pM[it],b1M[it]
  ;print,it,Ret[it],ttau[it],qM[it],pM[it],b1M[it]
  wait,w
endfor
;
@postproc_pq_loc
;
cwd,run
;print,"$mv idl.ps ~/tex/tina/"+texdir+"/fig/pq_"+run+".ps"
print,"$mv idl.ps ~/tex/mhd/"+texdir+"/fig/pq_"+run+".ps"
END