<div dir="ltr">Hi Rick,<div>Thanks for trying.</div><div>I am using the same shapefile_utils.ncl across the two machines.</div><div>As you suggest, I have tried the command in both machines and it looks like x/y coordinates are swapped in maschine #2 as shown below.</div><div><br></div><div><b>machine #1</b></div><div>(0)     -0.1916935366010736 51.34212550232069<br>(1)     -0.1918504471635003 51.34216369759741<br>(2)     -0.1919426195340767 51.3421861363549<br>.</div><div>.</div><div>.</div><div>(25817) -0.1912990774713476 51.34202950288345<br>(25818) -0.1916318559031973 51.34211049447309<br>(25819) -0.1916935366010736 51.34212550232069<br></div><div><br></div><div><b>machine #2</b></div><div>(0)     51.34212549581116 -0.1916934955276336<br>(1)     51.34216369108803 -0.1918504060899556<br>(2)     51.34218612984562 -0.1919425784604706<br></div><div>.</div><div>.</div><div>.</div><div>(25817) 51.34202949637353 -0.1912990363981704<br>(25818) 51.34211048796352 -0.1916318148297985<br>(25819) 51.34212549581116 -0.1916934955276336<br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 5, 2021 at 7:11 PM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Well this is truly bizarre. There's nothing I can think of that would or should account for the lat/lon coordinate swapping between the two machines (unless the shapefile_utils.ncl files are different across the two machines). Things work as expected for me. Unfortunately, without access to something like machine #2 that exhibits the behavior, I personally have no way of drilling down further to find out what the problem is. <br></div><div><br></div><div>Not that this helps solve the problem, but I'm curious if the x/y coordinates are really juxtaposed; on machine #2:</div><div><br></div><div>     f=addfile("file_gor_2011.shp","r")</div><div>     print(f->x + " " + f->y)</div><div><br></div><div>Do those values look like they are swapped?</div><div><br></div><div>I wish I had a better answer for you :-/</div><div><br></div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 5, 2021 at 8:37 AM S Br <<a href="mailto:sbr.climate@gmail.com" target="_blank">sbr.climate@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Rick,</div><div>The below are the details.</div><div><br></div><div>uname -a</div><div>Computer#1</div><div>Linux sbclim 4.15.0-128-generic #131-Ubuntu SMP Wed Dec 9 06:57:35 UTC 2020 x86_64 x86_64 x86_64 GNU/Linux<br></div><div>Computer#2<br></div><div>Linux sbclim 3.10.0-1127.19.1.el7.x86_64 #1 SMP Tue Aug 25 17:23:54 UTC 2020 x86_64 x86_64 x86_64 GNU/Linux<br></div><div><br></div><div>ls -l file_gor_2011.shp<br></div><div>Computer#1<br></div><div>-rw-r--r-- 1 climdir env 413276 Dec  8 16:11 file_gor_2011.shp<br></div><div>Computer#2<br></div><div>-rw-r--r-- 1 climdir users 413276 Jan  5 10:54 file_gor_2011.shp<br></div><div><br></div><div>Please find attached here the shape files.</div><div><br></div><div>Thanks</div><div>SB</div><div><br></div><div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 5, 2021 at 2:37 PM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi,</div><div><br></div><div>That result is very strange indeed. What are computers 1 and 2, i.e., what is the output of:</div><div><br></div><div>   uname -a</div><div><br></div><div>And just to be certain, what is the output of this on each computer:</div><div><br></div><div>  cd ...to-directory-where-file_gor_2011.shp-resides....</div><div>  ls -l file_gor_2011.shp</div><div><br></div><div>Finally, can you send me file_gor_2011.shp (recall that means the .sph, .shx, .dbf files)?</div><div><br></div><div>Rick<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jan 5, 2021 at 5:53 AM S Br via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Dave,</div><div>Thanks for your reply. </div><div>I have now compared the diagnostic results between the two computers and the below are the results.</div><div><b>Computer#1:</b></div><div>(0)      ======================================================================<br>(0)     Filename: "file_gor_2011.shp"<br>(0)       Geometry type: polygon<br>(0)    # of features: 1<br>(0)          Min/max lat:     51.29/  51.69<br>(0)         Min/max lon:     -0.51/   0.33<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)              label : string<br>(0)          name : string<br>(0)           code : string<br>(0)    ======================================================================<br></div><div><b>Computer#2:</b><br></div><div>(0)       ======================================================================<br>(0)     Filename: "file_gor_2011.shp"<br>(0)       Geometry type: polygon<br>(0)    # of features: 1<br>(0)          Min/max lat:     -0.51/   0.33<br>(0)         Min/max lon:     51.29/  51.69<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)              label : string<br>(0)          name : string<br>(0)           code : string<br>(0)    ======================================================================<br></div><div><br></div><div>I could see in one computer (computer#2 above) that the latitude values are replacing the longitude values and vice versa, and that may be the reason the country outlines are not displaying here.</div><div><br></div><div>Could you please suggest how to solve this issue? I am using the same data, scripts and only difference is the computer.</div><div><br></div><div>Thank you.</div><div>SB</div><div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 9, 2020 at 5:30 PM Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div>It would be better to apply debugging techniques to narrow down the problem, rather than to speculate as to the cause.  Use <b>addfile</b> to see if NCL can read the shapefile.  Use diagnostic print statements to examine the shapefile structure in memory.  Use <b>gsn_add_polygon</b> to see if NCL can draw a simple outline.</div><div><br></div><div>To understand the structure, refer to the shapefile documentation at the top of the NCL shapefiles example page.  That mentions a special procedure <b>print_shapefile_info</b> which will print additional useful details.  Compare diagnostic results between your two computers.  Good luck.</div><div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 9, 2020 at 1:16 AM S Br via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi All,<div>I am trying to plot a figure using the shapefile. Everything is working fine with the NCL script but it doesn't draw the country outline given as a shapefiles. However, the same script works fine in another computer and draws the required country outline using the shapefile.<br><div>What could be the possible reason that the NCL script works fine but fails to draw the country outline using shapefiles without showing any error?</div><div><br></div><div>Both the computers are installed with the NCL 6.6.2.<br></div><div><br></div><div>Thank you.</div><div>SB</div></div></div></blockquote></div></div></div>
</blockquote></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div></div>
</blockquote></div>
</blockquote></div>