firsttime=0
if not firsttime then goto,step2
current_zoom=.5
; for title
font_size=24
font_style=1 ; bold
insert_colorbar=1 ; if [2] then bottom left corner  of color bar (normalize coord).
rgb_table=33; blue-red .. could be a 3,256 vector
dlim=500
xrange=[-dlim,dlim]
yrange=[-dlim,dlim]
xtickunits=''
ytickunits=''
xtickvalues=0
ytickvalues=0
xticklayout=0
yticklayout=0
xtitle='South [feet]'
ytitle='North [feet]'
el='West'
eR='East'
zoom_on_resize=1
view_zoom=.25 ; start low so we can see
;
;title='11mar20: radial errors (measured - fit) 3 param fit. x,y center offset set t0 0.'
lmaxErr='5cm'
title='11mar20: radial errors (measured - fit) 3 param fit. Keep points < 5 cm error.'
if 0 then begin
 	lmaxErr='1.2cm'
	title='11mar20: radial errors (measured - fit) 3 param fit. maxRadialErrs<1.2cm'
endif
; notes.. axis is not showing up..
;         in the itools visualization properties, you need to 
;         1. click on axex
;         2. in table click style and then box show
;         
;iimage,imgClip,xout,yout,$
print,"Call iimage"
iimage,img,xout,yout,$
;   current_zoom=current_zoom,$
    /zoom_on_resize,$
    font_size=font_size,$
    font_style=font_style,$
    insert_colorbar=insert_colorbar,$
    rgb_table=rgb_table,$
    identifier=img4id,$
    title=title,$
	view_zoom=view_zoom,$
    xrange=xrange,yrange=yrange,$
    xtickvalues=xtickvalues,ytickvalues=ytickvalues,$
    xtickunits=xtickunits,ytickunits=ytickunits,$
    xticklayout=xticklayout,yticklayout=yticklayout


	imgid=igetid('image')
;igetproperty,imgid,_registered=imgpropall
;print,"Reposition the window size to what you want the enter to continue"
;key=checkkey(/wait)
;
; axis
;
print,"start axis setup"
axisid=igetid("axes")
isetproperty,axisid,style=2
print,"center image in window then restart with firsttime=0"
stop
step2:
radius=500.
transparency=60
ls=0 ; solid line
; bombs using imgid
oplotid=1
for ix=-475,475,25 do begin &$
    thT=acos(ix/radius) &$
    ymax=abs(sin(thT)*radius) &$
    xx=[ix,ix] &$
    yy=[-ymax,ymax] &$
    iplot,xx,yy,overplot=oplotid,xmajor=0,ymajor=0,xgrid=6,ygrid=6,transparency=transparency,linestyle=ls &$
endfor
;
; note after  drawing the vertical lines axis0 is down ,axis 1 is up, axis2 is left, axis 3 is right
;
axIdAr=strarr(4)
for i=0,3 do begin &$
	nm=string(format='("Axis ",i1)',i) &$
	axIdAr[i]=igetid(nm) &$
endfor
; so we switch index 1,2 
tmp=axIdar[1]
axidar[1]=axidar[2]
axidAr[2]=tmp
igetproperty,axidar[0],_registered=axis1propall
font_size=24
tickdir=1  ; out  ; for x,y
ticklen=.03
majorticks=6
for i=0,3 do begin &$
	if i lt 2 then begin &$
		print,i," font_size" &$
		isetproperty,axIdAr[i],font_size=font_size &$
		print,i," font_style" &$
    	isetproperty,axIdAr[i],font_style=font_style &$
		print,i," tickdir" &$
		isetproperty,axIdAr[i],tickdir=tickdir &$
	endif else begin &$
		print,i," tickdir" &$
		isetproperty,axIdAr[i],tickdir=0 &$
	endelse &$
	print,i," ticklen" &$
    isetproperty,axidar[i],ticklen=ticklen &$
	print,i," major" &$
    isetproperty,axidar[i],major=majorticks &$
endfor
;
; generate grid of main cables every 25 feet
;
; colar bar
print,"start colorbar"
colBarId=igetid('colorbar')
igetproperty,colbarid,_registered=cbpropall
font_size=24
border=1
tickdir=1
ticklen=.2
subticklen=.4
orientation=1  ; verticaal
textpos=1  ; right above
cbtitle='radial errors [cm]'
isetproperty,colbarid,orientation=orientation
isetproperty,colbarid,textpos=textpos
isetproperty,colbarid,font_size=font_size
isetproperty,colbarid,font_style=font_style
isetproperty,colbarid,border_on=1
isetproperty,colbarid,tickdir=tickdir
isetproperty,colbarid,ticklen=ticklen
isetproperty,colbarid,subticklen=subticklen
isetproperty,colbarid,axis_title=cbtitle
;-----------------------------------------------------------
;    label the edges still had to move them around a bit
print,"label the edges"
frac=20
text='East'
x=radius+frac
y=0
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
igetproperty,txtid,_registered=txtpropall
text='Y [feet]'
x=-radius -  6*frac
y=0
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
;
x=0 - 3*frac
y=-radius   - 4*frac
text='(West) X [feet]  (East)'
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
x=-1*frac 
y=radius + 1*frac
text='North'
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
;-----------------------------------------------------------
; text on right side
print,"start text on right side"
yeps=30
x=radius + frac
y0=470
text='Radius Used: 265.176 [m]'
itext,text,x,y0,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
;
y=y0-yeps
; cent set in doit center used
text=string(format='("Center:[",f6.4,",",f6.4,",",f8.4,"]")',cent)
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
;
y=y0-2*yeps
text='Keep points with error < ' + lmaxerr
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
;
y=y0-3*yeps
text='Dark blue is no data'
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
;
y=y0-4*yeps
text='grid 2000x2000'
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
;
y=y0-5*yeps
text='inverseDistance with 40pnts'
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
;
; label the 100 foot main cables
;
nmain=dishtarmain(cblI,/tdonly)
ok=intarr(nmain)
for i=0,nmain-1 do ok[i]=strmid(cblI[i].nm,1,1) eq '0'
ii=where( ok eq 1,cnt)
;
;now add the labels
;
print,"start labeling the main cables"
ctmp=cblI[ii]
jj=sort(ctmp.xth)
ctmp=ctmp[jj]
print,ctmp.xth
y0=-240
yeps=15
fs=16
for i=0,cnt-1 do begin &$
	x=ctmp[i].xth &$
	y=y0 - (i mod 2)*yeps &$
	print,x,y0," ",ctmp[i].nm  &$
	itext,ctmp[i].nm,x,y,target_identifier=imgid,/data,/visualization,font_style=font_style,font_size=fs &$
endfor
;
; label upper lower edge of color bar
;
x= radius + 3.5*frac
y=263
text='below dish'
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size
y=-500
text='above dish'
itext,text,x,y,target_identifier=imgid,/data,/visualization,identifier=txtid,font_style=font_style,font_size=font_size

print,"you need to manually position the color bar and above,below labels"
;
end