if !d.name eq 'PS' then begin
  ysize=10
  device,xsize=18,ysize=ysize,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
@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,0
default,pq_it2,1e9
default,pq_it3,1
it1=pq_it1 & it2=pq_it2 & it3=pq_it3
;
;  first plot
;
q=[0,.67]
if yr_pq_loc[1] le 2. then q2=[0,.58] else q2=q
default,yr_pq_loc,[0,2]
default,pq_label,'!6'
!x.title='!8q!6(!8t!6)'
!y.title='!8p!6(!8t!6)  and  !8p!6!dLS!n(!8t!6)'
plot,q,2*(1-2*q),yr=yr_pq_loc
bet=0. & oplot,q,q*(1.+bet),li=1 ;(beta=0 line)
bet=1. & oplot,q,q*(1.+bet),li=1 ;(beta=1 line)
bet=2. & oplot,q2,q2*(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,10*q-2.,li=2
default,xoutLS,.38
if yr_pq_loc[1] gt 2. then xyouts,xoutLS,2.4,'!8p!6!dLS!n'
;oplot,q,q*0+.57,li=2 ;(horizontal p=1/2 line)
xyouts,.22,1.45,'!7b!6=4',siz=siz
xyouts,.49,1.85,'!7b!6=3',siz=siz*.8
xyouts,.55,1.52,'!7b!6=2',siz=siz*.8
xyouts,.53,1.22,'!7b!6=1',siz=siz
xyouts,.53,0.62,'!7b!6=0',siz=siz
;xyouts,.88,1.80,pq_label,siz=siz
default,xpq_label,0.91
default,ypq_label,2.05
xyouts,xpq_label,ypq_label,pq_label,siz=siz
!x.title='!6'
;
print,'NOTE: pq.sav must have been written by pkft.pro'
restore,'pq.sav'
;
default,w,0.
nt=n_elements(qM)
it2=nt-1
for it=it1,it2,it3 do begin
  psiz=.4+1.2*float(it-it1)/float(it2-it1)
  circ_sym,psiz,0
  circ_sym,psiz,1
  oplot,[1,1]*qM[it],[1,1]*pM[it],col=122,ps=8
  oplot,[1,1]*qM[it],[1,1]*pL[it],col= 55,ps=8
  print,it,exp(lnt(it))*urms0*k0,exp(lnt(it))*vA0*k0,qM[it],pM[it],pL[it]
  ;print,it,Ret[it],ttau[it],qM[it],pM[it],b1M[it]
  wait,w
endfor
stop
;
default,read_other,0
if read_other then begin
  restore,'../'+dir_other+'/pq.sav'
  nt=n_elements(qM)
  it2=nt-1
  for it=it1,it2,it3 do begin
    psiz=.4+1.2*float(it-it1)/float(it2-it1)
    circ_sym,psiz,0
    oplot,[1,1]*qM[it],[1,1]*pM[it],col=122,ps=8
    oplot,[1,1]*qM[it],[1,1]*pL[it],col= 55,ps=8
    print,it,exp(lnt(it))*urms0*k0,exp(lnt(it))*vA0*k0,qM[it],pM[it],pL[it]
    wait,w
  endfor
endif
print,'qM[it],pM[it],pL[it]'
;
;@postproc_pq_loc
;
cwd,run
print,"$mv idl.ps ~/tex/mhd/HallMHD/fig/pq_"+run+".ps"
END