<div dir="ltr">Hello,<div><br></div><div>I am plotting a Lambert Conformal map using res@mpLimitMode set to "Corners". The map limits are configured to zoom in on a subset of the full domain of the gridded dataset in the x-y dimensions (lat-lon). I am then calling gsn_csm_contour_map to plot the data.</div><div><br></div><div>What I would like to do, is use the "Corners" resources, to create a polygon so that I can perform statistics on the data only within the view shown in the map that is drawn by gsn_csm_contour_map. The problem is, I have not figured out how gsn_csm_contour_map determines what the bounds are for the domain drawn when using the "Corners" resources.</div><div><br></div><div>My first attempt was to create a polygon from the "Corners" resources (res@mpLeftCornerLonF, res@mpLeftCornerLatF, res@mpRightCornerLonF, and res@mpRightCornerLatF) and then use gc_inout to identify which points from the gridded dataset are within the map limits. When I plot the resultant boolean mask however, the points "inside" the polygon do not directly match the view shown in the bounding box from gsn_csm_contour_map (see attached figure, Mask1.png).</div><div><br></div><div>The end goal is to make sure the statistics being computed are created from exactly those grid cells shown within the bounding region created by gsn_csm_contour_map.</div><div><br></div><div>Any suggestions or ideas are appreciated. Thanks!</div><div><br></div><div>-Dan</div><div><br></div><div><br></div></div>