<div dir="ltr"><div class="gmail_default" style="font-size:small">Vanúcia,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">You didn't say exactly what the problem was, but write_matrix doesn't append data. So, if you call it repeatedly with the same file, it's just going to overwrite that file with the most recent data you give it.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">In the "k" do loop, you are setting the output file name based on the counter from the "i" do loop. This means that the file name is the same for every iteration of k, and hence will get overwritten each time.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I suggest looking at write_table instead, which allows you to append data.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/write_table.shtml">http://www.ncl.ucar.edu/Document/Functions/Built-in/write_table.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">There are some examples at:</div><div class="gmail_default"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/write_ascii.shtml">http://www.ncl.ucar.edu/Applications/write_ascii.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">In particular, look at write_ascii_4.ncl.</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jan 16, 2018 at 6:52 AM, Vanúcia Schumacher via ncl-talk <span dir="ltr"><<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div dir="ltr">
<div id="m_-1804434117919121639divtagdefaultwrapper" style="font-size:14pt;color:#000000;font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p style="margin-top:0;margin-bottom:0"></p>
<div>
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span>Hi NCLs</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span><br>
</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span>I would like to create an array of 36 x 12 (time x months)</span><span> for each lat and lon values from file (91 values), generating 91
</span><span id="m_-1804434117919121639result_box" class="m_-1804434117919121639short_text" lang="en"><span>separate files</span></span>.<span class=""><br>
<span>I tried the script below, but I am not able to generate a loop for what I want to get, if anyone can give a help I appreciate.</span></span></span></p><span class="">
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span><br>
</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span></span></span></p>
<div> lats  = asciiread("lat.txt",(/91,3/),"<wbr>float")<br>
 lons = asciiread("lon..txt",(/91,3/),<wbr>"float")<br>
</div>
<div> lat    = lats(:,1)   <br>
 lon   = lons(:,2) </div>
<p style="margin-top:0;margin-bottom:0"></p>
</span><div>nlat   = dimsizes(lat)   <br>
nlon  = dimsizes(lon)</div>
<p></p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0"><span> list = systemfunc("ls /cru.*.nc")  ;
<span id="m_-1804434117919121639x_result_box" class="m_-1804434117919121639x_short_text" lang="en"><span>are 12 files separated by months</span></span></span></p>
<p style="margin-top:0;margin-bottom:0"><span><span> </span></span></p>
<div>var_matriz = new( (/36,dimsizes(list)/), float )<br>
      do j = 0,dimsizes(list)-1<br>
<br>
do i=0,nlat(0)-1<br>
do k=0,nlon(0)-1<br>
 <br>
       a = addfile(list(j), "r")    <br>
       var_matriz(:,j) = a->pr(:,i,k) <br>
<br>
    <span>arq=True</span>        <br>
   arq@fout="ex."+i+".txt"       <wbr>   <br>
   write_matrix( var_matriz, dimsizes(list)+"f10.2", arq ) </div>
</div>
<br>
<p></p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_-1804434117919121639divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>De:</b> Vanúcia Schumacher <<a href="mailto:vanucia-schumacher@hotmail.com" target="_blank">vanucia-schumacher@hotmail.<wbr>com</a>><br>
<b>Enviado:</b> terça-feira, 16 de janeiro de 2018 11:11:09<br>
<b>Para:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<b>Assunto:</b> help loop</font>
<div> </div>
</div><div><div class="h5">

<div dir="ltr">
<div id="m_-1804434117919121639x_divtagdefaultwrapper" dir="ltr" style="font-size:14pt;color:#000000;font-family:Calibri,Helvetica,sans-serif">
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span>Hi NCLs</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span><br>
</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span>I would like to create an array of 36 x 12 (time x months), but I can only generate one file at a time.</span>
<span>Also calculate a file (36 x 12) for each lat and lon values from another file (91 values).</span><br>
<span>I tried the script below, but I am not able to generate a loop for what I want to get, if anyone can give a help I appreciate.</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span><br>
</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span id="m_-1804434117919121639x_result_box" lang="en"><span></span></span></p>
<div> lats  = asciiread("lat.txt",(/91,3/),"<wbr>float")<br>
 lons = asciiread("lon..txt",(/91,3/),<wbr>"float")<br>
</div>
<div> lat    = lats(:,1)   <br>
 lon    = lons(:,2) </div>
<p></p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0"><span> f = systemfunc("ls /cru.*.nc")  ;
<span id="m_-1804434117919121639x_result_box" class="m_-1804434117919121639x_short_text" lang="en"><span>are 12 files separated by months</span></span></span></p>
<p style="margin-top:0;margin-bottom:0"><span><span> N = dimsizes( f )</span></span></p>
<p style="margin-top:0;margin-bottom:0"><span></span></p>
<div> do i = 0,N-1<br>
<div>a = addfiles(f(i), "r")    <br>
var = a[:]->pr(:,{lat},{lon}) </div>
<div> asciiwrite ("prec.cru."+i+".txt", sprintf("%9.2f", var))    <br>
 end do<br>
<br>
</div>
<br>
</div>
<br>
<p></p>
</div>
</div>
</div></div></div>

<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>