;
;   Simple IDL procedure to plot 1/2 degree netCDF UWM/COADS 
;   surface marine fluxes.
;
;   Usage:  uwmplot5,time,file_name
;
;
;   Input:
;
;     time       for the monthly climatology, enter a number
;                  from 1 to 12 (jan=1, feb=2,...,dec=12)
;                for the weekly climatology, enter a number
;                  from 1 to 52 for each week of the year.
;   file_name    netCDF file name (e.g., 'sst_clm.nc')
;
;   Example:
;
;   IDL> uwmplot5,a,12,'sst_clm.nc'
;
; See also: uwm_plot.pro to plot the 1x1 degree UWM/COADS files.
;
;
;   Example:
;
;   IDL> uwm_plot,12,'sst_clm.nc'
;
;......................................................................

PRO uwmplot5,time,filen

; Read the data
; -------------
  uwmread5,a,time,filen


; Create longitude/latitude
; -------------------------
;  nx = 720
;  ny = 360
;  lat = FINDGEN(ny)*0.5 - 89.75
;  lon = FINDGEN(nx)*0.5 + 0.25


; Define contour levels
; ---------------------
  aland = 1e10
  a(where(a eq -1e10)) = aland
  nct = 30
  cmin = min(a(where(a ne aland)))
  cmax = max(a(where(a ne aland)))
  help,cmin,cmax
  ci = ( cmax - cmin ) / (nct - 1)
  lev = FINDGEN(nct) * ci + cmin  


; Plot contours
; -------------
  LOADCT,3
  ni = 800
  nj = 400
  WINDOW,/free,xsize=ni,ysize=nj

  CONTOUR,a,/cell_fill,levels=lev,$
          max_value=aland,$ 
          ystyle=4,xstyle=4,$
          xmargin=7,ymargin=2

; All done
; --------
  end