<div dir="ltr">OK, I think I understand now. That particular shapefile, Glaciar, contains only point features (see dump below) with just 5 features in it. I believe the shapefile must contain polygonal regions against which to mask.<br><br><br>ncl 0> f=addfile("Glaciares.shp","r")<br>ncl 1> print(f)<br><br>Variable: f<br>Type: file<br>filename: Glaciares<br>path: Glaciares.shp<br> file global attributes:<br> layer_name : Glaciares<br> geometry_type : point<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 = 5 // unlimited<br> num_segments = 5<br> num_points = 5<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> integer Id ( num_features )<br><br> string Nombre ( num_features )<br><br>ncl 2> print(f->Nombre)<br><br><br>Variable: Nombre (file variable)<br>Type: string<br>Total Size: 40 bytes<br> 5 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [num_features | 5]<br>Coordinates: <br>Number Of Attributes: 0<br>(0) Glaciar Cipreses<br>(1) Glaciar Palomo<br>(2) Glaciar Cortaderal<br>(3) Glaciar Universidad<br>(4) Glaciar Cotón<br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jan 23, 2018 at 2:36 PM, Vanúcia Schumacher <span dir="ltr"><<a href="mailto:vanucia-schumacher@hotmail.com" target="_blank">vanucia-schumacher@hotmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_8841829022642877653divtagdefaultwrapper" style="font-size:14pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p style="margin-top:0;margin-bottom:0">Hi <br>
</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0"><span id="m_8841829022642877653result_box" lang="en"><span>The ascii file is a mask with missing values and values in some pixels, as in the attached figure.</span><br>
<span>The problem is that when I try to create another mask using the shapefile (to set only pixels in a specific region-
<span id="m_8841829022642877653result_box" class="m_8841829022642877653short_text" lang="en"><span>smaller region of the ascii file</span></span>) it returns only missing values.</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_8841829022642877653result_box" lang="en"><span><br>
</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_8841829022642877653result_box" lang="en"><span></span></span></p>
<div><span id="m_8841829022642877653result_box" class="m_8841829022642877653short_text" lang="en"><span>See printVarSummary</span></span>:<br>
<br>
Variable: data<br>
Type: float<br>
Total Size: 80496 bytes<br>
20124 values<br>
Number of Dimensions: 2<br>
Dimensions and sizes: [lat | 156] x [lon | 129]<br>
Coordinates: <br>
lat: [-34.69..-34.51]<br>
lon: [-70.36..-70.28]<br>
Number Of Attributes: 3<br>
lat2d : <ARRAY of 20124 elements><br>
lon2d : <ARRAY of 20124 elements><br>
_FillValue : -9999<br>
(0) <span style="color:rgb(255,0,0)"> min=0 max=14.9 </span><br>
==============================<wbr>==============================<wbr>==========<br>
(0) Filename: "Glaciares.shp"<br>
(0) Geometry type: point<br>
(0) # of features: 5<br>
(0) Min/max lat: -34.69/ -34.51<br>
(0) Min/max lon: -70.36/ -70.28<br>
(0) Variable names and their types:<br>
(0) geometry : integer<br>
(0) segments : integer<br>
(0) x : double<br>
(0) y : double<br>
(0) Id : integer<br>
(0) Nombre : string<br>
(0) =============================<wbr>==============================<wbr>===========<br>
(0) <span style="color:rgb(255,0,0)">min=-9999 max=-9999</span><br>
</div>
<br>
<p></p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_8841829022642877653divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>De:</b> Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
<b>Enviado:</b> terça-feira, 23 de janeiro de 2018 19:16:28<br>
<b>Para:</b> Vanúcia Schumacher<br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<b>Assunto:</b> Re: [ncl-talk] shapefiles problem</font>
<div> </div>
</div><div><div class="h5">
<div>
<div dir="ltr">
<div>
<div>
<div>Hi,<br>
<br>
</div>
Looking at your ascii file, I see its either all missing_values, or a handful of zeros. There's nothing to contour in a uniform scalar field, and NCL produced the contour plot I would expect in that case.<br>
<br>
</div>
Hope that helps...<br>
</div>
Rick<br>
<div class="m_8841829022642877653x_gmail_extra"><br>
<div class="m_8841829022642877653x_gmail_quote">On Tue, Jan 23, 2018 at 1:14 PM, Vanúcia Schumacher via ncl-talk
<span dir="ltr"><<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>></span> wrote:<br>
<blockquote class="m_8841829022642877653x_gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div id="m_8841829022642877653x_m_1693159819235043211divtagdefaultwrapper" dir="ltr" style="font-size:14pt;color:#000000;font-family:Calibri,Helvetica,sans-serif">
<p style="margin-top:0;margin-bottom:0"><span id="m_8841829022642877653x_m_1693159819235043211result_box" lang="en"><span>I was able to open the shapefile, but I can't plot the data_mask, because all the values are missing.</span>
<span>I do not know if I'm doing something wrong when using the shapefile.</span></span></p>
<p style="margin-top:0;margin-bottom:0"><br>
<span id="m_8841829022642877653x_m_1693159819235043211result_box" lang="en"><span></span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_8841829022642877653x_m_1693159819235043211result_box" lang="en"><span><span id="m_8841829022642877653x_m_1693159819235043211result_box" class="m_8841829022642877653x_m_1693159819235043211short_text" lang="en"><span>I've attached the files I'm using to facilitate</span></span></span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_8841829022642877653x_m_1693159819235043211result_box" lang="en"><span><br>
</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_8841829022642877653x_m_1693159819235043211result_box" lang="en"><span>I appreciate any help</span></span><br>
</p>
</div>
</div>
<br>
______________________________<wbr>_________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br>
</blockquote>
</div>
<br>
</div>
</div>
</div>
</div></div></div>
</blockquote></div><br></div>