;******************
;  Get the cband radar
;******************
do_cband_distance='no'
if do_cband_distance eq 'yes' then begin
  print,'reading cband'
  dbz_grid_cband=!null
  read_cband,julian_day_ship[result_ship[0]],dbz_grid_cband,dbz_mask_cband,x_out,y_out,$
    nbins_cband,rstart,rscale,julian_day_cband,elangle_cband
  if dbz_grid_cband ne !null then begin
    ;  Plot reflectivity
    pnum=1
    data_var=dbz_grid_cband
    ;result=where(data_var ne -9999,count)
    dmin=-20.0;min(data_var[result])
    dmax=50.0;max(data_var[result])
    data_image=bytscl(data_var,top=top_color,min=dmin,max=dmax)
    ;result=where(data_var eq -999.,count)
    result=where(data_var eq -9999,count)
    if count gt 0 then data_image[result]=253
    result=where(dbz_mask_cband eq -8888,count)
    if count gt 0 then data_image[result]=255
    p1=image(data_image,/current,image_dimensions=[isx,isy],position=pos[pnum,*],$
      rgb_table=mytable)
    c1=contour(data_var,x_out/1000.0,y_out/1000.0,ytitle='Y (km)',xtitle='X (km)',$
      xstyle=1,ystyle=1,font_size=14,/nodata,/current,position=pos[pnum,*],$
      irregular=0,xtickdir=1,ytickdir=1,font_name='Helvetica',axis_style=2)
    ;  plot range rings
    j=0
    for r=0,nbins_cband[j]-1,150 do begin
      print,'range ring',r
      ;  radius in km
      range_ring=(rstart[j]+float(r)*rscale[j])/1000.0
      result=circle(0,0,range_ring)
      p2=plot(result,/current,/overplot)
    endfor
    c2=colorbar(orientation=0,rgb_table=mycbtable,range=[dmin,dmax],title='Z (dBZ)',$
      position=reform(cbpos[pnum,*]),font_size=14,font_name='Helvetica',tickdir=1,border_on=1)
    caldat,julian_day_cband[0],mm,dd,yy,hh,mi,ss
    date_time=string(yy,format='(I04)')+string(mm,format='(I02)')+string(dd,format='(I02)')+$
      ' '+string(hh,format='(I02)')+':'+string(mi,format='(I02)')+':'+string(ss,format='(I02)')
    scan_str=string(elangle_cband[i],format='(F05.2)')
    t1=text(pos[pnum,0]+dx,pos[pnum,3]+2*dy,'el='+scan_str,font_size=fs1)
    t1=text(pos[pnum,0]+14*dx,pos[pnum,3]+2*dy,date_time,font_size=fs1)
    stop
  endif else begin
    pnum=1
    t1=text(pos[pnum,0]+dx,pos[pnum,3]-2*dy,'no cband file',font_size=fs1)
  endelse
  ;p0.save,'temp_map.png',height=pydim
endif  ;end of do cband

do_cband_latlon='no'
if do_cband_latlon eq 'yes' and dbz_grid_cband ne !NULL then begin

  ;  Done earlier for range rings
  ;  read_cband_latlon,julian_day_ship[result_ship[0]],dbz_grid_cband,colon_cband,lat_cband,$
  ;    nbins_cband,rstart,rscale,julian_day_cband,elangle_cband,ulat,llat,llon,rlon,clon,clat

  ;  Plot reflectivity
  pnum=1
  data_var=dbz_grid_cband
  dmin=dmin_radar;-20.0;min(data_var[result])
  dmax=dmax_radar;50.0;max(data_var[result])
  data_image=bytscl(data_var,top=top_color,min=dmin,max=dmax)
  ;result=where(dbz_grid_cband gt -8000 and dbz_grid_cband le -7000,count)
  result=where(dbz_grid_cband gt -8000 and dbz_grid_cband lt dmin,count)
  if count gt 0 then data_image[result]=255  ;gray
  result=where(dbz_grid_cband gt -10000 and dbz_grid_cband le -9000,count)
  if count gt 0 then data_image[result]=253  ;white

  p0=image(data_image,colon_cband,lat_cband,limit=[llat,lcolon,ulat,rcolon],$
    /current,dimensions=[pxdim,pydim],grid_units='degrees',$
    center_longitude=(rcolon-lcolon)/2.0,$
    rgb_table=mytable,$;,min_value=dmin,max_value=dmax,$
    map_projection='Mercator',$
    position=pos[pnum,*],font_size=fs1)

  p1=mapcontinents(color='black')
  p0.mapgrid.box_axes=1
  p0.mapgrid.box_color='black'
  p0.mapgrid.color='black'
  p0.mapgrid.label_color='black'
  p0.mapgrid.linestyle='dot'
  p0.mapgrid.label_position=0

  ;  plot range rings
  lon_cir=fltarr(360)
  lat_cir=fltarr(360)
  for r=1,4 do begin            ;radius of earth in meters
    arc_dist=(float(r)*25000.0)/(6371D*1000D)  ;m
    for s=0,359 do begin
      result=ll_arc_distance([clon,clat],arc_dist,s,/degrees)
      lon_cir[s]=result[0]
      lat_cir[s]=result[1]
    endfor
    p7=plot(lon_cir,lat_cir,/current,/overplot)
  endfor

  c2=colorbar(orientation=1,rgb_table=mycbtable,range=[dmin,dmax],title='Z (dBZ)',$
    position=reform(cbpos[pnum,*]),font_size=14,font_name='Helvetica',tickdir=1,border_on=1)
  caldat,julian_day_cband[0],mm,dd,yy,hh,mi,ss
  date_time=string(yy,format='(I04)')+string(mm,format='(I02)')+string(dd,format='(I02)')+$
    ' '+string(hh,format='(I02)')+':'+string(mi,format='(I02)')+':'+string(ss,format='(I02)')
  scan_str=string(elangle_cband[0],format='(F05.2)')
  t1=text(pos[pnum,0]+dx,pos[pnum,1]-2*dy,'el='+scan_str,font_size=fs1)
  t1=text(pos[pnum,0]+14*dx,pos[pnum,1]-2*dy,date_time,font_size=fs1)

  ;p0.save,'temp_map.png',height=pydim

endif