<div dir="ltr">Hello:<div><br></div><div>Have you considered to use get_ij.f</div><div><br></div><div>It will give 1 or end of domain values for values outside the domain.</div><div><br></div><div>





<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>i_loc := new(num_obs,integer)</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space"> </span>j_loc := new(num_obs,integer)</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space">   </span>ii_loc = 0</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space">   </span>jj_loc = 0</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space">   </span>do i_sta = 0,num_obs-1</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space">     </span>GET_IJ::</span><span class="gmail-s2" style="font-variant-ligatures:no-common-ligatures;background-color:rgb(224,228,9)">get_ij</span><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures">(xlat2d,xlon2d,lat_o(i_sta),lon_o(i_sta),ii_loc,jj_loc,dims2d(0),dims2d(1))</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space">     </span>i_loc(i_sta) = ii_loc-1</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space">     </span>j_loc(i_sta) = jj_loc-1</span></p>
<p class="gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:22px;line-height:normal;font-family:"Courier New";color:rgb(0,0,0)"><span class="gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span class="gmail-Apple-converted-space">   </span>end do</span></p></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Sep 11, 2019 at 7:32 PM Ehsan Taghizadeh via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">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><div class="gmail-m_-6490407475172893496ydpd4ce19beyahoo-style-wrap" style="font-family:"bookman old style","new york",times,serif"><div style="font-size:16px"></div>
        <div dir="ltr" style="font-size:16px">Yes, exactly.<br></div><div dir="ltr" style="font-size:16px">I want to interpolate WRF data on station locations, which have been included in "stations.dat" file.</div><div dir="ltr" style="font-size:16px">The problem is that some of the stations in this file don't locate in WRF domain, but most of them get values. If the station doesn't locate in the domain, how it could get interpolated values. For example, 40899 and 40890 are out of domain 2, and they got bad value "<span>9969209968386869046778552952102584320.0</span>", correctly. On the other hand, 40745, also is out of domain 2, but it got value! Logically, It must have "<span>9969209968386869046778552952102584320.0</span>", too.</div><div style="font-size:16px"><br></div><div dir="ltr"><i><font size="2">Sincerely</font></i></div><div dir="ltr"><i><font size="2">Ehsan</font></i></div>
        
        </div><div id="gmail-m_-6490407475172893496yahoo_quoted_9146187504" class="gmail-m_-6490407475172893496yahoo_quoted">
            <div style="font-family:"Helvetica Neue",Helvetica,Arial,sans-serif;font-size:13px;color:rgb(38,40,42)">
                
                <div>
                    On Wednesday, September 11, 2019, 07:16:56 PM GMT+4:30, Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:
                </div>
                <div><br></div>
                <div><br></div>
                <div><div id="gmail-m_-6490407475172893496yiv7289053809"><div><div dir="ltr">Hello:<div><br clear="none"></div><div>DId you mean that you want to interpolate WRF data to station locations?</div></div><br clear="none"><div class="gmail-m_-6490407475172893496yiv7289053809gmail_quote"><div class="gmail-m_-6490407475172893496yiv7289053809yqt3435483563" id="gmail-m_-6490407475172893496yiv7289053809yqt21933"><div class="gmail-m_-6490407475172893496yiv7289053809gmail_attr" dir="ltr">On Wed, Sep 11, 2019 at 3:08 PM Ehsan Taghizadeh via ncl-talk <<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br clear="none"></div><blockquote class="gmail-m_-6490407475172893496yiv7289053809gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="gmail-m_-6490407475172893496yiv7289053809gmail-m_-1845127014997727369ydp855ff780yahoo-style-wrap"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr" style="font-size:16px"><font face="times new roman, new york, times, serif">(Pardon me for incomplete previous email)</font></div><div style="font-size:16px"><font face="times new roman, new york, times, serif">Dears,</font></div><div style="font-size:16px"><font face="times new roman, new york, times, serif">I hope this email finds you well.</font></div><div style="font-size:16px"><font face="times new roman, new york, times, serif">I have a script for interpolating WRF data on some stations, using "rcm2points". The input file for stations includes stations which some of them don't locate on wrf out domain 2. I thought these stations must get missing value or so. But all stations, except a few ones, got values!</font></div><div style="font-size:16px"><font face="times new roman, new york, times, serif">May I ask you how it could possible?</font></div><div style="font-size:16px"><font face="times new roman, new york, times, serif">I've attached the script (wrfncl.ncl), the input file for stations (stations.dat), an image showing domain 2 (wrfd2.png), also an image showing all stations locations (stations.png). The wrfout file is about 2G. Should I put this file in ftp?</font></div><div dir="ltr" style="font-size:16px"><font face="times new roman, new york, times, serif">Some information of my system is as below:</font></div><div dir="ltr"><font face="courier new, courier, monaco, monospace, sans-serif" size="2">$ uname -a<br clear="none"></font></div><div dir="ltr"><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Linux localhost.localdomain 3.10.0-514.2.2.el7.x86_64 #1 SMP Tue Dec 6 23:06:41 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux</font><br clear="none"></div><div dir="ltr"><font face="courier new, courier, monaco, monospace, sans-serif" size="2"><span>$ ncl -V</span><br clear="none"></font></div><div dir="ltr"><font face="courier new, courier, monaco, monospace, sans-serif" size="2"><span><span>6.6.2</span></span></font></div><div style="font-size:16px"><font face="times new roman, new york, times, serif"><br clear="none"></font></div><div dir="ltr" style="font-size:16px"><font face="times new roman, new york, times, serif">Also, about some variables:</font></div><div dir="ltr"><div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Variable: wlat</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Type: float</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Total Size: 383616 bytes</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">            95904 values</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Number of Dimensions: 2</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Dimensions and sizes:   [south_north | 333] x [west_east | 288]</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Coordinates:</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Number Of Attributes: 6</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">  FieldType :   104</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">  MemoryOrder : XY</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">  description : LATITUDE, SOUTH IS NEGATIVE</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">  units :       degree_north</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">  stagger :</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">  coordinates : XLONG XLAT</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">(0)     LATITUDE, SOUTH IS NEGATIVE (degree_north) : min=31.2974   max=40.4085</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">(0)     LONGITUDE, WEST IS NEGATIVE (degree_east) : min=43.4514   max=53.5448</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2"><br clear="none"></font></div><div dir="ltr"><div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Variable: TC2</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Type: float</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Total Size: 8055936 bytes</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">            2013984 values</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Number of Dimensions: 3</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Dimensions and sizes:   [21] x [333] x [288]</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Coordinates:</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">Number Of Attributes: 1</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2">  units :       C</font></div><div><font face="courier new, courier, monaco, monospace, sans-serif" size="2"><br clear="none"></font></div><div dir="ltr"><font face="courier new, courier, monaco, monospace, sans-serif" size="2"></font><div><div>Variable: TC2_st</div><div>Type: float</div><div>Total Size: 32172 bytes</div><div>            8043 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:   [21] x [383]</div><div>Coordinates:</div><div>Number Of Attributes: 1</div><div>  _FillValue :  9.96921e+36</div><div><br clear="none"></div></div><div><div>Variable: SM</div><div>Type: float</div><div>Total Size: 8055936 bytes</div><div>            2013984 values</div><div>Number of Dimensions: 3</div><div>Dimensions and sizes:   [Time | 21] x [south_north | 333] x [west_east | 288]</div><div>Coordinates:</div><div>Number Of Attributes: 6</div><div>  FieldType :   104</div><div>  MemoryOrder : XYZ</div><div>  description : SOIL MOISTURE</div><div>  units :       m3 m-3</div><div>  stagger :     Z</div><div>  coordinates : XLONG XLAT XTIME</div><div><br clear="none"></div><div>Variable: SM_st</div><div>Type: float</div><div>Total Size: 32172 bytes</div><div>            8043 values</div><div>Number of Dimensions: 2</div><div>Dimensions and sizes:   [21] x [383]</div><div>Coordinates:</div><div>Number Of Attributes: 1</div><div>  _FillValue :  9.96921e+36</div><div><br clear="none"></div></div></div><div style="font-size:16px">I'll be thankful if I hear from you.<br clear="none"></div></div></div></div></div><div style="font-size:16px"><br clear="none"></div></div></div><div dir="ltr" style="font-size:16px"><i><font size="2">Sincerely</font></i></div><div dir="ltr" style="font-size:16px"><i><font size="2">Ehsan</font></i></div></div><br clear="none"></div></div></div>_______________________________________________<br clear="none">
ncl-talk mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br clear="none">
List instructions, subscriber options, unsubscribe:<br clear="none">
<a rel="nofollow" shape="rect" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div></div><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div class="gmail-m_-6490407475172893496yiv7289053809gmail_signature" dir="ltr"><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 clear="none" style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br clear="none" style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br clear="none" style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br clear="none" 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 clear="none"></span></div><div>C.E.O, Weather It Is, LTD<br clear="none">Weather and Climate Focus<br clear="none"><a rel="nofollow" shape="rect" href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br clear="none">Jerusalem, Israel<br clear="none">Local: 02 930 9525<br clear="none">Cell: 054 7 231 170<br clear="none">Int-IS: x972 2 930 9525<br clear="none"><br clear="none"></div></div></div></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_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>