Dear Muhammad,<div><br></div><div>I am sorry it took me so long to answer. Thank you for your answer on point 2. </div><div><br></div><div>Regarding point 1, opening .nc files. I have attached some files. </div><div>The images from file1.jpg and file2.jpg are grads images (printim) from nc files produced with RegCM4, one is gxout contour and the other is gxout shaded.</div>
<div>The image from file3.jpg is a grads image (printim) gxout shaded for an nc file from the ESRL database. </div><div>The last image was produced with the PANOPLY data viewer from the same file that produced file3.jpg</div>
<div><br></div><div>As you can see the data is not being plotted with grads as with normal .ctl files (which are working properly).</div><div>I cannot seem to understand why it is not reading it properly. Maybe there are some definitions that I have to set? How should I check if all the settings are correct?</div>
<div><br></div><div>Regarding point 3, I think I can do them as a normal mathematical expression, but I'll be trying that later on. </div><div><br></div><div>Regards,</div><div><br></div><div>James<br><br><div class="gmail_quote">
On 11 December 2010 12:24, James Ciarlo` <span dir="ltr"><<a href="mailto:james.ciarlo@physics.org">james.ciarlo@physics.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Dear Muhammad,<div><br></div><div>I just had a look at it this morning. Unfortunately I haven't yet had the time to put it to use. Thank you very much for your help. I will get back to you about it for sure. Thanks again</div>
<div><br></div><div>Regards,</div><div><br></div><div><font color="#888888">James</font><div><div></div><div class="h5"><br><br><div class="gmail_quote">On 11 December 2010 09:59, Muhammad Rahiz <span dir="ltr"><<a href="mailto:muhammad.rahiz@ouce.ox.ac.uk" target="_blank">muhammad.rahiz@ouce.ox.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi James,<br>
I was wondering if you got my reply to the query you posted because I'm not sure if it gets delivered successfully.<br>
<br>
Anyway, if you did not, these are what I suggested<br>
<br>
1. display error<br>
- Try set gxout to something other than the default e.g. set gxout grfill<br>
- set clevs<br>
- check units<br>
<br>
2. print to ascii<br>
what you need is the <a href="http://fprintf.gs" target="_blank">fprintf.gs</a> script available at <a href="http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file" target="_blank">http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file</a><br>
<br>
Copy the code and name the file as <a href="http://fprintf.gs" target="_blank">fprintf.gs</a>. Then do the following,<br>
<br>
sdfopen <a href="http://input.nc" target="_blank">input.nc</a><br>
d var # where var = variable fprintf var dat.txt # prints data of var<br>
fprintf lon lon.txt # prints corresponding lon<br>
fprintf lat lat.txt # prints corresponding lat<br>
<br>
3. statistical analysis<br>
I know how to do so in R. If you're interested in using R for this purpose, I can help.<div><br>
<br>
HTH,<br>
<br>
--<br>
Muhammad Rahiz<br>
Researcher & DPhil Candidate (Climate Systems & Policy)<br>
School of Geography & the Environment<br>
University of Oxford<br>
<br>
On Fri, 10 Dec 2010, James Ciarlo` wrote:<br>
<br>
</div><div><div></div><div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear all,<br>
I am using GrADS 2.0.a9 on LINUX. I am trying to compare measured and modelled data by using GrADS. However since I am new with the program, there are<br>
some issues that I cannot seem to figure out, or find a solution to understand.<br>
I am experiencing:<br>
<br>
1) Difficulty with displaying data from NetCDF files<br>
I have opened the file with the command 'sdfopen'<br>
When I tried 'd prmsl' I got the following error:<br>
Operation error: Invalid dimension environment<br>
Min longitude > max longitude: 320 144<br>
So I used ' set lon -40 145' and tried the display command again<br>
I got the following line and the image attached (ncimage.jpg):<br>
Contouring: 99000 to 102300 interval 300<br>
And I am getting similar results with other nc files<br>
<br>
2) Difficulty with finding a command to output text files, or dat, or anything similar.<br>
I have used set lat, set lon, and set lev to obtain one value, and set t to get a long time interval, I have used d psa to display the time series, I<br>
can use printim file.jpg to get an image, but I cannot find a command to produce a file.txt so that I can make use of the actual numbers for<br>
statistical purposes.<br>
<br>
3) Difficulty with finding a command to work out comparative statistics, particularly correlation and significance.<br>
<br>
Regards,<br>
<br>
James Ciarlo`<br>
<br>
</blockquote>
</div></div></blockquote></div><br></div></div></div>
</blockquote></div><br></div>