In my case the correct map (with the correct coodrinates) 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>But in that case the data file must be the usual fortran unformatted and direct output file. 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.<br><br>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 I could do that?<br><br>Brian, in case you want to use this method be careful. If I remember correctly there is a number of ways to describe your grid and "lcc" (in the pdef line) is only one of them<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_1199788065_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<br><br><br><br><br><br><p> 
<hr size=1>Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. <a href="http://us.rd.yahoo.com/evt=51733/*http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ "> Try it now.</a>