Thanks Charles i go to try this... thanks... Regards...<br><br><div><span class="gmail_quote">2007/5/22, Charles Seman &lt;<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>&gt;:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Rafael,<br><br>Please find attached a short GrADS script called &quot;test_xy2w.gs&quot; (and the<br>resulting output file &quot;test_xy2w.txt&quot;) with code to plot some &quot;i&quot;, &quot;j&quot;,<br>&quot;lonval&quot;, and &quot;latval&quot; values into an output file &quot;test_xy2w.txt&quot; with a
<br>specified format...<br><br>see also:<br><a href="http://grads.iges.org/grads/gadoc/mathfunctions.html">http://grads.iges.org/grads/gadoc/mathfunctions.html</a><br>for info on formatting numbers for output (I couldn&#39;t get the &quot;prnopts&quot;
<br>command to do this)...<br><br>I hope this helps,<br>Chuck<br><br>Rafael Mundaray wrote:<br>&gt; hi GrADS Users, Anybody here know how i write a file with the<br>&gt; coordinates of my grid...<br>&gt; example:<br>&gt; x y lon&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;lat
<br>&gt; 1 1 -67.3&nbsp;&nbsp; 3.1<br>&gt; .&nbsp;&nbsp;.&nbsp;&nbsp;&nbsp;&nbsp; .&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;.<br>&gt; .&nbsp;&nbsp;.&nbsp;&nbsp;&nbsp;&nbsp; .&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;.<br>&gt; .&nbsp;&nbsp;.&nbsp;&nbsp;&nbsp;&nbsp; .&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;.<br>&gt; .&nbsp;&nbsp;.&nbsp;&nbsp;&nbsp;&nbsp; .&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;.<br>&gt; n n&nbsp;&nbsp;&nbsp;&nbsp; n&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n<br>&gt;<br>&gt; Can you help me with this???<br>&gt; I know that &#39;q xy2w x y&#39; corvets the x and y dimension to reals
<br>&gt; coordinates (Lon,Lat), but i want that save this in one file with all<br>&gt; my points of my grids... I was trying this ( bold type) but didn&#39;t<br>&gt; function...<br>&gt;<br>&gt; &#39;set gxout print&#39;<br>
&gt; &#39;set prnopts %1.3f 1 1 u&#39;<br>&gt;<br>&gt; &#39;q dims&#39;<br>&gt;<br>&gt; rec=sublin(result,2)<br>&gt; Imin=subwrd(rec,11)<br>&gt; Imax=subwrd(rec,13)<br>&gt; say &#39;Imin &#39;Imin&nbsp;&nbsp; &#39;Imax &#39; Imax<br>
&gt; rec=sublin(result,3)<br>&gt; Jmin=subwrd(rec,11)<br>&gt; Jmax=subwrd(rec,13)<br>&gt; say &#39;Jmin &#39;Jmin&nbsp;&nbsp; &#39;Jmax &#39; Jmax<br>&gt;<br>&gt; I = 1<br>&gt; while (I =&lt; Imax)<br>&gt; &#39;q xy2w &#39;I &#39; &#39; J
<br>&gt; rec=sublin(result,1)<br>&gt; I=subwrd(rec,3)<br>&gt; J=subwrd(rec,6)<br>&gt; endwhile<br>&gt;&nbsp;&nbsp; while (J =&lt; Jmax)<br>&gt;&nbsp;&nbsp; &#39;q xy2w &#39;I &#39; &#39; J<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;rec=sublin(result,1)<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;I=subwrd(rec,3)
<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;J=subwrd(rec,6)<br>&gt; rc = write (&#39;/root/Desktop/LatLon.txt&#39;,result)<br>&gt; rc = close (&#39;/root/Desktop/LatLon.txt&#39;)<br>&gt; endwhile<br>&gt; --<br>&gt; Ing. Rafael E. Mundaray M.<br><br>--<br>
<br>Please note that <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a> should be considered my NOAA<br>email address, not <a href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a>.<br><br>********************************************************************
<br> Charles Seman&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;&nbsp;&nbsp;&nbsp;&nbsp;<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br> U.S. Department of Commerce / NOAA / OAR<br> Geophysical Fluid Dynamics Laboratory&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; voice: (609) 452-6547
<br> 201 Forrestal Road&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;&nbsp;&nbsp;fax: (609) 987-5063<br> Princeton, NJ&nbsp;&nbsp;08540-6649&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<a href="http://www.gfdl.noaa.gov/~cjs/">http://www.gfdl.noaa.gov/~cjs/</a><br>********************************************************************
<br><br>&quot;The contents of this message are mine personally and do not reflect<br>any position of the Government or NOAA.&quot;<br><br><br>&#39;reinit&#39;<br>&#39;sdfopen /net/stg/netcdf/isabel_ras.nc&#39;<br>*<br>*&nbsp;&nbsp;define scaling environment by displaying any variable (otherwise an error
<br>*&nbsp;&nbsp;results...)<br>*<br>&#39;set gxout shaded&#39;<br>&#39;d ps&#39;<br>*<br>*&nbsp;&nbsp;open output file and print header...<br>*<br>rc = write(&#39;test_xy2w.txt&#39;,&#39;i j&nbsp;&nbsp;lonval latval&#39;)<br>i=1<br>while ( i &lt;= 2 )
<br>&nbsp;&nbsp;j=1<br>&nbsp;&nbsp;while ( j &lt;= 2 )<br>&nbsp;&nbsp;&nbsp;&nbsp;&#39;q xy2w &#39;i&#39; &#39;j<br>&nbsp;&nbsp;&nbsp;&nbsp;rec = sublin(result,1)<br>&nbsp;&nbsp;&nbsp;&nbsp;lonv = subwrd(rec,3)<br>&nbsp;&nbsp;&nbsp;&nbsp;latv = subwrd(rec,6)<br>&nbsp;&nbsp;&nbsp;&nbsp;say &#39;&#39;i&#39; &#39;j&#39; &#39;lonv&#39; &#39;latv
<br>&nbsp;&nbsp;&nbsp;&nbsp;lonval = math_format(&#39;%6.3f&#39;,lonv)<br>&nbsp;&nbsp;&nbsp;&nbsp;latval = math_format(&#39;%6.3f&#39;,latv)<br>&nbsp;&nbsp;&nbsp;&nbsp;rc = write(&#39;test_xy2w.txt&#39;,&#39;&#39;i&#39; &#39;j&#39; &#39;lonval&#39; &#39;latval&#39;&#39;)<br>&nbsp;&nbsp;&nbsp;&nbsp;j=j+1
<br>&nbsp;&nbsp;endwhile<br>&nbsp;&nbsp;i=i+1<br>endwhile<br>rc = close(&#39;test_xy2w.txt&#39;)<br><br>i j&nbsp;&nbsp;lonval latval<br>1 1 279.617 13.130<br>1 2 279.617 17.090<br>2 1 284.369 13.130<br>2 2 284.369 17.090<br><br></blockquote></div><br>
<br clear="all"><br>-- <br>Ing. Rafael E. Mundaray M.