pro show_pps_iss,rr

if(n_params le 0) then begin
    print,'show_pps_iss,rr'
    print,'display Voyager tau, I_lit, I_unlit<
    print,'for the range rr=[rmin,rmax] (in 1000 km)'
return
endif

nypix=100

;PPS optical depth
restore,'pps_tau_vs_r.sav'
r=r_pps/1000.
ind=where(r ge rr(0) and r le rr(1),nxpix)
ima=fltarr(nxpix,nypix)
for i=0,nypix-1 do begin
    ima(*,i)=tau_pps(ind)
endfor



;imaging
restore,'vgr_iss.save'
r_lit=r_lit/1000.
ind_lit=where(r_lit ge rr(0) and r_lit le rr(1),nxpix_lit)
ima_lit=fltarr(nxpix_lit,nypix)
for i=0,nypix-1 do begin
    ima_lit(*,i)=if_lit(ind_lit)
endfor

r_unlit=r_unlit/1000.
ind_unlit=where(r_unlit ge rr(0) and r_unlit le rr(1),nxpix_unlit)
ima_unlit=fltarr(nxpix_unlit,nypix)
for i=0,nypix-1 do begin
    ima_unlit(*,i)=if_unlit(ind_unlit)
endfor

ima_lit=congrid(ima_lit,nxpix,nypix)
ima_unlit=congrid(ima_unlit,nxpix,nypix)

nwin
tile_to_pos,1,4,page=[.05,.05,.95,.95],pos,fsize=0.8
!p.position=pos(0,*)
plot,r(ind),tau_pps(ind),ytitle='tau',yr=[0,2],ys=1
!p.position=pos(1,*)
tvplot2,ima,xx0=r(ind),/noerase,/nonew,ytitle='tau',ycharsize=0.01
!p.position=pos(2,*)
tvplot2,ima_lit,xx0=r(ind),/noerase,/nonew,ytitle='I/F lit side',ycharsize=0.01
!p.position=pos(3,*)
tvplot2,ima_unlit,xx0=r(ind),/noerase,/nonew,ytitle='I/F unlit side',xtitle='radius (1000 km)',ycharsize=0.01
!p.multi=0
!p.position=0.
end