<div dir="ltr"><div>Hi Mary, <br></div><div><br></div><div>I took your advice again and checked the shapefile with
print_shapefile_info
and that was the problem. The shapefile were on different coordinates, my bad. <br></div><div><br></div><div>Thank you very much and so to Rick that helped me a lot understanding the function. <br></div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Oct 2, 2018 at 3:16 PM Mary Haley <<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">I forgot to mention : please respond back to ncl-talk with all follow-up messages.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thanks,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 2, 2018 at 2:15 PM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</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 dir="ltr"><div class="gmail_default" style="font-size:small">The way shapefile_mask_data works is that it makes a copy of your original variable, and then it sets values in that new variable to missing at locations where the lat/lon location falls outside the mask area of interest. So, it is correct that what this function gives back is a variable of the same size.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I think what's going on is that it didn't find any lat/lon pairs from your original data that fell inside the lat/lon area defined in your shapefile. This might happen because your two areas do not overlap, or if your shapefile longitudes go from 0 to 360 and your data longitudes go from -180 to 180 (or vice versa).</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I would first suggest that you make sure your shapefile area and your data lat/lon area actually intersect. You can do this by using the "print_shapefile_info" procedure mentioned at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/shapefiles.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/shapefiles.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">and then simply using printMinMax to print your data's lat/lon min and max:</div><div class="gmail_default"><br></div><div class="gmail_default">printMinMax(var@lat2d,0)</div><div class="gmail_default">printminMax(var@lon2d,0)</div><div class="gmail_default"><br></div><div class="gmail_default">If the lat/lon areas don't intersect, then you will get back all missing values.</div><span class="m_2998005738679421870HOEnZb"><font color="#888888"><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default"><br></div><div class="gmail_default"><br></div></font></span></div></div><div class="m_2998005738679421870HOEnZb"><div class="m_2998005738679421870h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 2, 2018 at 11:00 AM, Christian Adrián Álvarez Báez <span dir="ltr"><<a href="mailto:christian.alvarez813@gmail.com" target="_blank">christian.alvarez813@gmail.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 dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hi Mary, <br></div><div><br></div><div>I took your advice and that made the code advance. I have tried before with setting "opt" on
<div style="color:rgb(212,212,212);background-color:rgb(30,30,30);font-family:Consolas,"Courier New",monospace;font-weight:normal;font-size:14px;line-height:19px;white-space:pre-wrap"><div><span style="color:rgb(212,212,212)"> data_keep = shapefile_mask_data(var,shp_filename,opt)</span></div></div>
</div><div><br></div><div>but that sent me and error: <br></div><div>---- shapefile_mask_data : Error: The given variable name to use does not exist on the given shapefile.</div><div>My only options are:</div><div>
<div style="color:rgb(212,212,212);background-color:rgb(30,30,30);font-family:Consolas,"Courier New",monospace;font-weight:normal;font-size:14px;line-height:19px;white-space:pre-wrap"><div><span style="color:rgb(106,153,85)">;---Set masking options for shapefile_mask_data function</span></div><div><span style="color:rgb(212,212,212)"> opt = </span><span style="color:rgb(86,156,214)">True</span></div><div><span style="color:rgb(212,212,212)"> opt@shape_var = </span><span style="color:rgb(206,145,120)">"NAME_1"</span></div><div><span style="color:rgb(212,212,212)"> opt@shape_names = </span><span style="color:rgb(206,145,120)">"GOLFO"</span></div><div><span style="color:rgb(212,212,212)"> opt@debug = </span><span style="color:rgb(86,156,214)">True</span></div><div><span style="color:rgb(212,212,212)"> opt@keep = </span><span style="color:rgb(86,156,214)">True</span></div><div><span style="color:rgb(212,212,212)"> data_keep = shapefile_mask_data(var,shp_filename,opt)</span></div></div>
</div><div><br></div><div>So I tried with just "True" and it worked, but I don't think that's ok because when I print min and max of the temperature on the data_keep, it says:<br>(0) Temperature (C) : min=9.96921e+36 max=9.96921e+36<br></div><div><br></div><div>A print min and max over var would send me:</div><div>(0) Temperature (C) : min=1.26712 max=20.9307</div><div><br></div><div>When I print a var summary on data_keep it let me know that shapefile_mask_data is not working (al least not how I expected) since it has the same size as the variable "var" before applying the mask even though the values are different as I said above:</div><div>---------------------------------<br></div><div>Variable: var<br>Type: float<br>Total Size: 47124 bytes<br> 11781 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [south_north | 119] x [west_east | 99]<br>Coordinates:<br>Number Of Attributes: 8<br> description : Temperature<br> coordinates : XLONG XLAT XTIME<br> stagger :<br> units : C<br> MemoryOrder : XYZ<br> FieldType : 104<br> lon2d : <ARRAY of 11781 elements><br> lat2d : <ARRAY of 11781 elements><br>--------------------------------<br>Variable: data_keep<br>Type: float<br>Total Size: 47124 bytes<br> 11781 values<br>Number of Dimensions: 2<br>Dimensions and sizes: [south_north | 119] x [west_east | 99]<br>Coordinates:<br>Number Of Attributes: 9<br> description : Temperature<br> coordinates : XLONG XLAT XTIME<br> stagger :<br> units : C<br> MemoryOrder : XYZ<br> FieldType : 104<br> lon2d : <ARRAY of 11781 elements><br> lat2d : <ARRAY of 11781 elements></div><div><br></div><div>I think it has to do with the not using "opt" instead of just "True" but then I need to solve that "shapefile_mask_data : Error: The given variable name to use does not exist on the given shapefile" error first. Sorry for so much detail, but I just wanted to give the pieces of the puzzle.</div><div><br></div><div>Thank you for your attention!<br></div></div></div></div></div></div></div><div class="m_2998005738679421870m_-8798404079134532226HOEnZb"><div class="m_2998005738679421870m_-8798404079134532226h5"><br><div class="gmail_quote"><div dir="ltr">On Mon, Oct 1, 2018 at 3:06 PM Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:small">If I may jump in. This error:</div><div class="gmail_default" style="font-size:small"><br></div></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div dir="ltr"><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">fatal:Argument type mismatch on argument (1) of (shapefile_mask_data) can not coerce</span></div></div></blockquote><div dir="ltr"><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">is indicating that the second argument (argument count starts at 0) is not the correct expected type.</span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">This argument is expecting the shapefile that's already been opened by addfile, and not the name of the shapefile itself.</span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div class="gmail_default"><span style="font-size:12.8px">For example:</span></div><div class="gmail_default"><span style="font-size:12.8px"><br></span></div><div class="gmail_default"><div class="gmail_default"><span style="font-size:12.8px"><font face="monospace, monospace"> shp_file = "ne_10m_geography_regions_polys.shp"</font></span></div><div class="gmail_default"><span style="font-size:12.8px"><font face="monospace, monospace"> a = addfile(shp_file,"r")</font></span></div><div class="gmail_default"><span style="font-size:12.8px"><font face="monospace, monospace"> opt = True</font></span></div><div class="gmail_default"><font face="monospace, monospace"><span style="font-size:12.8px"><i> <...set possible "opt" options here...></i></span></font></div><div class="gmail_default"><span style="font-size:12.8px"><font face="monospace, monospace"> data_mask = shapefile_mask_data(data,shp_file,opt)</font></span></div><div class="gmail_default"><font face="monospace, monospace"><span style="font-size:12.8px"><br></span></font>--Mary</div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Sep 28, 2018 at 1:38 PM, Rick Brownrigg <span dir="ltr"><<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</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>Hmmm...that all looks about right. Have you done a printVarSummary(var) to make sure that variable is what you think it is?</div><div><br></div><div>Rick<br></div></div><br><div class="gmail_quote"><div><div class="m_2998005738679421870m_-8798404079134532226m_623612568931246892m_-8709309445418978871h5"><div dir="ltr">On Fri, Sep 28, 2018 at 1:01 PM Christian Adrián Álvarez Báez <<a href="mailto:christian.alvarez813@gmail.com" target="_blank">christian.alvarez813@gmail.com</a>> wrote:<br></div></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_2998005738679421870m_-8798404079134532226m_623612568931246892m_-8709309445418978871h5"><div dir="ltr"><div dir="ltr"><div dir="ltr">Hi, I'm triying to extract data with a mask from an WRF output using the example shapefiles_19.ncl. I had many errors so looking I found some modifications, including the shapefile_mask_data: Error: not a valid rectilinear, curvilinear, or unstructured grid that I think it was solved by using the copy_VarMeta function. But I'm stuck with the error:</div><div dir="ltr">fatal:Argument type mismatch on argument (1) of (shapefile_mask_data) can not coerce</div><div>caused on line 60 of my script. <br></div><div>
<div style="color:rgb(212,212,212);background-color:rgb(30,30,30);font-family:Consolas,"Courier New",monospace;font-weight:normal;font-size:14px;line-height:19px;white-space:pre-wrap"><div><span style="color:rgb(212,212,212)">tc = </span><span style="color:rgb(220,220,170)">wrf_user_getvar</span><span style="color:rgb(212,212,212)">(a,</span><span style="color:rgb(206,145,120)">"tc"</span><span style="color:rgb(212,212,212)">,it) </span><span style="color:rgb(106,153,85)"></span><span style="color:rgb(212,212,212)"> </span></div><div><span style="color:rgb(212,212,212)">var=tc(</span><span style="color:rgb(181,206,168)">0</span><span style="color:rgb(212,212,212)">,:,:)</span></div><div><span style="color:rgb(220,220,170)">copy_VarMeta</span><span style="color:rgb(212,212,212)">(tc,var)</span></div></div>
</div><div><br></div><div>Line 60<br></div><div>
<div style="color:rgb(212,212,212);background-color:rgb(30,30,30);font-family:Consolas,"Courier New",monospace;font-weight:normal;font-size:14px;line-height:19px;white-space:pre-wrap"><div><span style="color:rgb(212,212,212)"> data_keep = shapefile_mask_data(var,shp_filename,opt)</span></div><div><span style="color:rgb(212,212,212)"></span></div></div>
</div><div>
<div style="color:rgb(212,212,212);background-color:rgb(30,30,30);font-family:Consolas,"Courier New",monospace;font-weight:normal;font-size:14px;line-height:19px;white-space:pre-wrap"><div><span style="color:rgb(212,212,212)"><br></span></div><div><span style="color:rgb(212,212,212)"></span></div></div>
</div><div><br></div><div> The script it's a little messy because of many modifications, I hope I made sense. I don't know if the problem could be the shape, it's just a simple polygon test projected on Lambert conformal conical as the WRF output. <br></div><div><br></div><div>I'm sending my script and my shape files if there is a clue.</div><div><br></div><div>Thank you very much,</div><div>Christian<br></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><br>
</blockquote></div>
<br>_______________________________________________<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><br>
<br></blockquote></div><br></div>
</blockquote></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</blockquote></div>