<div dir="ltr">There appears to be a dimension mismatch between your two files if you're getting that error message. Post the PDEF, XDEF, YDEF, ZDEF, and TDEF (and if it exists, EDEF) lines of each control file for comparison. We can help you with that.<div>
<br></div><div>Jeff Duda</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Sep 18, 2013 at 6:23 PM, Goodson,Ron [Edm] <span dir="ltr"><<a href="mailto:Ron.Goodson@ec.gc.ca" target="_blank">Ron.Goodson@ec.gc.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
slight correction ... meant to say "uave2" below. But .. I guess in<br>
your case it doesn't really matter since you a not comparing the two<br>
files<br>
<br>
(and I see you have d uave in both sections .. but assume that was a<br>
typo?)<br>
<span class="HOEnZb"><font color="#888888"><br>
ron<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
<br>
-----Original Message-----<br>
From: Goodson,Ron [Edm]<br>
Sent: September 18, 2013 5:15 PM<br>
To: 'GrADS Users Forum'<br>
Subject: RE: [gradsusr] Differencing values from two files in GrADS<br>
<br>
When I'm fiddling with two files with different times (or time steps),<br>
and I get confused by var.1(t=XX) or var.2(t=YY) at least, for debugging<br>
purposes, I will explicitly do a<br>
<br>
set dfile 1<br>
<br>
.. work with file 1<br>
<br>
set dfile 2<br>
<br>
.. work with file 2<br>
<br>
This way I know what t = 1, N really stands for. So, in your case I'd<br>
probably open both file then<br>
<br>
set dfile 2<br>
.. calculate uave<br>
<br>
set dfile 1<br>
.. compare variables in dfile 1 to uave<br>
<br>
<br>
ron<br>
<br>
<br>
<br>
-----Original Message-----<br>
From: <a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a><br>
[mailto:<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] On Behalf Of Phil Klotzbach<br>
Sent: September 18, 2013 4:57 PM<br>
To: <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
Subject: [gradsusr] Differencing values from two files in GrADS<br>
<br>
</div></div><div class="HOEnZb"><div class="h5">Hello,<br>
<br>
I've made a bit of progress since yesterday, but I'm still running into<br>
trouble. Effectively, the problem seems that I am trying to manipulate<br>
two files at the same time, but I'm not sure how to get around this,<br>
given that I have values for 2013 as well as a 30-year long-term<br>
climatology:<br>
<br>
What I am trying to do is as follows, generate a height/time<br>
cross-section of zonal wind anomalies from June 1 - September 15, 2013.<br>
<br>
I have two files that each have the same lat/lon grids. I can<br>
manipulate each file individually (with no other files open) and<br>
generate a figure, but what is giving me trouble is doing these<br>
calculations when both files are open at the same time. What I want to<br>
do is basically difference the 2013 zonal wind from the zonal wind<br>
climatology to get the anomaly. To generate the attached figures (with<br>
only one file open), here are the commands that I used:<br>
<br>
For 2013 Zonal Wind File:<br>
<br>
sdfopen <a href="http://uwnd.2013.nc" target="_blank">uwnd.2013.nc</a><br>
set t 152 258<br>
set z 1 16<br>
define uave = aave (uwnd, lon=-70, lon=-20, lat=0, lat=30) set x 1 set y<br>
1 d uave draw ylab height draw xlab date draw title 2013 Zonal Wind<br>
printim c:\2013uwind.jpg<br>
<br>
For 1981-2010 Climatology<br>
<br>
sdfopen <a href="http://uwnd.day.1981-2010.ltm.nc" target="_blank">uwnd.day.1981-2010.ltm.nc</a><br>
set t 152 258<br>
set z 1 16<br>
define uave2 = aave (uwnd.2, lon=-70, lon=-20, lat=0, lat=30) set x 1<br>
set y 1 d uave draw ylab height draw xlab date draw title 1981-2010<br>
Zonal Wind Climatology printim c:\zonalwindclimatology.jpg<br>
<br>
It seems like the problem arises when I try to do the same calculations<br>
with both files open at the same time. Even if I explicitly specify the<br>
t and z dimensions for both files before making the calculations, it<br>
appears that GrADS forgets(?) the dimensions in the file not being used<br>
for the specific calculation. For example, if I do the calculation for<br>
2013 first and then attempt the long-term climatology, the first<br>
calculation will be done correctly, and it will tell me that the second<br>
request is completely outside of file limits. If I change the active<br>
file and then make the calculation, the value in the first variable is<br>
lost.<br>
<br>
Any ideas? Thanks so much!<br>
<br>
Phil<br>
<br>
--<br>
----------------------------------------------------<br>
Phil Klotzbach, Ph.D.<br>
Research Scientist<br>
Department of Atmospheric Science<br>
Colorado State University<br>
Email: <a href="mailto:philk@atmos.colostate.edu">philk@atmos.colostate.edu</a><br>
----------------------------------------------------<br>
<br>
<br>
</div></div><div class="HOEnZb"><div class="h5">_______________________________________________<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><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
</div>