In my case the correct map can be plotted by gradsc using a .ctl file with the line:<br><br>pdef 98 78 lcc 44.00 15.00 49.00 39.00 55.00 31.00 15.00 60000. 60000.<br><br>I know that for sure because the usual output of the model I use (RegCM3) is like that. This time I have problems because I used the PostProcess routines to interpolate the CRU database to the same grid and unfortunately the output is an .nc file. If I could somehow change the format of the file from .nc to the usual direct and unformatted fortran file I would be OK for now. Do you know how <br>I could do that?<br><br>Kostas<br><br>>Date: Mon, 7 Jan 2008 05:58:07 -0800<br>>From: brian vant-hull <<span style="border-bottom: 1px dashed rgb(0, 102, 204); background: transparent none repeat scroll 0% 50%; cursor: pointer; -moz-background-clip: -moz-initial; -moz-background-origin: -moz-initial; -moz-background-inline-policy: -moz-initial;" class="yshortcuts"
id="lw_1199787245_0">brianvanthull@YAHOO.COM</span>><br>>Subject: Re: wrong coordinates when plotting from an .nc file<br>><br>>Hi Costas;<br>>I hope to do something similar in the near future, but<br>>haven't started work on it yet. What I want to do is<br>>feed in the data plus somewhat irregular (but<br>>gridlike) lat and lon, and get out the data in a<br>>projection that GrADS will plot.<br>><br>>It looks like the projections described in the user<br>>manual *might* do this, but I haven't figured out how<br>>to feed in the original unprojected lat lon data. I<br>>fear the lat lon grids are expected to be already<br>>mercator projected. If so I'll have to use IDL or<br>>MatLab to warp my grid to the projection first, or<br>>write a rebinning program and perhaps lose some<br>>resolution.<br>><br>>Anyone else out there: do we need to warp our data to<br>>a given projection outside of
GrADS, or will GrADS do<br>>this for us? I see no provision for feeding arbitrary<br>>lat/lon grids into grads.<br>><br>>(unless plotting as station data - though it doesn't<br>>look like this will provide contoured output.)<br>><br>><br>>--- Costas Douvis <cdouvis@YAHOO.COM> wrote:<br>><br>>> Thank you for your response Brian<br>>><br>>> My grid has an interval of 50km. So when I use the<br>>> usual "rectilinear" map projection (where the<br>>> coordinates are latitude and longitude and a<br>>> "square" of 10x10 degrees does appear as a square)<br>>> my domain should look stretched towards the pole,<br>>> narrow towards the equator and curved. But it<br>>> doesn't: it appears rectangular<br>>><br>>> I'm running SUSE 10.0 - 32bit<br>>> and Grads version 1.8SL11<br>>><br>>> Please let me know if this information is still not<br>>>
enough<br>>><br>>> Kostas<br>>><br>>>> Date: Fri, 4 Jan 2008 07:18:23 -0800<br>>>> From: brian vant-hull<br>>>><br>>>> Subject: Re: wrong coordinates when plotting from an<br>>>> .nc file<br>>>><br>>>> You'll need to be more precise...for example if the<br>>>> coordinates are in equal increments of lat and lon<br>>>> then in an equal distance mapping they'd get<br>>>> narrower<br>>>> towards the poles, but the typical lat-lon plot<br>>>> would<br>>>> stretch it. You must have something different, but<br>>>> hard to tell what it is.<br>>>><br>>>> I may not be able to answer the question, but for<br>>>> someone else to answer it they'll need to know alot<br>>>> more details!<br>>>><br>>>> -Brian.<br>>>><br>>>> > ---
Costas Douvis wrote:<br>>>> ><br>>>> > Hello and happy new year to everyone<br>>>> ><br>>>> > I have problem plotting the data from an .nc file.<br>>>> I<br>>>> > believe that the values of the field are correct<br>>>> but<br>>>> > the coordinates are "stretched" so that they fit<br>>>> in<br>>>> > a rectangle. The real grid comes from an RCM model<br>>>> > and is getting narrower towards the equator.<br>>>> ><br>>>> > Thanks in advance for any help<br>>>> > Kostas Douvis<p> 
<hr size=1>Looking for last minute shopping deals? <a href="http://us.rd.yahoo.com/evt=51734/*http://tools.search.yahoo.com/newsearch/category.php?category=shopping">
Find them fast with Yahoo! Search.</a>