<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 class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>Djordje,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>I haven’t used a GrADS-DODS server in recent years (shocking, I know) and so don’t know quite how it works in detail...but it looks to me as if the command you give will extract the data and I assume that you then write it to the file </span><span lang=EN style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>RTmavglobal(1961-1990).grd using further GrADS commands (e.g. ‘set gxout fwrite’ then “display “). If so, I assume you can write a subset of the data by using ‘set lon’ or ‘set x’ before using “display”. If the actual process you use to write the file is very different then, sorry, I don’t know what to do.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN style='font-size:11.0pt;font-family:"Calibri","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>(And I guess that your expression returns the same values regardless of whether you use 360 or 357.5 because in each case GrADS (or more likely the DODS server in this case....I told you I’m not sure of details!) decides that this covers the gridpoints at all longitudes.)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN style='font-size:11.0pt;font-family:"Calibri","sans-serif"'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN style='font-size:11.0pt;font-family:"Calibri","sans-serif"'>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><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 11:21<br><b>To:</b> GrADS Users Forum<br><b>Subject:</b> Re: [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 Doug, <o:p></o:p></p></div><div><p class=MsoNormal> <o:p></o:p></p></div><div><p class=MsoNormal>Now I see... I made a slip during the calculation. Namely, I've used NCEP/NCAR reanalysis data to calculate relative topography, and to do that I've used the command: <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 class=MsoNormal>sdfopen <a href="http://nomad1.ncep.noaa.gov:9090/dods/_expr_%7breanalyses/reanalysis-1/6hr/pgb/pgb%7d%7bhgtprs(lev=500)-hgtprs(lev=1000)%7d%7b0">http://nomad1.ncep.noaa.gov:9090/dods/_expr_{reanalyses/reanalysis-1/6hr/pgb/pgb}{hgtprs(lev=500)-hgtprs(lev=1000)}{0</a>:<span style='background:#CC0000'>360</span>,-90:90,1000:1000,00z1jan1961:00z1jan1991}<o:p></o:p></p></div><div><p class=MsoNormal>***************************<o:p></o:p></p></div><div><p class=MsoNormal>However, the longitude in the original file (on nomad server) goes only up to 357.5. BUT, even if I replace 360 with 357.5 in the above command I get the same result---WHY? <o:p></o:p></p></div><div><p class=MsoNormal>Dear Doug, since I have more than a month of work behind me and all of that work was done usung the above command (with lon=360) and the result is the same as it would be 357.5 instead of 360, but the difference is in writing the result in a file, is it possible to change something in the CTL file and to correct a mistake? You were right, if in the previous command is 357.5, than I am able to set lon from -30 to 60. Since eveything else is the same, except this ploting issue, please is it possible to correct that so that I don't have to do the same tremandeus work from the beginning? <o:p></o:p></p></div><div><p class=MsoNormal> <o:p></o:p></p></div><div><p class=MsoNormal>Regards,<o:p></o:p></p></div><div><p class=MsoNormal>Djordje<o:p></o:p></p></div><div><p class=MsoNormal><br><br> <o:p></o:p></p></div><div><p class=MsoNormal>2011/5/9 Clark, Douglas B. <<a href="mailto:dbcl@ceh.ac.uk">dbcl@ceh.ac.uk</a>><o:p></o:p></p><div><div><p><span lang=EN>Djordje,</span><o:p></o:p></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”.</span><o:p></o:p></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.)</span><o:p></o:p></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.</span><o:p></o:p></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.</span><o:p></o:p></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.</span><o:p></o:p></p><p><span lang=EN>Doug</span><o:p></o:p></p><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><span style='font-size:11.0pt;color:#1F497D'> </span><o:p></o:p></p><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'><b><span lang=EN-US style='font-size:10.0pt'>From:</span></b><span lang=EN-US style='font-size:10.0pt'> <a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>] <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</span><o:p></o:p></p></div><div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>Dear Users,<o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>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 style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'> <o:p></o:p></p></div><div><p class=MsoNormal style='mso-margin-top-alt:auto;mso-margin-bottom-alt:auto'>***************************************************************************<o:p></o:p></p></div><div><p><span lang=EN>dset ^RTmavglobal(1961-1990).grd</span><o:p></o:p></p><p><span lang=EN>undef 9.999e+20</span><o:p></o:p></p><p><span lang=EN>xdef 145 linear 0 2.5</span><o:p></o:p></p><p><span lang=EN>ydef 73 linear -90 2.5</span><o:p></o:p></p><p><span lang=EN>zdef 1 levels 1000</span><o:p></o:p></p><p><span lang=EN>tdef 43709 linear 0Z31jan1961 360mn</span><o:p></o:p></p><p><span lang=EN>vars 1</span><o:p></o:p></p><p><span lang=EN>RT 1 t,z,y,x **Relative topography 500/1000 [hPa]</span><o:p></o:p></p><p><span lang=EN>endvars</span><o:p></o:p></p><p><span lang=EN>**************************************************************************</span><o:p></o:p></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: </span><o:p></o:p></p><p><span lang=EN>***********************************************************************</span><o:p></o:p></p><p><span lang=EN>set lon -30 60</span><o:p></o:p></p><p><span lang=EN>set lat 0 90</span><o:p></o:p></p><p><span lang=EN>d rt</span><o:p></o:p></p><p><span lang=EN>************************************************************************</span><o:p></o:p></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). </span><o:p></o:p></p><p><span lang=EN>Please, can someone help me with this issue? What is causing this? </span><o:p></o:p></p><p><span lang=EN>Kind regards,</span><o:p></o:p></p><p><span lang=EN>Djordje</span><o:p></o:p></p></div></div></div></div><p class=MsoNormal><span style='color:#888888'><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.</span><o:p></o:p></p></div><p class=MsoNormal style='margin-bottom:12.0pt'><br>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><o:p></o:p></p></div><p class=MsoNormal><o:p> </o:p></p></div></body></html>