On Thu, Nov 24, 2011 at 3:31 AM, Francesco Benincasa <span dir="ltr">&lt;<a href="mailto:francesco.benincasa@bsc.es">francesco.benincasa@bsc.es</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi all,<br>
I&#39;m using pygrads to draw some maps from netcdf files, </blockquote><div><br></div><div>Can you give us some more details here? Is this a global or regional dataset? Also, which version of pygrasd is this? Platform?</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">and I&#39;ve noticed that<br>
in some cases the map is drawn correctly only if I use the &quot;transf&quot; option set<br>
to True (the projection is cylindric).<br>
<br></blockquote><div><br></div><div>In principle, you should not have to set the &quot;transf&quot; attribute. This is an internal parameter to signal the package that a map transformation is necessary, &quot;under the hood&quot; kind of stuff.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
I&#39;ve taken a look at the code in _contourf method in the galab.py module:<br>
<br>
#       Evaluate GrADS expression<br>
#       -------------------------<br>
        if self.transf: self.cmd(&#39;set lon -180 180&#39;,Quiet=True)<br>
        Z = self.exp(expr)<br>
<br>
        g = Z.grid<br>
        if self.transf: self.cmd(&#39;set x %s %s&#39;%dh.x,Quiet=True)<br>
<br>
<br>
But I don&#39;t understand why it does that and why in some models it needs<br>
that options.<br>
<br>
Any idea?<br><br></blockquote><div><br></div><div>Here is the main issue that drives this special handling. GrADS is very flexlible in its addressing of longitudes: [-180,180],  [0,360], [20,380], all is fair game. The Basemap toolkit (and the proj4 library it relies on) less so, it wants longitudes in the range [-180,180]. Therefore, anytime a map transformation is necessary the longitude domain of the input grid get changed to the range [-180,180], a map transform is applied, and the the original grid arrangement is restored after the transform.</div>
<div><br></div><div>  I am not exactly clear what your particular problem is. If you send me a small sample file I can try it out and fix whatever bugs may be in there. I always work with global files, so it is possible that I may have overlooked some important detail for regional files (if indeed you have a regional file).</div>
<div><br></div><div>   Arlindo</div><div><br></div><div><br></div></div>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a><br>