
PRO IC,a,ic,text,clevel,clabel,csz,min_img,max_img,noerase=noerase

 sza = size(a)

 xb = sza(1)/csz
 yb = sza(2)/csz
 
 IF (sza(1) MOD csz EQ 0) AND (sza(2) MOD csz EQ 0) THEN BEGIN 
   b = REBIN(a,xb,yb)
 ENDIF ELSE BEGIN
   b = CONGRID(a,FIX(xb),FIX(yb),/INTERP)				
 ENDELSE
 
 xb = csz*FINDGEN(FIX(xb))
 yb = csz*FINDGEN(FIX(yb))

 ;if max(b) lt 1 then clevel = clevel * 10.0^(-9.0)

 CASE ic OF
  0: BEGIN	; Image & contours


	;plot,[0,0],[0,1] ;& erase
	px = !x.window * !d.x_vsize
	py = !y.window * !d.y_vsize                 
	sx = px(1) - px(0)+1
	sy = py(1) - py(0)+1
		   

	TV,POLY_2D(BYTSCL(a>min_img<max_img),[[0,0],[sza(1)/sx,0]],$
	   [[0,sza(2)/sy],[0,0]],0,sx,sy),px(0),py(0)
        print,!p.multi(0)
        !P.MULTI(0) = !p.multi(0)-1
print,!p.multi(0)
	contour,b,xb,yb, $
	        xstyle = 1,ystyle = 1, $
	        levels = clevel, c_labels = clabel, $
		title = text, c_color = [100,200], $
		/spline,/noerase;,POSITION=[0,0,1,1]
        !P.MULTI(0) = !p.multi(0)+2
print,!p.multi(0)
     END                 
  1: BEGIN	; Image	
  	       
;	plot,[0,0],[0,0] & erase
	px = !x.window * !d.x_vsize
	py = !y.window * !d.y_vsize
	sx = px(1) - px(0)+1
	sy = py(1) - py(0)+1
	
	TV,POLY_2D(BYTSCL(a>min_img<max_img),[[0,0],[sza(1)/sx,0]],$
	   [[0,sza(2)/sy],[0,0]],0,sx,sy),px(0),py(0)
	   
	contour,b,xb,yb, $
	        xstyle=1,ystyle=1, $
	        levels=[-100000], $
	        title=text, $
		/noerase,/spline

     END
  2: BEGIN
;	erase
;	contour,b,xb,yb, $
;	        xstyle = 1,ystyle = 1, $
;	        levels = clevel, c_labels = clabel, $
;		title = text, c_color = [100,200], $
;		/noerase,/spline
        IF NOT KEYWORD_SET(noerase) THEN erase
        contour,b,xb,yb, $
	        xstyle = 1,ystyle = 1, $
	        levels = clevel, c_labels = clabel, $
		title = text, $
		/spline

     END     
     ELSE: donothing=0
 ENDCASE
		 
RETURN
END
