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
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
end
;
@../Hf0_t2em5_k180a/parameters.pro
siz=1.2
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
it1=20
it3=1
;
;  first plot
;
q=[0,.67]
yr_pq_loc=[0.,2.]
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)'
plot,q,2*(1-2*q),yr=yr_pq_loc,xr=[0.,.63]
oplot,q,2*(1-q),li=3
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=1.7& oplot,q,q*(1.+bet),li=1,col=122
bet=1.5& oplot,q,q*(1.+bet),li=1,col=55
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,10*q-2.,li=2  ;(forgot what this was)
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)
siz=1.8
si2=1.5
xyouts,.22,1.45,'!7b!6=4',siz=siz
xyouts,.49,1.85,'!7b!6=3',siz=siz;*.8
xyouts,.48,1.60,'!7b!6=2',siz=siz;*.8
xyouts,.53,1.17,'!7b!6=1',siz=siz
xyouts,.45,1.38,'!7b!6=1.7',siz=si2,col=122
xyouts,.37,0.85,'!7b!6=1.5',siz=si2,col=55
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'
dir='ADsig0_k180a'
dir='ADsig0_k180b'
dir='ADsig0_k180c'
dir='ADsig0_k180d'
restore,'../'+dir+'/pq.sav'
;
w=.1
default,w,0.
nt=n_elements(qM)
it2=nt-1
print,'1: ',nt
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]
  ;print,it,Ret[it],ttau[it],qM[it],pM[it],b1M[it]
  wait,w
endfor
;
;dir_other='Hf0_t1em5_k180a_Hosk' & read_other=1
dir_other='Hf0_t2em5_k180a_Hosk' & read_other=1
default,read_other,0
if read_other then begin
  restore,'../'+dir_other+'/pq.sav'
  nt=n_elements(qM)
  it2=nt-1
print,'2: ',nt
  for it=it1,it2,it3 do begin
    psiz=.4+.7*float(it-it1)/float(it2-it1)
    circ_sym,psiz,1
    oplot,[1,1]*qM[it],[1,1]*pM[it],col=55,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]'
;
dir_other='Hf0_t2em5_k180a_Hosk_k2' & read_other=1
default,read_other,0
if read_other then begin
  restore,'../'+dir_other+'/pq.sav'
  nt=n_elements(qM)
  it2=nt-1
print,'2: ',nt
  for it=it1,it2,it3 do begin
    psiz=.4+.7*float(it-it1)/float(it2-it1)
    circ_sym,psiz,1
    oplot,[1,1]*qM[it],[1,1]*pM[it],col=155,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]'
;
p=10./13.
q=4./13.
oplot,[1,1]*q,[0.,p],li=3
oplot,[0.,q],[1,1]*p,li=3
xyouts,.314,.08,'!8q!6=4/13'
xyouts,.012,.84,'!8p!6=10/13'
;
;@postproc_pq_loc
;
cwd,run
print,"$mv idl.ps ~/P/HallCascade/Nordita/pq_"+run+".ps"
print,"$mv idl.ps ~/tex/mhd/HallNonhel/fig/pq_comp.eps"
END