pro isotropized_angles,flux,sizes,err
; Plots q_precip vs loss-cone angles for the case where the pitch-angles are
;  isotropized on time scales less than the bounce time.

n=100
alp1=(findgen(n)/n)*180.*!dtor+1.d-5

alpha1=alp1#replicate(1,n)
alpha2=transpose(alpha1)

q_prec = (1+cos(alpha1))*(3-cos(alpha2))/(4.-(1-cos(alpha1))*(1-cos(alpha2)))
xra=[0,180]
lev=findgen(20)/10.
window,0

xt='ALPHA1' & YT='ALPHA2'
tit='ISOTROPIZED PRECIPITATION FRACTION AS FCN OF LOSS-CONE ANGLES'
contour,q_prec,alp1/!dtor,alp1/!dtor,xra=xra,yra=xra,/xst,/yst,lev=lev,xtit=xt,ytit=yt,C_annotation=string(lev,format='(20f4.2)'),TIT=TIT,/iso


window,1
tit='FLUX RATIO VS LOSS-CONE ANGLES'
flux_ratio=(1+cos(alpha1))*(3-cos(alpha2))/(1+cos(alpha2))*(3-cos(alpha1))
lev=[.1,0.2,0.3,0.5,1,2.0,3.0,5.0,10,100.,1000.]
contour,(flux_ratio),alp1/!dtor,alp1/!dtor,lev=lev,xra=xra,yra=xra,/xst,/yst,/iso,xtit=xt,ytit=yt,C_annotation=string(lev,format='(20f8.3)'),TIT=TIT

!p.charsize=1
yt='WIDTH RATIO' & xT='FLUX RATIO'
ti='FOOTPOINT ASYMMETRIES FOR ISOTROPIZED LOSS-CONE ANGLES'
loadct,2

if !d.name eq 'X' then color=[254,25,50,75,100,125,150,175,200,225] $
else color=[1,25,50,75,100,125,150,175,200,225]
 n=45
for j=1,9 do begin  ; constant-alpha2 curves
  alpha2=j*10*!dtor + 1.D-5
  width_ratio=sin(alpha1)/sin(alpha2)
  flux_ratio=(1+cos(alpha1))*(3-cos(alpha2))/(1+cos(alpha2))*(3-cos(alpha1))
  w=where(flux_ratio ge 0.0001)
  if j eq 1 then  plot,flux_ratio(w),width_ratio(w),xtit=xt,ytit=yt,YRA=[0.1,1],XRA=[0.1,1],color=color(0),tit=ti,PSYM=-6,symsize=1,/xst,/yst;,/xlog,/ylog
  oplot,flux_ratio(w),width_ratio(w),color=color(j),psym=-6,symsize=1
	oploterr,flux,sizes,err,psym=5,type=3

  if j eq 1 then xyouts,0.31,0.15,'!7a!3!D2!N=',/norm,color=color(0),charsize=0.9
  xyouts,0.30+0.05*j,0.15,strcompress(10*j),/norm,color=color(j),charsize=0.9

endfor

end