<div dir="ltr">Hi,<div><br></div><div>The problem arises from the fact that two files are opened at same time and their time steps are different. Your first file starts from 1st Jan 2013. This means time step t=1 is 2013 01 01 and goes on increasing. But your second file starts from 1981. When your default file is 1, t=1 for second file will not be equal to for 1981 01 01. 't' value for 1981 01 01 will be something like </div>
<div>-(2013-1981)*365 (minus sign). I hope it is clear.</div><div><br></div><div>Either create a dummy file from 1981-2012 or separate 2013 from the climatology file. This should solve your issue. I hope.</div><div><br></div>
<div>-rupak</div><div><br></div><div>Rupak Rajbhandari</div><div>Department of Meteorology</div><div>Tri-Chandra College</div><div>Kathmandu, NEPAL</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Thu, Sep 19, 2013 at 7:03 AM, Jeff Duda <span dir="ltr"><<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<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"><div><div class="h5"><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><font color="#888888"><br>
ron<br>
</font></span><div><div><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" target="_blank">gradsusr-bounces@gradsusr.org</a><br>
[mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">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" target="_blank">gradsusr@gradsusr.org</a><br>
Subject: [gradsusr] Differencing values from two files in GrADS<br>
<br>
</div></div><div><div>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" target="_blank">philk@atmos.colostate.edu</a><br>
----------------------------------------------------<br>
<br>
<br>
</div></div><div><div>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">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></div></div><span class="HOEnZb"><font color="#888888">-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>
Center for Analysis and Prediction of Storms<br>
</font></span></div>
<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><br>
<br></blockquote></div><br></div>