pro pltce
  file = 'for_rsi.dat' 
  openr,lu,file,/get_lun, /swap_if_little_endian
;
  nlon = 0l
  readu,lu,nlon
  lons = fltarr(nlon)
  readu,lu,lons
;
  nlat = 0l
  readu,lu,nlat
  lats = fltarr(nlat)
  readu,lu,lats
;
  data = fltarr(nlon,nlat)
  readu,lu,data
  free_lun,lu
;
; print,'lons=' & print,lons
; print,'lats=' & print,lats
; print,' '
; print,'data min,max=',min(data),max(data)
;
  window,xsize=500,ysize=500
  loadct,5
; loadct,6
;
; Set map position (normalized coords):
;
  xll = .15 & yll = .48
  xur = .95 & yur = yll+.5*(xur-xll)
  mappos = [xll,yll,xur,yur]
  dlat = lats[1]-lats[0] ; assume constant delta lat
  poledeg = (dlat/2.)/180.*(yur-yll) ; distance between top lat and pole
  conpos = [xll,yll+poledeg,xur,yur-poledeg]
;
  xtitle = 'LONGITUDE (DEGREES GEOGRAPHIC)'
  ytitle = 'LATITUDE (DEGREES GEOGRAPHIC)'
  xticlabs = ['-180','-120','-60','0','60','120','180']
  yticlabs = ['-90','-60','-30','0','30','60','90']

  map_set,/cylindrical,position=mappos,/noborder
  image = bytscl(data,top=!d.table_size-3,min=fmin,max=fmax)+1

; imscale = 8
; image = congrid(image,nlon*imscale,nlat*imscale,/interp)
; image = map_image(image,startx,starty,xsize,ysize,$
;   latmin=lats[0],latmax=lats[nlat-1],$
;   lonmin=lons[0],lonmax=lons[nlon-1],compress=1)

; lons_patch = fltarr(nlon-1)
; lons_patch[*] = lons[0:nlon-2]
; lons_patch = findgen(nlon-1)/(nlon-2)*360-180
  
  image_patch = congrid(image,nlon-1,nlat)
  image = map_patch(image_patch,dummy,lats,$
    xstart=startx,ystart=starty,xsize=xsize,ysize=ysize)

  tv,image,startx,starty,xsize=xsize,ysize=ysize
  contour,data,lons,lats,/noerase,position=mappos,/nodata,$
    xtitle=xtitle,ytitle=ytitle,xstyle=1,ystyle=1,ticklen=-.015,$
    xticks=n_elements(xticlabs)-1,xminor=3,xtickname=xticlabs,$
    yticks=n_elements(yticlabs)-1,yminor=3,ytickname=yticlabs
  contour,data,lons,lats,/noerase,position=conpos,xstyle=5,ystyle=5,$
    /follow,nlevels=15
end