pro imagecon

file = 'lonslice.dat'
openr,lu,file,/get_lun, /swap_if_little_endian
 
nlat = 0l
readu,lu,nlat
lats = fltarr(nlat)
readu,lu,lats
 
nlev = 0l
readu,lu,nlev
levs = fltarr(nlev)
readu,lu,levs
 
data = fltarr(nlat,nlev)
readu,lu,data
free_lun,lu

;print,'nlev=',nlev,' levs = ' & print,levs
;print,'nlat=',nlat,' lats = ' & print,lats
;print,'data min,max=',min(data),max(data)

xur = .90 & yur = .90
xll = .15 & yll = .31
pos = [xll,yll,xur,yur]

xtitle = 'LATITUDE'
ytitle = 'ZP (LN(P0/P))'

window,xsize=500,ysize=500
loadct,5

image = bytscl(data,top=!d.table_size)

xyll = convert_coord([xll],[yll],/norm,/to_device)
xll_dev = xyll(0,0) & yll_dev = xyll(1,0)
xyur = convert_coord([xur],[yur],/norm,/to_device)
xur_dev = xyur(0,0) & yur_dev = xyur(1,0)
image=congrid(image,xur_dev-xll_dev,yur_dev-yll_dev,/interp,/minus_one)
tv,image,xll,yll,/norm,xsize=xur-xll,ysize=yur-yll
clevels=[160.,170.,180.,190.,200.,210.,220.,230.,240.,250.,260.,270.,280.]
contour,data,lats,levs,/noerase,/nodata,position=pos,$
  xtitle=xtitle,ytitle=ytitle,xstyle=1,ystyle=1,ticklen=-.015
contour,data,lats,levs,/noerase,position=pos,/follow,xstyle=5,ystyle=5,$
  levels=clevels

;px = !x.window * !d.x_vsize
;py = !y.window * !d.y_vsize
;sx = px[1] - px[0]+1
;sy = py[1] - py[0]+1
;image = congrid(image,sx,sy)
;tv,image,px[0],py[0]
;clevels=[160.,170.,180.,190.,200.,210.,220.,230.,240.,250.,260.,270.,280.]
;contour,data,lats,levs,/noerase,position=pos,/follow,xstyle=1,ystyle=1,$
;  levels=clevels

;px = !x.window * !d.x_vsize
;py = !y.window * !d.y_vsize
;image = rebin(image(*,*),nlat*10,nlev*15)
;sz = size(image)
;tv,image,px[0],py[0]
;clevels=[160.,170.,180.,190.,200.,210.,220.,230.,240.,250.,260.,270.,280.]
;contour,data,lats,levs,/noerase,/follow,xstyle=1,ystyle=1,levels=clevels,$
;  position=[px[0],py[0],px[0]+sz[1]-1,py[0]+sz[2]-1],/device

end