<div dir="ltr"><div>var.txt: The 1st column of numbers is the line number [added for clarity].<br></div><div> ln lat lon<br></div><div> 1 <b>10.0 </b> <b>78.0 </b> 12.3 28.4 15.1 0.4 29.3 78.2 <===lat ; longitudes starts<br></div><div> 2 10.0 78.2 16.9 29.5 15.3 0.1 30.2 81.0 <br></div><div>...</div><div>...</div><div> 10 10.0 79.8 19.5 29.4 21.7 0.2 32.8 85.5<br> 11 <b>10.0 80.0 </b> 20.0 28.7 21.9 0.2 34.5 85.5<br> 12 <b>10.2</b> <b>78.0 </b> 12.5 27.9 13.9 0.3 25.7 84.7 <== new lat row; lon restarts</div><div> 13 10.2 78.2 15.3 28.9 14.2 0.2 27.1 82.9 <br></div><div> 14 10.2 78.4 15.3 29.6 14.9 0.4 27.5 83.4<br> 15 10.2 78.6 15.1 29.9 15.6 0.6 27.0 83.2<br> 16 10.2 78.8 16.0 30.3 16.5 0.6 26.3 84.3<br> 17 10.2 79.0 17.6 30.7 17.6 0.5 26.0 88.3<br> 18 10.2 79.2 18.3 30.9 18.5 0.2 26.0 92.8<br> 19 10.2 79.4 18.8 30.8 20.0 0.2 26.6 89.3<br> 20 10.2 79.6 19.4 30.5 21.5 0.2 28.5 85.1<br> 21 10.2 79.8 18.9 30.3 21.5 0.2 29.8 85.0<br> <b>22 10.2 80.0 </b> 19.4 29.4 21.7 0.2 32.1 85.2<br> <b>23 10.4</b> <b>78.0</b> 12.3 27.6 13.7 0.2 19.3 87.6 <== new lat row; lon restarts</div><div>...</div><div> <b>33 10.4 </b> 80.0 18.9 30.1 21.5 0.2 29.5 84.8<br> <b>34 10.6 </b> <b>78.0 </b> 14.0 28.7 14.2 0.2 14.4 85.9 <== new lat row; lon restarts</div><div>...</div><div><b>120 12.0 </b> 79.8 14.8 31.0 18.9 0.2 22.6 85.1<br><b>121 12.0</b> 80.0 17.7 30.6 20.3 0.1 22.4 86.3</div><div><br></div><div>It becomes amatter of book-keeping to plave the values onto a grid.</div><div>A rectilinear grid is *necessary* for <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/linint2_points.shtml" target="_blank"><b>linint2_points</b></a> or <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/linint2_points_Wrap.shtml" target="_blank"><b>linint2_points_Wrap <br></b></a></div><div><br></div><div>See attached. There are many print statements. Please look at the output.</div><div><br></div><div>%> ncl <b>jagan_TNAU.linint2.ncl</b></div><div><br></div><div>Good luck</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Aug 17, 2019 at 1:16 PM Barry Lynn via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@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 Jagan:<div><br></div><div>Do you understand the error message, but are not sure why you are getting it?</div><div><br></div><div>The message means that the parameters within your interpolation line are not dimensionally consistent. If you look at the page that describes the function, you might be able to figure out which of your variables does not conform dimensionally. </div><div><br></div><div>If you've already gone past this point of trying to understand the error, then please write back.</div><div><br></div><div>Barry</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Aug 17, 2019 at 8:06 PM jagan TNAU via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@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 dir="ltr"><div dir="ltr"><div dir="ltr">Dear Users,<div><br></div><div>I was trying to interpolate gridded data to some station locations using linint2_points function with the following script.</div><div><div>;----------------------------------------------------------------------</div><div>; lin2points.ncl</div><div>;----------------------------------------------------------------------</div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div><br></div><div>begin</div><div>;************************************************</div><div>; read ascii input file</div><div>;************************************************</div><div> values = asciiread("/Users/Jag/Desktop/cdz-regrid/var.txt",-1,"float")</div><div> nvals = dimsizes(values)</div><div> ncols = 8</div><div> nrows = ncols/nvals</div><div> lat = values(0::ncols)</div><div> lon = values(1::ncols) </div><div> srad = values(2::ncols) </div><div> tmax = values(3::ncols)</div><div> tmin = values(4::ncols) </div><div> rain = values(5::ncols)</div><div> wind = values(6::ncols) </div><div> rhum = values(7::ncols) </div><div>;************************************************</div><div><br></div><div>; Lat*lon array 11*11 = 121 rows read from the files</div><div>; lat = (/10.0, 10.2, 10.4, 10.6, 10.8, 11.0, 11.2, 11.4, 11.6, 11.8, 12.0/)</div><div>; lon = (/78.0, 78.2, 78.4, 78.6, 78.8, 79.0, 79.2, 79.4, 79.6, 79.8, 80.0/)</div><div><br></div><div>; points for which interpolation required totally 30</div><div> nlat=(/10.38, 10.40/)</div><div> nlon=(/79.85, 79.49/)</div><div><br></div><div> ntmax = linint2_points(lat,lon,tmax,False,nlat,nlon,0)</div><div><br></div><div>; Print only tmax as a sample</div><div> print(ntmax)</div><div><br></div><div>; nsrad = linint2_points(lat,lon,srad,False,nlat,nlon,0)</div><div>; ntmin = linint2_points(lat,lon,tmin,False,nlat,nlon,0)</div><div>; nrain = linint2_points(lat,lon,rain,False,nlat,nlon,0)</div><div>; nwind = linint2_points(lat,lon,wind,False,nlat,nlon,0)</div><div>; nrhum = linint2_points(lat,lon,rhum,False,nlat,nlon,0)</div><div><br></div><div>; alist = [/nsrad,ntmax,ntmin,nrain,nwind,nrhum/]</div><div>; write_table("tmp11.txt", "w", alist, "%6.1f%6.1f%6.1f%7.1f%6.1f%6.1f")</div><div><br></div><div>end</div><div>-----------------------------------------------</div><div>But I encountered the following error </div><div><p style="margin:0px;font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><font size="4">fatal:linint2_points: The rightmost </font>dimensions of fi must be nyi x nxi, where nyi and nxi are the lengths of yi and xi respectively</span></p>
<p style="margin:0px;font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">fatal:["Execute.c":8637]:Execute: Error occurred at or near line 32 in file lin2points.ncl</span></p><p style="margin:0px;font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures"><br></span></p><p style="margin:0px;font-stretch:normal;line-height:normal;font-family:Menlo;color:rgb(0,0,0)"><span style="font-variant-ligatures:no-common-ligatures">Will you please help me to solve the problem. Sample file attached</span></p></div><div><br></div><div><br></div>-- <br><div dir="ltr" class="gmail-m_853049223304161494gmail-m_-6018463958552543820gmail-m_5886783158412746250gmail_signature"><div dir="ltr"><div>With regards<br><br>Dr.R.Jagannathan <br>Professor & Former Dean<br>Tamil Nadu Agricultural University <br>Coimbatore - 641 003 India <br><br>PHONE: Mob: +91 94438 89891<br><br>DO NOT PRINT THIS E-MAIL UNLESS NECESSARY. THE ENVIRONMENT CONCERNS US ALL.</div></div></div></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></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail-m_853049223304161494gmail-m_-6018463958552543820gmail_signature"><div dir="ltr"><div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Associate Scientist, Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><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></blockquote></div>