plotting subset of a grid

Hi all,
I try to extract data over senegal (country) but I cannot plot the data,
any help ?
here is my script and the content of the variable rr
==

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
in = addfile("EnsECMWFoutput2011_.grib2","r");
rr=in->tp_P11_L1_GLL0_acc(1,92,1,{12.5:17},{342:349})
wks = gsn_open_wks("ps","test") ;
gsn_define_colormap(wks,"rainbow")
res = True
res@gsnAddCyclic=False
plot = gsn_csm_contour_map(wks,rr,res)

==
ncl 11> print(rr)
Variable: rr
Type: float
Total Size: 600 bytes
            150 values
Number of Dimensions: 2
Dimensions and sizes: [lat_0 | 10] x [lon_0 | 15]
Coordinates:
            lat_0: [12.5..17]
            lon_0: [342..349]
Number Of Attributes: 14
  forecast_time0 : 48
  initial_time0_hours : 1854684
  ensemble0 : 1
  center : European Center for Medium-Range Weather Forecasts
  production_status : TIGGE Operational products
  long_name : total precipitation
  units : kg/m^2
  _FillValue : 1e+20
  grid_type : Latitude/longitude
  parameter_discipline_and_category : Meteorological products, Moisture
  parameter_template_discipline_category_number : ( 11, 0, 1, 52 )
  level_type : Ground or water surface
  type_of_statistical_processing : Accumulation
  statistical_process_duration : initial time to forecast time
===

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Jul 1 16:34:40 2012