so should i change the Proj_units file to:<br>
<br>
unit: degree<br>
units: degrees<br>
meters:&nbsp; 111300 * cos(y())<br>
<br>
-vishal<br>
<br>
<br><div><span class="gmail_quote">On 8/25/06, <b class="gmail_sendername">Glynn Clements</b> &lt;<a href="mailto:glynn@gclements.plus.com">glynn@gclements.plus.com</a>&gt; wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>Vishal Mehta wrote:<br><br>&gt; &gt; This problem can be solved to an extent by using the value from the<br>&gt; &gt; PROJ_UNITS file to convert ground distances to metres (although that<br>&gt; &gt; requires that the PROJ_UNITS file actually contains the correct
<br>&gt; &gt; value).<br>&gt;<br>&gt; I did'nt fully understand this comment Glynn. Could you please elaborate<br>&gt; this idea?<br>&gt; FYI,<br>&gt; GRASS 6.0.1 (cru3):~/CRU &gt; g.proj -w<br>&gt; GEOGCS[&quot;wgs84&quot;,
<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; DATUM[&quot;WGS_1984&quot;,<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; SPHEROID[&quot;wgs84&quot;,6378137,298.257223563],<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; TOWGS84[0.000,0.000,0.000]],<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; PRIMEM[&quot;Greenwich&quot;,0],<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp; UNIT[&quot;degree&quot;,
0.0174532925199433]]<br>&gt;<br>&gt; And my Proj_units file in the PERMANENT mapset is:<br><br>&gt; unit: degree<br>&gt; units: degrees<br>&gt; meters: 1.0<br><br>The &quot;meters&quot; value in the PROJ_UNITS file is supposed to be a scale
<br>factor which converts map units to metres (i.e. 1.0 for metres, 0.3048<br>for feet). For degrees, it should be roughly 111300 * cos(latitude).<br><br>r.slope.aspect uses that value to convert map units to metres; if it's
<br>set correctly (and the zfactor= option is set correctly if your<br>elevations aren't in metres), your slope map should be correct. There<br>might still be problems with r.sun itself, though.<br><br>&gt; I'm not sure what the last line of 
g.proj -w means, what does unit there<br>&gt; refer to, as opposed to the proj_units 'meters' field?<br><br>0.0174532925199433 is pi/180, i.e. it represents the size of the<br>(angular) units in radians.<br><br>--<br>Glynn Clements &lt;
<a href="mailto:glynn@gclements.plus.com">glynn@gclements.plus.com</a>&gt;<br></blockquote></div><br>