<div dir="ltr"><div>Ajay,</div><div><br></div><div>Thanks for sending the files. I experimented with this for a good while and learned something I was not aware of. First I'd point out that when working with irregularly spaced point data, NCL must employ a triangulation of the data, and it is not uncommon to see strange-ish linear contour lines along the outer edges of a resultant contour plot, particularly if the point data has a concave shape, which your's does.<br></div><div><br></div><div>So what I realized is that shapefile_mask_data() works in *data* space, not in *graphics* space (i.e., contour lines). All of the station points in your dataset are obviously inside the boundaries of the shapefile, so it returns a mask that is all 1's (inside). Then the contouring algorithm runs against that result, which in effect is the entire original concave dataset.  The inevitable "strange-ish" contour lines form, but they are in graphics space, and hence, they are not clipped against the shapefile.</div><div><br></div><div>So...I'm not sure what to tell you. One possibility would be to interpolate a regular grid from the station data (I plotted the individual stations and they seem to be quite dense, so that the resultant grid might make a good representation). That said, I don't know offhand what facilities are available for doing that -- perhaps that's a separate post to ncl-talk.</div><div><br></div><div>Wish I had a better answer, but I hope that helps explain what it going on.</div><div><br></div><div>Rick</div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jan 8, 2020 at 8:03 PM Ajay Bankar <<a href="mailto:ajaybankar123@gmail.com">ajaybankar123@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Rick,<div>  Â  Â  Â  Â  Â The same shapefile works fine with GPM data and wrf output data but only having issue with station data. Please find attached shapefiles. I have attached one plot in which I have used same shapefile for wrf model output and GPM data. Here plots (a,b,c) represents wrf plots and (d) represent GPM data plot and (e) represent station data plot.  </div><div><br></div><div>This is the output of ncl_filedump District_Boundary.shp</div><div><br></div><div>Variable: f<br>Type: file<br>filename:        District_Boundary<br>path:        District_Boundary.shp<br>  Â file global attributes:<br>  Â  Â  layer_name : District_Boundary<br>  Â  Â  geometry_type : polygon<br>  Â  Â  geom_segIndex : 0<br>  Â  Â  geom_numSegs : 1<br>  Â  Â  segs_xyzIndex : 0<br>  Â  Â  segs_numPnts : 1<br>  Â dimensions:<br>  Â  Â  geometry = 2<br>  Â  Â  segments = 2<br>  Â  Â  num_features = 30 Â // unlimited<br>  Â  Â  num_segments = 33<br>  Â  Â  num_points = 164066<br>  Â variables:<br>  Â  Â  integer geometry ( num_features, geometry )<br><br>  Â  Â  integer segments ( num_segments, segments )<br><br>  Â  Â  double x ( num_points )<br><br>  Â  Â  double y ( num_points )<br><br>  Â  Â  double PERIMETER ( num_features )<br><br>  Â  Â  double AREA ( num_features )<br><br>  Â  Â  int64 ID ( num_features )<br><br>  Â  Â  string DIST_NAME ( num_features )<br><br>  Â  Â  double SQ_KMS ( num_features )</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 9, 2020 at 4:50 AM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Can you either i) provide the shapefile (and I'll point out that means the .shp, .shx, .dbf files), or perhaps show us the results of <br></div><div><br></div><div>   ncl_filedump District_Boundary.shp</div><div><br></div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jan 8, 2020 at 7:24 AM Ajay Bankar via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear NCL Users,<div>  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  I'm trying to make spatial plot of rainfall from station data available in csv file. I tried different options in resources but data is not masking with shapefile. </div><div>I have attached plot and data file for reference.</div><div><br></div><div>Thanks for any help!</div><div><br></div><div>Here is my script  Â Â </div><div><br></div><div>load "/shp.ncl"  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â Â </div><div><br><div><div>begin <br> ;-- shapefile name<br> shp_filename = "/scratch/Administrative_Boundary_Headquarters/District_Boundary.shp" <br><br></div><div> infile = "station_data.csv"<br> lines = asciiread(infile, -1, "string")<br><br>rain = tofloat( str_get_field(lines(1:), 7, ","))<br>rain@lat1d Â = tofloat( str_get_field(lines(1:), 5, ","))<br>rain@lon1d Â = tofloat( str_get_field(lines(1:), 6, ","))<br>;print(rain)<br> <br>;-- define sub-region, here for Karnataka<br>  Â  Â minlat = Â 11<br>  Â  Â maxlat = Â 19<br>  Â  Â minlon = Â 73.5<br>  Â  Â maxlon = Â 79.1<br> <br>;-- open workstation<br>  Â  Â wks_type Â  Â  Â  Â  Â = "pdf"<br>  Â  Â wks Â = gsn_open_wks(wks_type, "station_plot") <br> ;-- resource settings<br>  Â  Â res Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â = Â True<br>  Â  Â res@gsnDraw Â  Â  Â  Â  Â  Â  Â  Â = Â False Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ;-- don't draw plot yet<br>  Â  Â res@gsnFrame Â  Â  Â  Â  Â  Â  Â  = Â False Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ;-- don't advance frame yet<br>  Â  Â res@gsnMaximize Â  Â  Â  Â  Â  Â = Â True Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ;-- maximize plot in frame<br> <br>  ; Â res@tfDoNDCOverlay Â  Â  Â  Â = True<br><br>  Â  cmap Â  Â  := read_colormap_file("BlAqGrYeOrReVi200")<br>  Â  cmap(0,:) = (/0,0,0,0/) Â  Â ; make first color fully transparent<br>  Â  res@cnFillOn Â  Â  Â  Â  Â  Â  Â  = Â True<br>  Â  res@cnLinesOn Â  Â  Â  Â  Â  Â  Â = Â False<br>  Â  res@cnLineLabelsOn Â  Â  Â  Â  = Â False<br>  Â  res@cnFillPalette Â  Â  Â  Â  Â = cmap Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ;"BlueYellowRed"<br>  Â  res@cnLevelSelectionMode Â  = "ExplicitLevels" Â  Â  Â  Â  Â  Â  Â  ;-- manually set contour levels<br>  Â  res@cnLevels  Â  Â  Â  = (/5,20,35,50,65,80,95,110,125,140,155/)<br>  Â Â res@gsnAddCyclic Â  Â  Â  Â  Â  = False<br>  <br>  Â  res@mpGridLineColor Â  Â  Â  Â = Â 0 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ;"grey40" <br>  Â  res@mpFillOn Â  Â  Â  Â  Â  Â  Â  = Â False<br>  Â  res@mpOutlineOn Â  Â  Â  Â  Â  Â = Â False<br>  Â  res@mpGeophysicalLineColor = "black"<br>  Â  res@mpLimitMode Â  Â  Â  Â  Â  Â = "LatLon"<br>  Â  res@mpMinLatF Â  Â  Â  Â  Â  Â  Â = minlat<br>  Â  res@mpMaxLatF Â  Â  Â  Â  Â  Â  Â = maxlat<br>  Â  res@mpMinLonF Â  Â  Â  Â  Â  Â  Â = minlon<br>  Â  res@mpMaxLonF Â  Â  Â  Â  Â  Â  Â = maxlon<br>  Â  ;res@mpAreaMaskingOn Â  Â  Â  Â = True<br><br>  Â res@lbBoxMinorExtentF Â  Â  Â = Â 0.1 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ;-- decrease height of labelbar boxes<br>  Â res@pmLabelBarOrthogonalPosF = 0.08 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â ;-- move labelbar to the left side of plot<br>  Â res@lbBoxMinorExtentF Â  Â  Â = Â 0.2 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  ;-- decrease height of labelbar boxes<br><br>  Â res@pmTickMarkDisplayMode Â = "Always"<br>  Â res@tmXTOn Â  Â  Â  Â  Â  Â  Â  Â  = Â False<br>  Â res@tiMainString Â  Â  Â  Â  Â  = "STATION DATA PLOT"<br>  <br> ;-- create plot<br>  Â  plot = gsn_csm_contour_map(wks,rain,res)<br>  Â  var_mask = shapefile_mask_data(rain,shp_filename,True)<br> <br> ;-- create contours of masked data<br>  Â  plot_mask = gsn_csm_contour_map(wks,var_mask,res)<br> <br> ;-- add plot_mask to plot<br>  Â  Â lnres Â  Â  Â  Â  Â  Â  Â  Â  Â  Â = Â True<br>  Â  Â lnres@gsLineColor Â  Â  Â  Â = "black"<br>  Â  Â lnres@gsLineThicknessF Â  = Â 1.0<br> <br>  Â  Â plot@lines = gsn_add_shapefile_polylines(wks, plot_mask, shp_filename, lnres)<br><br>  ;-- draw the plot and advance the frame<br>  Â  Â draw(plot_mask)<br>  Â  Â frame(wks)<br> end<br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div><b>Thanks & Regards,<br></b></div><b>Ajay<br></b></div><b><br></b></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div><b>Thanks & Regards,<br></b></div><b>Ajay Bankar.<br></b></div><b><br></b></div></div>
</blockquote></div>