Ranyere,<br><br>&nbsp; As an alternative, here is how you would do his using PyGrADS:<br><br>---------------------------------------------------------------------------------------------------<br>import grads<br>ga = grads.GaNum(Bin=&#39;gradsnc&#39;,Echo=False,Window=False)&nbsp; # start grads<br>
fh = ga.open(&quot;../data/model.nc&quot;)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # open file<br>ts = ga.exp(&#39;ts&#39;)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; # export a variable to python<br>print &quot;&nbsp;&nbsp; Lon&nbsp;&nbsp;&nbsp;&nbsp; Lat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ts&quot;<br>print &quot;-------- -------- --------&quot;<br>
for j in range(fh.ny):<br>&nbsp;&nbsp;&nbsp; for i in range(fh.nx):<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; print &quot;%8.3f %8.3f %8.2f&quot;%(ts.grid.lon[i],ts.grid.lat[j],ts[j,i])<br>----------------------------------------------------------------------------------------------------------<br>
which produces:<br><br>&nbsp; Lon&nbsp;&nbsp;&nbsp;&nbsp; Lat&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ts<br>-------- -------- --------<br>&nbsp;&nbsp; 0.000&nbsp; -90.000&nbsp;&nbsp; 258.49<br>&nbsp;&nbsp; 5.000&nbsp; -90.000&nbsp;&nbsp; 258.49<br>&nbsp; 10.000&nbsp; -90.000&nbsp;&nbsp; 258.49<br>&nbsp; 15.000&nbsp; -90.000&nbsp;&nbsp; 258.49<br>&nbsp; 20.000&nbsp; -90.000&nbsp;&nbsp; 258.49<br>
&nbsp; 25.000&nbsp; -90.000&nbsp;&nbsp; 258.49<br>&nbsp; ...<br>&nbsp;230.000&nbsp; -18.000&nbsp;&nbsp; 300.49<br>&nbsp;235.000&nbsp; -18.000&nbsp;&nbsp; 299.99<br>&nbsp;240.000&nbsp; -18.000&nbsp;&nbsp; 299.49<br>&nbsp;245.000&nbsp; -18.000&nbsp;&nbsp; 298.99<br>...<br><br>An equivalent script using the Perl interface would be very similar, since the methods open() and exp() are also available there.<br>
<br>&nbsp;&nbsp; Arlindo<br><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>