;  routine to display loop of averaged sobel-filtered images

;  coordinates for active part of image
avimg = 1.

restore,'/swing/d/lites/eclipse98/peppi/sh_genseq10.save'
nxx = sizeof(avimg,1) & nyy = sizeof(avimg,2) 
numimgs = sizeof(avimg,3)


;;; Dimensions of the final "movv" images are 300 by 228
;;; Min is 0, and max is ~263 -- movv is a float array.

mn= 0. & mx=263.  & xsz=300 & ysz=228  & snn=' '
mn= 0. & mx=225.  & xsz=300 & ysz=228  & snn=' '

mxC=180
set_plot,'z'
device,set_resolution=[xsz,ysz+60],set_colors=mxC,z_buffering=0

    loadct,0
    tvlct,red,green,blue,/get


    ;  loop over frames
    ii=1
    for k = 0,numimgs-1 do begin ;{

	gncor = avimg(*,*,k)

	;  block out region above moon's limb to avoid confusion
	for i = 0,nxx-1 do gncor(i,moonlim(i):nyy-1) = 0

	; reduce size of image
	movv = gncor(3:nxx-4,*)




	if     ( ii lt 10  ) then snn='0'+string(ii) $
	else                      snn='' +string(ii)
	snn=strcompress(snn,/remove_all)



	erase
	tv,bytscl(movv,min=mn,max=mx,top=mxC)

	xyouts, 5 , 14*3+ysz , 'February 26, 1998 eclipse, Curacao.',$
            /device,charsize=1.0,color=mxC
	xyouts, 5 , 14*2+ysz , 'Peppi image '+snn+', averaged, enhanced.',$
	    /device,charsize=1.0,color=mxC
	xyouts, 5 ,    5+ysz , 'High Altitude Observatory, NCAR/UCAR.',$
	    /device,charsize=1.1,color=mxC
 
	zzz=tvrd()

	
	write_gif,'Images/peppi98_'+snn+'.gif',zzz,red,green,blue

	ii=ii+1

     
    endfor ;}
     

device,/close
set_plot,'x'

tv,zzz

end