Hi, I think that you can just do:<br><br><div><p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">ncl 0&gt; A = (/ 1, 2, 3, 4, 5, 6, 7, 8, 9/)</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">ncl 1&gt; B = (/ 2, 14, 7, 1/)</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">ncl 2&gt; &nbsp;</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">ncl 3&gt; C = A(ind(.not.ismissing(get1Dindex(B,A))))</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">ncl 4&gt; print(C)</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); min-height: 16px;"><br></p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0); min-height: 16px;"><br></p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">Variable: C</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">Type: integer</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">Total Size: 12 bytes</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 3 values</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">Number of Dimensions: 1</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">Dimensions and sizes: &nbsp; [3]</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">Coordinates:&nbsp;</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">(0) &nbsp; &nbsp; 1</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">(1) &nbsp; &nbsp; 2</p>
<p style="margin: 0px; font-size: 12px; font-family: Monaco; color: rgb(245, 245, 245); background-color: rgb(0, 0, 0);">(2) &nbsp; &nbsp; 7</p><br>Good luck!</div><div><br></div><div>Yongjie Huang<br><br>&gt;&nbsp;Date:&nbsp;Fri,&nbsp;19&nbsp;Dec&nbsp;2014&nbsp;08:39:55&nbsp;-0700<br>&gt;&nbsp;From:&nbsp;Dennis&nbsp;Shea&nbsp;&lt;shea@ucar.edu&gt;<br>&gt;&nbsp;Subject:&nbsp;Re:&nbsp;[ncl-talk]&nbsp;How&nbsp;to&nbsp;find&nbsp;the&nbsp;same&nbsp;elements&nbsp;in&nbsp;two&nbsp;matrixes<br>&gt;&nbsp;To:&nbsp;Lin&nbsp;Wang&nbsp;&lt;wanglin@post.iap.ac.cn&gt;<br>&gt;&nbsp;Cc:&nbsp;ncl-talk&nbsp;&lt;ncl-talk@ucar.edu&gt;<br>&gt;&nbsp;Message-ID:<br>&gt;&nbsp;        &lt;CAOF1d_5DLMe=Tv794H9ekb598XgNKcuf9AtES=9X8ABhj9pFWg@mail.gmail.com&gt;<br>&gt;&nbsp;Content-Type:&nbsp;text/plain;&nbsp;charset="utf-8"<br>&gt;&nbsp;<br>&gt;&nbsp;I&nbsp;do&nbsp;not&nbsp;think&nbsp;you&nbsp;can&nbsp;do&nbsp;it&nbsp;without&nbsp;a&nbsp;loop.<br>&gt;&nbsp;<br>&gt;&nbsp;function&nbsp;get_dup_elements(a[*],&nb
 sp;b[*])<br>&gt;&nbsp;begin<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;dima&nbsp;=&nbsp;dimsizes(a)<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;dimb&nbsp;=&nbsp;dimsizes(b)<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;ndup&nbsp;=&nbsp;-1<br>&gt;&nbsp;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;if&nbsp;(dima.le.dimb)&nbsp;then<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;dimd&nbsp;=&nbsp;new(&nbsp;dima,&nbsp;typeof(a))<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;do&nbsp;n=0,dima-1<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;ii&nbsp;:=&nbsp;ind(b.eq.a(n))<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;if&nbsp;(.not.ismissing(ii(0)))&nbsp;then<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;ndup&nbsp;=&nbsp;&nbsp;ndup+1<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;dimd(ndup)&nbsp;=&nbsp;a(n)<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n
 bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;end&nbsp;if<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;end&nbsp;do<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;else<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;dimd&nbsp;=&nbsp;new(&nbsp;dimb,&nbsp;typeof(a))<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;do&nbsp;n=0,dimb-1<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;ii&nbsp;:=&nbsp;ind(a.eq.b(n))<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;if&nbsp;(.not.ismissing(ii(0)))&nbsp;then<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;ndup&nbsp;=&nbsp;ndup+1<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;dimd(ndup)&nbsp;=&nbsp;b(n)<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;end&nbsp;if<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;end&nbsp;do<br>&gt;&nbsp;<br>&gt;&nbsp;
 &nbsp;&nbsp;&nbsp;end&nbsp;if<br>&gt;&nbsp;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;return(dimd(0:ndup))<br>&gt;&nbsp;end<br>&gt;&nbsp;<br>&gt;&nbsp;A&nbsp;=&nbsp;(/&nbsp;1,&nbsp;2,&nbsp;3,&nbsp;4,&nbsp;5,&nbsp;6,&nbsp;7,&nbsp;8,&nbsp;9/)<br>&gt;&nbsp;B&nbsp;=&nbsp;(/&nbsp;2,&nbsp;14,&nbsp;7,&nbsp;1/)<br>&gt;&nbsp;<br>&gt;&nbsp;dup&nbsp;=&nbsp;get_dup_elements(A,B)<br>&gt;&nbsp;print(dup)<br>&gt;&nbsp;<br>&gt;&nbsp;;;&nbsp;sort(dup)&nbsp;&nbsp;;&nbsp;ascending&nbsp;order<br>&gt;&nbsp;<br>&gt;&nbsp;<br>&gt;&nbsp;<br>&gt;&nbsp;<br>&gt;&nbsp;On&nbsp;Wed,&nbsp;Dec&nbsp;17,&nbsp;2014&nbsp;at&nbsp;1:59&nbsp;AM,&nbsp;Lin&nbsp;Wang&nbsp;&lt;wanglin@post.iap.ac.cn&gt;&nbsp;wrote:<br>&gt;&nbsp;<br>&gt;&nbsp;&gt;&nbsp;Dear&nbsp;all,<br>&gt;&nbsp;&gt;<br>&gt;&nbsp;&gt;&nbsp;Suppose:<br>&gt;&nbsp;&gt;&nbsp;A&nbsp;=&nbsp;(/&nbsp;1,&nbsp;2,&nbsp;3,&nbsp;4,&nbsp;5,&nbsp;6,&nbsp;7,&nbsp;8,&nbsp;9/)<br>&gt;&nbsp;&gt;&nbsp;B&nbsp;=&nbsp;(/&nbsp;2,&nbsp;14,&nbsp;7,&nbsp;1/)<br>&gt;&nbsp;&gt;<br>&gt;&nbs
 p;&gt;&nbsp;If&nbsp;don't&nbsp;use&nbsp;loop,&nbsp;how&nbsp;to&nbsp;find&nbsp;the&nbsp;same&nbsp;elements&nbsp;between&nbsp;A&nbsp;and&nbsp;B&nbsp;such&nbsp;as<br>&gt;&nbsp;&gt;&nbsp;(/1,&nbsp;2,&nbsp;7/)?&nbsp;Thanks.<br>&gt;&nbsp;&gt;<br>&gt;&nbsp;&gt;&nbsp;Best&nbsp;wishes,<br>&gt;&nbsp;&gt;&nbsp;Lin&nbsp;Wang<br>&gt;&nbsp;&gt;<br><span><br><div><font face="Comic Sans MS">-----------------------------------------------------------------------------------</font></div><font face="Comic Sans MS">HUANG&nbsp;Yong-jie<br>Key&nbsp;Laboratory&nbsp;of&nbsp;Cloud-Precipitation&nbsp;Physics&nbsp;and&nbsp;Severe&nbsp;Storms&nbsp;(LACS),&nbsp;<br>Institute&nbsp;of&nbsp;Atmospheric&nbsp;Physics,&nbsp;Chinese&nbsp;Academy&nbsp;of&nbsp;Sciences<br>Email&nbsp;Address:&nbsp;huangyj@mail.iap.ac.cn&nbsp;OR&nbsp;huangynj@gmail.com&nbsp;</font><div><span style="font-family: 'Comic Sans MS';">-----------------------------------------------------------------------------------</span><font face="C
 omic Sans MS"><br></font><br></div></span></div><br><br><br>