;*************************************************
; NCL Graphics: lsm_2.ncl
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" 
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("lsmh_01153.nc","r")
;************************************************
; read in soil moisture
;************************************************
  h2o = a->H2OSOI(0,0,:,:) 
; note, in LSM files, it is safer to copy the missing value from the
; file than try to write it out, since the number of decimal points on the
; file may not be the same number NCL reads. (9.99999e29)

  h2o@_FillValue = a->TAM@missing_value
  h2o@long_name=h2o@Title                   ; no long_name on this file
  h2o!0="lat"                               ; coordinate vars are called
  h2o!1="lon"                               ; "latitude" etc. plots want "lat"
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps","lsm")                ; open a ps file

  res                      = True               ; plot mods desired
  res@mpMinLonF            = -140               ; set min lon
  res@mpMaxLonF            =  -65               ; set max lon
  res@mpMinLatF            =  10.               ; set min lat
  res@mpFillColors         = (/0,9,28,25/)      ; colors to draw (default map)
  plot = gsn_csm_contour_map_ce(wks,h2o, res)   ; create plot


  end