<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" 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="x_result_box" class="" lang="en"><span class="">Hi NCLs</span></span></p>
<p style="margin-top:0; margin-bottom:0"><span id="x_result_box" class="" lang="en"><span class=""><br>
</span></span></p>
<p style="margin-top:0; margin-bottom:0"><span id="x_result_box" class="" lang="en"><span class="">I would like to create an array of 36 x 12 (time x months)</span><span class=""> for each lat and lon values from file (91 values), generating 91
</span><span id="result_box" class="short_text" lang="en"><span class="">separate files</span></span>.<br>
<span class="">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="x_result_box" class="" lang="en"><span class=""><br>
</span></span></p>
<p style="margin-top:0; margin-bottom:0"><span id="x_result_box" class="" lang="en"><span class=""></span></span></p>
<div> lats = asciiread("lat.txt",(/91,3/),"float")<br>
lons = asciiread("lon..txt",(/91,3/),"float")<br>
</div>
<div> lat = lats(:,1) <br>
lon = lons(:,2) </div>
<p style="margin-top:0; margin-bottom:0"></p>
<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="x_result_box" class="x_short_text" lang="en"><span class="">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" <br>
write_matrix( var_matriz, dimsizes(list)+"f10.2", arq ) </div>
</div>
<br>
<p></p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>De:</b> Vanúcia Schumacher <vanucia-schumacher@hotmail.com><br>
<b>Enviado:</b> terça-feira, 16 de janeiro de 2018 11:11:09<br>
<b>Para:</b> ncl-talk@ucar.edu<br>
<b>Assunto:</b> help loop</font>
<div> </div>
</div>
<style type="text/css" style="display:none">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_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="x_result_box" class="" lang="en"><span class="">Hi NCLs</span></span></p>
<p style="margin-top:0; margin-bottom:0"><span id="x_result_box" class="" lang="en"><span class=""><br>
</span></span></p>
<p style="margin-top:0; margin-bottom:0"><span id="x_result_box" class="" lang="en"><span class="">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 class="">Also calculate a file (36 x 12) for each lat and lon values from another file (91 values).</span><br>
<span class="">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="x_result_box" class="" lang="en"><span class=""><br>
</span></span></p>
<p style="margin-top:0; margin-bottom:0"><span id="x_result_box" class="" lang="en"><span class=""></span></span></p>
<div> lats = asciiread("lat.txt",(/91,3/),"float")<br>
lons = asciiread("lon..txt",(/91,3/),"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="x_result_box" class="x_short_text" lang="en"><span class="">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></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>
</span><br>
<p></p>
</div>
</div>
</body>
</html>