PRO findazza,i,c0,r0,az,za,out,debug=debug

az=0 & za=0
s=SIZE(i)
meancent,i,c,r,npix
IF ((2500 LT npix) AND (npix LT 10000)) THEN BEGIN
    dc=c0-c
    dr=r0-r
    d=SQRT(dc^2+dr^2)
    za=(BYTE(d+20)/40)*20
    IF za NE 0 THEN az=FIX(ROUND(!RADEG*ATAN(dc,dr)))
    IF KEYWORD_SET(debug) THEN BEGIN
	PRINT,c,FORM='("c= ",F8.2)'
	PRINT,r,FORM='("r= ",F8.2)'
	PRINT,dc,FORM='("dc=",F8.2)'
	PRINT,dr,FORM='("dr=",F8.2)'
	PRINT,d,FORM='("d= ",F8.2)'
	PRINT,az,FORM='("az=",I4)'
	PRINT,za,FORM='("za=",I4)'
	ENDIF
    out=STRING(za,az,FORMAT='(I3," za, ",I4," az")')
    ENDIF ELSE out="  no LBS image!"
IF KEYWORD_SET(debug) THEN PRINT,out
RETURN
END