;--------------------------------------------
;ring_panorama.pro
;--------------------------------------------

;1) displays the A-ring part of the Cassini panoram image with
;a brightness scale that shows the contrast reversal of the
;Mimas 5:3 bending wave region

ima_panorama=read_tiff('PIA06175_ring_panorama.tif')
nwin,xs=1200,ys=1000

loadct,3

;tvscl,ima(2,*,*)<125>113,-4530,0    ; good for showing asymmetry
tvscl,ima_panorama(2,*,*)<175>100,-4530,50     ; good for density waves

ots,"PIA06175.tif   12.12.2004  B=-4, B'=-23.2  alpha=42  (theta_sun=-40) scale=10.5km/pixel",/arp


;add the Voyager ISS profile to the image
;lit  =lit side I/F

rmax=137.
rmin=124.
show_iss,[rmin,rmax],xs=1200,ys=50,lit=lit,/noplot
!p.position=[0.,0.03,1.,0.1]
xx=rmin+findgen(1200)/1200.*(rmax-rmin)
tvplot2,sigrange(lit),xx=xx,/nonew,/noerase
!p.position=0

;-------------------------------------------------------
;the "ring_panorama.ps" was made by using wset to the above image
;and then calling psdump: 
;(in Wellesley this did not work on color, but probably I did not
; remember all the keywords needed!)
;psdump,'ring_panorama.ps',/color,/true,/copy
;-------------------------------------------------------






;also make a floating ruler
show_iss,[124,137.],xs=1200/0.95     ; ruler - VGR ISS



;------------------------------------------------------
;this is another Cassini image
;"RINGS FROM AFAR" = 572_1170_3.tif  

;just demonstrating how different the resonances can look depending
;on geometry!

ima2=read_tiff('572_1170_3.tif')
nwin,xs=699,ys=869
tvscl,ima2<200>100,0,-200
ots," 572_1170_3.tif   1.11.2004  B=-13.5, B'=-23.5  alpha=112 !c!c  (theta_sun=-120) scale=129km/pixel" 


end