<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Vanúcia,<br></div><div dir="ltr"><br></div><div>Example shapefiles_21 does not work the way I expected. It needs a shapefile with areas surrounding the region of interest, because it works by blanking out unwanted areas with white, leaving only the central area visible. NCL does not currently have the ability to show contour plot areas inside specified polygons. It can only blank out unwanted polygons in graphics mode.</div><div><br></div><div>The shapefile in website example shapefiles_21 contains surrounding areas, as seen in figure 1 below. When you specify the list of area names to display, what actually happens is that all of the areas NOT in the list are blanked out in white. Please notice the difference between this shapefile, and your original shapefile that contains only your wanted areas.</div><br>Your original program crashed because in "areas_of_interest" you listed ALL AREAS inside the shape file. This reduced the number of blanking areas to zero. This resulted in an unexpected zero length array which was not properly handled by the program. There is nothing basically wrong with the shapefile, it's just that example shapefiles_21 can't use it in this manner.</div><div dir="ltr"><br></div><div dir="ltr">To demonstrate, I removed one of the four areas from your area list. This avoided the program crash. The modified script with a few other changes is attached. In figure 2, notice two things. (1) The area Cuenca Coton is blanked to white, because I removed it from the area list. (2) All of the unwanted surrounding areas are VISIBLE in the contour plot, because they are NOT available in the shapefile; therefore they can not be blanked out. I hope this shows you how example 21 actually works.</div><div dir="ltr"><br></div><div dir="ltr">(A) Your best solution would be to find a larger shapefile which contains many basin outlines for this part of Chile. That should just drop in to the current script, and blank out the surrounding unwanted areas.</div><div dir="ltr"><br></div><div dir="ltr">(B) It is possible to construct exclusion masks using shapefile editing tools. Draw two rectangles which completely cover your plot frame, then use your original shapefile to cut out the desired visible areas, leaving only exclusion areas. You need two rectangles, not one, because NCL does not currently handle polygons with internal holes. There are many shapefile tools, but I don't have experience with them. Perhaps someone else could recommend a tool that can do this.<br></div><div dir="ltr"><br></div><div>(C) It is even possible to use NCL to manually extract polyline segments from your original shapefile, and combine them with rectangle vertices to construct these two exclusion masks. This is tricky and tedious, so the first two solutions are much better if you can do them. Good luck.</div><div><br></div><div>--Dave</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 31, 2019 at 5:58 AM Vanúcia Schumacher <<a href="mailto:vanucia-schumacher@hotmail.com">vanucia-schumacher@hotmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
Unfortunately I do not have another shapefile, just this one referring to my work.<br>
</div>
<div id="gmail-m_-4449798872923559840appendonsend"></div>
<br><b style="font-size:11pt;color:rgb(0,0,0);font-family:Calibri,sans-serif">De:</b><span style="font-size:11pt;color:rgb(0,0,0);font-family:Calibri,sans-serif"> Dave Allured - NOAA Affiliate <</span><a href="mailto:dave.allured@noaa.gov" target="_blank" style="font-size:11pt;font-family:Calibri,sans-serif">dave.allured@noaa.gov</a><span style="font-size:11pt;color:rgb(0,0,0);font-family:Calibri,sans-serif">></span><div id="gmail-m_-4449798872923559840divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000">
<b>Enviado:</b> terça-feira, 29 de janeiro de 2019 16:16<br></font><div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">Vanúcia,</div>
<div dir="ltr"><br>
</div>
<div>If you can find another shape file of the same region from a different source, give that a try. This might be quicker than waiting for analysis of the original shape file. If you find a working shape file, please let us know.</div>
<div><br>
</div>
<div>--Dave</div>
<div> <br>
</div>
<br>
<div class="gmail-m_-4449798872923559840x_gmail_quote">
<div dir="ltr" class="gmail-m_-4449798872923559840x_gmail_attr">On Mon, Jan 28, 2019 at 10:50 AM Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>> wrote:<br>
</div>
<blockquote class="gmail-m_-4449798872923559840x_gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div dir="ltr">Thanks. I can't look at this today. I will get back to you in a day or two. Perhaps someone else would like to check out the shape file and the script.</div>
<div dir="ltr">
<div><br>
</div>
</div>
<br>
<div class="gmail-m_-4449798872923559840x_gmail_quote">
<div dir="ltr" class="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail_attr">On Mon, Jan 28, 2019 at 10:45 AM Vanúcia Schumacher <<a href="mailto:vanucia-schumacher@hotmail.com" target="_blank">vanucia-schumacher@hotmail.com</a>> wrote:<br>
</div>
<blockquote class="gmail-m_-4449798872923559840x_gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span>Is attached.<br>
</span>
<div><br>
</div>
<span>Thank you for the advance</span><br>
</div>
<div id="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711appendonsend"></div>
<br><div id="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711divRplyFwdMsg" dir="ltr">
<font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>De:</b> Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>><br>
<b>Enviado:</b> segunda-feira, 28 de janeiro de 2019 15:37<br></font><div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">Vanúcia,</div>
<div dir="ltr"><br>
</div>
<div>Shape files come in different varieties. It looks like example shapefiles_21 may need an upgrade to be able to read your shape file. If you can send a copy of the shape file, and also your current script shape21.ncl, I will test it on my system. Please
do not include the scientific data file; I do not need it for testing the shape file. The NCL support team might need to fix this.</div>
<div><br>
</div>
<div>--Dave</div>
<div><br>
</div>
<br>
<div class="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711x_gmail_quote">
<div dir="ltr" class="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711x_gmail_attr">
On Mon, Jan 28, 2019 at 9:47 AM Vanúcia Schumacher <<a href="mailto:vanucia-schumacher@hotmail.com" target="_blank">vanucia-schumacher@hotmail.com</a>> wrote:<br>
</div>
<blockquote class="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711x_gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span style="font-size:12pt">Thanks for the tip.</span><span><br>
</span><span style="font-size:12pt">I'm trying to use the script from example</span><span style="font-size:12pt">
</span><span><span style="color:rgb(33,33,33);font-size:12pt;background-color:rgb(255,255,255);display:inline">shapefiles_21</span></span><span style="font-size:12pt">, b</span><span style="font-size:12pt">ut I'm having this error in function. I can not
fix it.</span><br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span><br>
</span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span>
<div style="margin:0px;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<span style="font-size:12pt">fatal:Subscript out of range, error in subscript #0</span><br>
</div>
<div style="margin:0px;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<span style="font-size:12pt">fatal:An error occurred reading features</span><br>
</div>
<div style="margin:0px;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<span style="font-size:12pt">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 56 in file shape21.ncl</span></div>
<div style="margin:0px;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<span style="font-size:12pt">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 250 in file shape21.ncl</span><br>
</div>
<br>
</span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span style="font-size:12pt">The lines correspond to: </span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span><span style="font-size:12pt">function get_areas_of_interest(shp_file_name,shp_var_name,opt[1]:logical)</span><span><br>
</span>
<div><span style="font-size:12pt">begin</span><br>
</div>
<div><span style="font-size:12pt">;---Open the shapefile</span><br>
</div>
<div><span style="font-size:12pt"> f = addfile(shp_file_name,"r")</span><br>
</div>
<div><span style="font-size:12pt"> features = f->$shp_var_name$</span><br>
</div>
<div><br>
</div>
<div><span style="font-size:12pt"> if(opt.and.isatt(opt,"areas_to_exclude")) then</span><br>
</div>
<div><span style="font-size:12pt"> features@_FillValue = default_fillvalue(typeof(features))</span><br>
</div>
<div><span style="font-size:12pt"> do na=0,dimsizes(opt@areas_to_exclude)-1</span><br>
</div>
<div><span style="font-size:12pt"> ii := ind(features.eq.opt@areas_to_exclude(na))</span><br>
</div>
<div><span style="font-size:12pt"> if(.not.any(ismissing(ii))) then</span><br>
</div>
<div><span style="font-size:12pt"> features(ii) = features@_FillValue</span><br>
</div>
<div><span style="font-size:12pt"> end if</span><br>
</div>
<div><span style="font-size:12pt"> end do</span><br>
</div>
<div><span style="font-size:12pt"> return(features(ind(.not.ismissing(features)))) ------> </span><span style="color:rgb(87,6,6);font-size:14pt"> line 56</span><br>
</div>
<div><span style="font-size:12pt"> else</span><br>
</div>
<div><span style="font-size:12pt"> return(features)</span><br>
</div>
<div><span style="font-size:12pt"> end if</span><br>
</div>
<span style="font-size:12pt">end</span><br>
</span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span><span><br>
</span></span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span style="font-size:12pt">opt = True</span><span><br>
</span>
<div><span style="font-size:12pt"> opt@areas_to_exclude = areas_of_interest</span><br>
</div>
<div><span style="font-size:12pt"> areas_to_fill = get_areas_of_interest(shp_filename1,shape_var_name,opt) ----></span><span style="font-size:14pt">
</span><span style="color:rgb(87,6,6);font-size:14pt">line 250</span><br>
</div>
<span></span><br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span style="font-size:12pt">Infos:</span><br>
</div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span><br>
</span></div>
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<span><span style="font-size:12pt">0) Terrain Height (m) : min=0 max=5070.39</span><span><br>
</span>
<div><span style="font-size:12pt">(0) ======================================================================</span><br>
</div>
<div><span style="font-size:12pt">(0) Filename: "cuencasx.shp"</span><br>
</div>
<div><span style="font-size:12pt">(0) Geometry type: polygon</span><br>
</div>
<div><span style="font-size:12pt">(0) # of features: 4</span><br>
</div>
<div><span style="font-size:12pt">(0) Min/max lat: -34.77/ -34.23</span><br>
</div>
<div><span style="font-size:12pt">(0) Min/max lon: -70.53/ -70.13</span><br>
</div>
<div><span style="font-size:12pt">(0) Variable names and their types:</span><br>
</div>
<div><span style="font-size:12pt">(0) geometry : integer</span><br>
</div>
<div><span style="font-size:12pt">(0) segments : integer</span><br>
</div>
<div><span style="font-size:12pt">(0) x : double</span><br>
</div>
<div><span style="font-size:12pt">(0) y : double</span><br>
</div>
<div><span style="font-size:12pt">(0) COD_CUEN : string</span><br>
</div>
<div><span style="font-size:12pt">(0) COD_SUBC : string</span><br>
</div>
<div><span style="font-size:12pt">(0) COD_SSUBC : string</span><br>
</div>
<div><span style="font-size:12pt">(0) NOMBRE : string</span><br>
</div>
<div><span style="font-size:12pt">(0) AREAKM2 : double</span><br>
</div>
<div><span style="font-size:12pt">(0) ======================================================================</span><br>
</div>
<div><br>
</div>
<div><span style="font-size:12pt">Variable: features</span><br>
</div>
<div><span style="font-size:12pt">Type: string</span><br>
</div>
<div><span style="font-size:12pt">Total Size: 32 bytes</span><br>
</div>
<div><span style="font-size:12pt"> 4 values</span><br>
</div>
<div><span style="font-size:12pt">Number of Dimensions: 1</span><br>
</div>
<div><span style="font-size:12pt">Dimensions and sizes: [num_features | 4]</span><br>
</div>
<div><span style="font-size:12pt">Coordinates: </span><br>
</div>
<div><span style="font-size:12pt">Number Of Attributes: 0</span><br>
</div>
<div><span style="font-size:12pt">(0) Cuenca Cipreses</span><br>
</div>
<div><span style="font-size:12pt">(1) Cuenca Coton</span><br>
</div>
<div><span style="font-size:12pt">(2) Cuenca Cortaderal</span><br>
</div>
<div><span style="font-size:12pt">(3) Cuenca Universidad</span><br>
</div>
<div><span style="font-size:12pt">(0) ==================================================</span><br>
</div>
<div><span style="font-size:12pt">(0) Shapefile: cuencasx.shp</span><br>
</div>
<div><span style="font-size:12pt">(0) Areas of interest: Cuenca Cipreses,Cuenca Coton,Cuenca Cortaderal,Cuenca Universidad</span><br>
</div>
<div><span style="font-size:12pt">(0) min_lat_chk: -37.4943</span><br>
</div>
<div><span style="font-size:12pt">(0) max_lat_chk: -31.9474</span><br>
</div>
<div><span style="font-size:12pt">(0) min_lon_chk: -73.6526</span><br>
</div>
<div><span style="font-size:12pt">(0) max_lon_chk: -66.7254</span><br>
</div>
<div><span style="font-size:12pt">(0) min_lat_data: -37.4943</span><br>
</div>
<div><span style="font-size:12pt">(0) max_lat_data: -31.9474</span><br>
</div>
<div><span style="font-size:12pt">(0) min_lon_data: -73.6526</span><br>
</div>
<div><span style="font-size:12pt">(0) max_lon_data: -66.7254</span><br>
</div>
<div><span style="font-size:12pt">(0) 4761 data values originally</span><br>
</div>
<div><span style="font-size:12pt">(0) Will keep data values inside given shapefile areas</span><br>
</div>
<div><span style="font-size:12pt">(0) ==================================================</span><br>
</div>
<div><span style="font-size:12pt">(0) 13 data values kept</span><br>
</div>
<div><span style="font-size:12pt">(0) shapefile_mask_data: elapsed time: 0.00758901 CPU seconds.</span><br>
</div>
<div><span style="font-size:12pt">(0) ==================================================</span><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><b style="font-size:11pt;font-family:Calibri,sans-serif">De:</b><span style="font-size:11pt;font-family:Calibri,sans-serif"> Dave Allured - NOAA Affiliate <</span><a href="mailto:dave.allured@noaa.gov" style="font-size:11pt;font-family:Calibri,sans-serif" target="_blank">dave.allured@noaa.gov</a><span style="font-size:11pt;font-family:Calibri,sans-serif">></span><br>
</div>
</span></div>
<div id="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711x_gmail-m_1098706355727375035divRplyFwdMsg" dir="ltr">
<font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>Enviado:</b> domingo, 27 de janeiro de 2019 18:23<br></font><div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">Vanúcia,</div>
<div dir="ltr"><br>
</div>
<div>Your two data files are on different grids. If you only need to make plots that show the same region on a map, then use graphics masking, not data masking. See the third plot in example shapefiles_21 on this page:</div>
<div><br>
</div>
<div><a href="https://www.ncl.ucar.edu/Applications/shapefiles.shtml" target="_blank">https://www.ncl.ucar.edu/Applications/shapefiles.shtml</a><br>
</div>
<div><br>
</div>
<div>If you need to perform calculations between the two data sets, then you will need to regrid the data from one file to match the other grid. Please see this documentation for regridding:</div>
<div><br>
</div>
<div><a href="https://www.ncl.ucar.edu/Applications/regrid.shtml" target="_blank">https://www.ncl.ucar.edu/Applications/regrid.shtml</a><br>
</div>
<div><a href="https://www.ncl.ucar.edu/Applications/ESMF.shtml" target="_blank">https://www.ncl.ucar.edu/Applications/ESMF.shtml</a><br>
</div>
<div><br>
</div>
<div>--Dave</div>
<div><br>
</div>
<br>
<div class="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711x_gmail-m_1098706355727375035x_gmail_quote">
<div dir="ltr" class="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711x_gmail-m_1098706355727375035x_gmail_attr">
On Sun, Jan 27, 2019 at 12:40 PM Vanúcia Schumacher <<a href="mailto:vanucia-schumacher@hotmail.com" target="_blank">vanucia-schumacher@hotmail.com</a>> wrote:<br>
</div>
<blockquote class="gmail-m_-4449798872923559840x_gmail-m_-1350803966612980561gmail-m_4020804067210245711x_gmail-m_1098706355727375035x_gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div style="font-family:Calibri,Helvetica,sans-serif;font-size:14pt;color:rgb(0,0,0)">
<div style="margin:0px;font-size:14pt;font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont;background-color:rgb(255,255,255)">
Hi all,</div>
<div style="margin:0px;font-size:14pt;font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:14pt;font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont;background-color:rgb(255,255,255)">
I need support to find the problem with my script, <span style="margin:0px">with the purpose of cutting two different data in a shapefile, </span><span style="margin:0px">but keeping the information such as time lat and<span> </span></span>lon.</div>
<div style="margin:0px;font-size:14pt;font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont;background-color:rgb(255,255,255)">
My script (attachment) is "cutting" the region from shapefile different for each input (data) (see Figure attachment).<br>
</div>
<div style="margin:0px;font-size:14pt;font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont;background-color:rgb(255,255,255)">
I'd like to select the shapefile region of the same size independent of the input file, and keep the time, lat, and lon information.<br>
</div>
<div style="margin:0px;font-size:14pt;font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:14pt;font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont;background-color:rgb(255,255,255)">
I appreciate any help</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</blockquote></div></div></div>