<div dir="ltr">It turns out the data in the wrfinput file was transcribed correctly.<div><br></div><div>So, the programs I sent before to read a variable and substitute data in a wrfinput file variable work fine.<br><div><br></div><div>The problem was in accessing the array within WRF.<br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 7, 2017 at 5:29 PM, Barry Lynn <span dir="ltr">&lt;<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Thank you for giving it a go.<br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Barry Lynn</b> <span dir="ltr">&lt;<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>&gt;</span><br>Date: Mon, Mar 6, 2017 at 6:55 PM<br>Subject: An NCL Mystery (regarding writing to wrfinput_d03)<br>To: ncl-talk &lt;<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>&gt;<br><br><br><div dir="ltr">Hello:<div><br></div><div>I wrote a program (feed_data.ncl -- with a bit of help from Dennis, Mary, etc) to read an ascii file and write the data to the variable &quot;l_obs&quot; in the <a href="http://wrfinput_d03.nc" target="_blank">wrfinput_d03.nc</a> file.</div><div><br></div><div>When I read the asci file and then write the data to <a href="http://wrfinput_d03.nc" target="_blank">wrfinput_d03.nc</a>, and then plot the data (using l_obs_map.ncl) it looks absolutely fine (attached). </div><div><br></div><div>However, when I access the data within the WRF itself it is not located at the correct indices.  For instance, points that should be at about latitude 30.95, longitude 36.6 are at latitude 28, longitude 35).</div><div><br></div><div>I realize that there is a switching of indices from the NCL to the WRF, but I assumed that this was done correctly/automatically​<div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:rgb(245,245,245);padding:5px;color:rgb(34,34,34);font-family:arial;font-style:normal;font-weight:bold;font-size:13px;border:1px solid rgb(221,221,221);line-height:1"><a href="https://drive.google.com/file/d/0B3E3GjUeKpBdMTJGOHd0R00ya0E/view?usp=drive_web" style="display:inline-block;max-width:366px;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0;border:none" target="_blank"><img style="vertical-align:bottom;border:none" src="https://ssl.gstatic.com/docs/doclist/images/icon_10_generic_list.png"> <span dir="ltr" style="color:rgb(17,85,204);text-decoration:none;vertical-align:bottom">lobs_output_1.3km.txt</span></a><img style="opacity:0.55;float:right;display:none"></div>​ as the size of the array in wrfinput_d03 is correct. </div><div><br></div><div>Just to reiterate, I can read the data with ncl from <a href="http://wrfinput_d03.nc" target="_blank">wrfinput_d03.nc</a> and map it (and it corresponds to what I see using a different mapping program).  However, when I access the array in the WRF, as l_obs(i,k_lev,j), where k_lev is the time increment, it is not in the correct location.</div><div><br></div><div>Thank you for your suggestions.</div><div><br></div><div>Barry </div><div><br></div><div><br></div><div><br clear="all"><div><br></div>-- <br><div class="m_-1182482505884624583m_-4288084641722880064gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior 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>US <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div></div>
</div><br><br clear="all"><div><br></div>-- <br><div class="m_-1182482505884624583gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior 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>US <a href="tel:(914)%20432-3108" value="+19144323108" target="_blank">914 432 3108</a><br></div></div></div>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior 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>US 914 432 3108<br></div></div></div>
</div></div></div></div>