;PV-Wave routine to produce a polar dial.
pro polar_dial_n,icglat,ihemi
;icglat = 0 ;0=ilat, 1= cglat, 2=glat labeling
!p.font=-5
pi = 3.141592654
csize=1.00
angle  = fltarr(721)
angle = 2.*pi*indgen(721)/720.
radius = replicate(10.,721) 
plot,/polar,radius,angle,xrange = [-47.,47.],xstyle=5,$
 yrange=[-47.,47.],ystyle=5,position = [0.2,0.1,0.8,0.9];,/ISOTROPIC
for i=20.,40.,10. do begin
  radius = replicate(i,721)
  oplot,/polar,radius,angle
endfor
oplot,[0.,0.],[-45.,45.]
oplot,[-45.,45.],[0.,0.]
for i=10,40,10 do begin
 tang = -45.*pi/180.
 if (ihemi eq 1) then begin
   txt = ' '+strtrim(string(90-i),2) + '!9%!5'
 endif else begin
   txt = ' '+strtrim(string(i-90),2) + '!9%!5'
 endelse
 if (i eq 40 and icglat eq 0) then txt = txt + ' ILAT'
 if (i eq 40 and icglat eq 1) then txt = txt + ' MLAT'
 if (i eq 40 and icglat eq 2) then txt = txt + ' GLAT'
 xyouts,float(i)*cos(tang),float(i)*sin(tang),txt,charsize=csize
endfor
xyouts,0.,-47.5,'0!CMLT',charsize=csize,alignment=0.5
xyouts,46.,0.,'6',charsize=csize,alignment=0.
xyouts,0.,46.,'12',charsize=csize,alignment=0.5
xyouts,-46.,0.,'18',charsize=csize,alignment=1.
return
end