<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-GB link=blue vlink=purple><div class=WordSection1><p><span lang=EN>Djordje,<o:p></o:p></span></p><p><span lang=EN>I’m not sure, but I suspect that the problem might be related to the fact that the xdef entry in the ctl file appears to specify a grid that covers more than 360degrees of longitude, in that 145x2.5=362.5. So I think GrADS has failed to identify that the grid should be cyclic in longitude or just can’t work out which data to use in the “overlap” area. I certainly couldn’t get anything to plot across 0E when I hacked up an example, even when I tried a few standard “tricks”.<o:p></o:p></span></p><p><span lang=EN>(The fact that the grid is rather large in latitude too, >180deg, probably doesn’t matter because there’s no need for GrADS to wrap round in latitude. Also note that for a “basic” gridded data file, the t,z,y,x bit in the var line is not used, as far as I know – I tend to put “99” in there instead.)<o:p></o:p></span></p><p><span lang=EN>IF you only had one time level in the data file (which only has a single variable on a single level in it), you might be able to change the xdef line to have nx=144 so that GrADS thinks it only has data for 360deg – BUT THAT WILL NOT work for your current example which has >1 time.<o:p></o:p></span></p><p><span lang=EN>I think you’ll probably have to trim the grid to 360deg before going into GrADS if you want to plot across 0E.<o:p></o:p></span></p><p><span lang=EN>If the data file RTmavglobal(1961-1990).grd was written using fwrite to output a subset of a larger file, then I think you will be able to proceed by explicitly setting ‘set x 1 144’ before you write the file. Otherwise GrADS will likely have added an extra column of data because of how it normally sets up for a global grid.<o:p></o:p></span></p><p><span lang=EN>Doug<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b><span lang=EN-US style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span lang=EN-US style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> gradsusr-bounces@gradsusr.org [mailto:gradsusr-bounces@gradsusr.org] <b>On Behalf Of </b>Djordje Romanic<br><b>Sent:</b> 09 May 2011 08:54<br><b>To:</b> GrADS Users Forum<br><b>Subject:</b> [gradsusr] LON-LAT for European sector<o:p></o:p></span></p></div><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>Dear Users,<o:p></o:p></p></div><div><p class=MsoNormal> <o:p></o:p></p></div><div><p class=MsoNormal>I have a *.grd file (relative topography data) covering the whole globe. His CTL file is as follows: <o:p></o:p></p></div><div><p class=MsoNormal> <o:p></o:p></p></div><div><p class=MsoNormal>***************************************************************************<o:p></o:p></p></div><div><p><span lang=EN>dset ^RTmavglobal(1961-1990).grd<o:p></o:p></span></p><p><span lang=EN>undef 9.999e+20<o:p></o:p></span></p><p><span lang=EN>xdef 145 linear 0 2.5<o:p></o:p></span></p><p><span lang=EN>ydef 73 linear -90 2.5<o:p></o:p></span></p><p><span lang=EN>zdef 1 levels 1000<o:p></o:p></span></p><p><span lang=EN>tdef 43709 linear 0Z31jan1961 360mn<o:p></o:p></span></p><p><span lang=EN>vars 1<o:p></o:p></span></p><p><span lang=EN>RT 1 t,z,y,x **Relative topography 500/1000 [hPa]<o:p></o:p></span></p><p><span lang=EN>endvars<o:p></o:p></span></p><p><span lang=EN>**************************************************************************<o:p></o:p></span></p><p><span lang=EN>The data on global map (LAT -90 90; LON 0 360) are attached (global.jpeg). Now, from the global file I want to extract the data for the European sector (longitude interval is from 30W to 60E and latitude interval is from 0 90N), but after the following commands: <o:p></o:p></span></p><p><span lang=EN>***********************************************************************<o:p></o:p></span></p><p><span lang=EN>set lon -30 60<o:p></o:p></span></p><p><span lang=EN>set lat 0 90<o:p></o:p></span></p><p><span lang=EN>d rt<o:p></o:p></span></p><p><span lang=EN>************************************************************************<o:p></o:p></span></p><p><span lang=EN>I get the image that does not have the data East of the 0 meridian(attached as file euro.jpeg). <o:p></o:p></span></p><p><span lang=EN>Please, can someone help me with this issue? What is causing this? <o:p></o:p></span></p><p><span lang=EN>Kind regards,<o:p></o:p></span></p><p><span lang=EN>Djordje<o:p></o:p></span></p></div></div><br />--
<br />This message (and any attachments) is for the recipient only. NERC
<br />is subject to the Freedom of Information Act 2000 and the contents
<br />of this email and any reply you make may be disclosed by NERC unless
<br />it is exempt from release under the Act. Any material supplied to
<br />NERC may be stored in an electronic records management system.</body></html>