Re: plotting subset of a grid

Hi,

Thanks for providing your short script and the output from printVarSummary.
I don't have quite enough information, though.

What kind of plot do you get when you run this script? Do you get any errors?
If you get errors, please cut-and-paste the full error and include that in
an email to nc-talk.

If you can also include the test.ps file, that would be helpful.

One thing you can try: if you are looking at a very small area on the map,
then you may need to set some resources to zoom in on the map.

The "easiest" way to do this is with these resources:

  res@mpMinLatF = min(rr&lat_0)
  res@mpMaxLatF = max(rr&lat_0)
  res@mpMinLonF = min(rr&lon_0)
  res@mpMaxLonF = max(rr&lon_0)

--Mary

On Jul 1, 2012, at 4:34 PM, ousmane ndiaye wrote:

> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jul 2 08:47:35 2012