<div dir="ltr"><div class="gmail_default" style="font-size:small">Jagan,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I don't recognize the script that you sent, so it may be an older version.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You may want to look at the newer "shapefile_mask_data.ncl" script, which allows you to mask various types of data against a shapefile outline.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">To see an example of using this function, go to:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style><a href="http://www.ncl.ucar.edu/Applications/shapefiles.shtml">http://www.ncl.ucar.edu/Applications/shapefiles.shtml</a><br></div><div class="gmail_default" style><br></div><div class="gmail_default" style>and look at examples 4, 11, 12, 14, or 16.</div><div class="gmail_default" style><br></div><div class="gmail_default" style>--Mary</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 12, 2014 at 6:04 PM, jagan TNAU <span dir="ltr"><<a href="mailto:jagan@tnau.ac.in" target="_blank">jagan@tnau.ac.in</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>Thanks Rick.<br></div><div><br></div><div>The entire script was provided to me by Mary. I have used it for producing 2x4 panel plot. Now I have planned to go for 3x4 panel plots and modified it. The script reads a shape file and select a particular region (districts) for plotting. I have not changed anything related to the shape file reading. </div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 12, 2014 at 10:08 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><div>Hi,<br><br></div>I don't fully understand your script, but at line 249, you are looping over the number of lat/lon values in the shapefile:<br><br> do iltn=0,nlatlon_shp-1<br> ;---gc_inout returns True if point is inside lat/lon area.<br> if(ismissing(tmax1_mask(iltn))) then<br><br></div><div>and indexing into tmax1_mask based upon "iltn". But tmax1_mask was dimensioned as:<br></div><div><br> values = asciiread("tmax-had45.txt",-1,"float")<br>...<br> ncols = 15<br>...<br> tmax1 = values(2::ncols)<br>...<br> tmax1_mask = new(dimsizes(tmax1),typeof(tmax1),tmsg)<br><br></div><div>So its dimensions would seem to be unrelated to the number of points in the shapefile (?) <br><br></div><div>In any case, the error message means what it says -- an attempt to access an array outside its bounds.<br><br></div><div>Hope that helps pinpoint your problem...<br></div><div>Rick<br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, Nov 12, 2014 at 8:22 AM, jagan TNAU <span dir="ltr"><<a href="mailto:jagan@tnau.ac.in" target="_blank">jagan@tnau.ac.in</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div>Hello,</div><div><br></div>Please find attached a shape file script for which I am receiving the error below. Kindly help me to correct it.<div>---------------</div><div><div> NCAR Command Language Version 6.0.0</div><div> The use of this software is governed by a License Agreement.</div><div> See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> for more details.</div><div>(0)<span style="white-space:pre-wrap"> </span></div><div>(0)<span style="white-space:pre-wrap"> </span>numFeatures = 132</div><div>(0)<span style="white-space:pre-wrap"> </span>numSegments = 202</div><div>(0)<span style="white-space:pre-wrap"> </span>size of lat/lon arrays in shapefile = 192128</div><div>(0)<span style="white-space:pre-wrap"> </span>minlat_reg = 12.2281</div><div>(0)<span style="white-space:pre-wrap"> </span>maxlat_reg = 13.5628</div><div>(0)<span style="white-space:pre-wrap"> </span>minlon_reg = 79.2832</div><div>(0)<span style="white-space:pre-wrap"> </span>maxlon_reg = 80.3488</div><div>(0)<span style="white-space:pre-wrap"> </span>Number of line segments to draw = 2</div><div>(0)<span style="white-space:pre-wrap"> </span>Largest line segment length = 8691</div><div>(0)<span style="white-space:pre-wrap"> </span>Number of areas to mask = 200</div><div>(0)<span style="white-space:pre-wrap"> </span>Larkest mask segment length = 10255</div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:An error occurred reading tmax1_mask</div><div>fatal:Execute: Error occurred at or near line 249 in file tmax-tn-mask.ncl</div><div>--------------------</div><span><font color="#888888"><div><br></div>-- <br><div>With regards<br><br>Dr.R.Jagannathan <br>Professor of Agronomy<br>Agro Climate Research Centre<br>Tamil Nadu Agricultural University <br>Coimbatore - 641 003 India <br><br>PHONE: Mob: <a href="tel:%2B91%2094438%2089891" value="+919443889891" target="_blank">+91 94438 89891</a><br><br>DO NOT PRINT THIS E-MAIL UNLESS NECESSARY. THE ENVIRONMENT CONCERNS US ALL.</div>
</font></span></div></div>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div>With regards<br><br>Dr.R.Jagannathan <br>Professor of Agronomy<br>Agro Climate Research Centre<br>Tamil Nadu Agricultural University <br>Coimbatore - 641 003 India <br><br>PHONE: Mob: <a href="tel:%2B91%2094438%2089891" value="+919443889891" target="_blank">+91 94438 89891</a><br><br>DO NOT PRINT THIS E-MAIL UNLESS NECESSARY. THE ENVIRONMENT CONCERNS US ALL.</div>
</div>
</div></div><br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>