;
function mkclevels,fmin,fmax,cmin,cmax,cint
;
; Set up contour min,max,interval. Args fmin,fmax are input,
; args cmin,cmax,cint are intent inout.
;
;print,'Enter mkclevels: fmin,max=',fmin,fmax,' cmin,max=',cmin,cmax,$
;  ' cint=',cint
;
; Check for valid fmin,fmax (e.g., all data missing)
;
if fmin ge fmax then begin
  print,'>>> WARNING mkclevels: fmin ge fmax: fmin=',fmin,' fmax=',fmax
  return,[0.]
endif
;
; Default number of levels:
;
;nlevdef = 12
nlevdef = 10
;
  nlevels = nlevdef
  levels = fltarr(nlevels)
  conint = (fmax-fmin) / (nlevels-1)
  for i=0,nlevels-1 do levels(i) = fmin+i*conint
  if cmin eq cmax then begin
    cmin=fmin & cmax=fmax
  endif
  cint = conint
; print,'mkclevels: cmin,max,int=',cmin,cmax,cint,' levels=' & print,levels
  return,levels
end