<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: arial,helvetica,sans-serif; font-size: 12pt; color: #000000'>Hi Aaron,<div><br></div><div>The anomaly plot looks realistic, and your script seems to be correct (there are a few redundant commands, but the result is OK).</div><div><br></div><div>Eric<br><br><hr id="zwchr"><b>From: </b>"Aaron @ LSC Perry" <Aaron.Perry@lsc.vsc.edu><br><b>To: </b>"GrADS Users Forum" <gradsusr@gradsusr.org><br><b>Sent: </b>Saturday, January 18, 2014 3:41:20 PM<br><b>Subject: </b>Re: [gradsusr] Serious Lterp Issue With GFS Data & Climatology Data<br><br>
<style>
<!--
p
{margin:0}
-->
</style><style id="owaParaStyle">P {margin-top:0;margin-bottom:0;}</style>
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Eric,<br>
<br>
Thank you very much for responding! I use xcbar to plot my colorbars.<br>
<br>
I've attached a sample of the anomaly graphic to this email, as well as the latest working script.<br>
<div><br>
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-family:Tahoma; font-size:13px"><font size="3"><font face="Arial"><font size="2" face="Tahoma"><font size="2"><br>
Onward and Upward,<br>
<br>
</font>Aaron D. Perry<br>
</font><font size="3"><font size="3"><font size="2" face="Tahoma">Boston, MA<br>
Lyndon State College '16<br>
Twitter: <a href="https://twitter.com/AaronPerryLSC" target="_blank">@AaronPerryLSC</a><br>
</font></font></font></font></font><font size="3"><font face="Arial"><font size="3"><font size="3"><font size="2" face="Tahoma">Mobile: 617-780-4312</font></font></font></font></font></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr>
<div style="direction: ltr;" id="divRpF775366"><font size="2" face="Tahoma" color="#000000"><b>From:</b> gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of Eric Altshuler [ela@cola.iges.org]<br>
<b>Sent:</b> Friday, January 17, 2014 11:56 PM<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] Serious Lterp Issue With GFS Data & Climatology Data<br>
</font><br>
</div>
<div></div>
<div>
<div style="font-family:arial,helvetica,sans-serif; font-size:12pt; color:#000000">
Hi Aaron,
<div><br>
</div>
<div>You need to include a color bar on your plot, showing numerical values for each color. Without it, there's no way to know quantitatively what the anomalies are. This grads script will make a color bar:</div>
<div><br>
</div>
<div>ftp://cola.gmu.edu/grads/scripts/cbarn.gs</div>
<div><br>
</div>
<div>Run the cbarn.gs script immediately after displaying the anomaly variable.</div>
<div><br>
</div>
<div>Best regards,</div>
<div><br>
<br>
<div><span></span>Eric L. Altshuler<br>
Research Scientist<br>
Center for Ocean-Land-Atmosphere Studies<br>
253 Research Hall, Mail Stop 6C5
<div>George Mason University</div>
<div>4400 University Drive<br>
Fairfax, VA 22030 USA<br>
<br>
E-mail: ela@cola.iges.org<br>
Phone: (703) 993-5725<br>
Fax: (703) 993-5770</div>
<span></span><br>
</div>
<br>
<hr id="zwchr">
<b>From: </b>"Aaron @ LSC Perry" <Aaron.Perry@lsc.vsc.edu><br>
<b>To: </b>"GrADS Users Forum" <gradsusr@gradsusr.org><br>
<b>Sent: </b>Friday, January 17, 2014 3:15:00 PM<br>
<b>Subject: </b>Re: [gradsusr] Serious Lterp Issue With GFS Data & Climatology Data<br>
<br>
Hi Jim,<br>
<br>
I was able to get the script to work with that you suggested! Thank you! This method seems to be the only method that really works. I've been experimenting with another method for getting the 850mb temp anoms but, I end up getting crazy library errors and then
it crashes GrADS.<br>
<br>
I attached an image of what I was able to generate for 850mb anoms. Does it look reasonable?<br>
<br>
<br>
Onward and Upward,<br>
<br>
Aaron D. Perry<br>
Boston, MA<br>
Lyndon State College '16<br>
Twitter: @AaronPerryLSC<br>
Mobile: 617-780-4312<br>
<br>
________________________________________<br>
From: gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of James T. Potemra [jimp@hawaii.edu]<br>
Sent: Friday, January 17, 2014 1:17 AM<br>
To: GrADS Users Forum<br>
Subject: Re: [gradsusr] Serious Lterp Issue With GFS Data & Climatology Data<br>
<br>
Hi Aaron,<br>
<br>
I'm a little confused as to the problem. What is error you are<br>
getting? When you load two different files you are<br>
potentially loading two different grids, including x/y/z/t. In these<br>
cases it's much easier to define lat/lon/z/t<br>
for each first. In your example below, you open one file, then the<br>
next, then define lon/lat. Try this instead:<br>
<br>
1. open first file<br>
2. define lat/lon/z/t range<br>
3. define new variable (in this case terp1)<br>
4. open second file<br>
5. set dfile 2<br>
6. define lat/lon/z/t range<br>
7. define new variable (in this case terp2)<br>
<br>
This will at least help narrow down whether you have mismatched grids<br>
and if so in what dimension. In your<br>
example file below the vertical dimension will default to the grid in<br>
the first file. I assume you want to interpolate<br>
lat/lon from one file to the lat/lon of the other? I don't think lterp<br>
works in the vertical.<br>
<br>
<br>
Jim<br>
<br>
On 1/16/14 1:10 PM, Perry, Aaron @ LSC wrote:<br>
> Good Evening Andrew,<br>
><br>
> I tried to incorporate the modify command into the script according to the modify example in the GrADS documentation but, no luck.<br>
><br>
><br>
> Onward and Upward,<br>
><br>
> Aaron D. Perry<br>
> Boston, MA<br>
> Lyndon State College '16<br>
> Twitter: @AaronPerryLSC<br>
> Mobile: 617-780-4312<br>
><br>
> ________________________________________<br>
> From: gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of Andrew Friedman [andfried@berkeley.edu]<br>
> Sent: Thursday, January 16, 2014 5:55 PM<br>
> To: GrADS Users Forum<br>
> Subject: Re: [gradsusr] Serious Lterp Issue With GFS Data & Climatology Data<br>
><br>
> Hi Aaron,<br>
> Have you tried the ‘modify’ command for making the climatology apply to the whole record?<br>
> Andrew<br>
><br>
> On Jan 16, 2014, at 7:57 PM, Perry, Aaron @ LSC <Aaron.Perry@lsc.vsc.edu> wrote:<br>
><br>
>> Hi Everyone,<br>
>><br>
>> I tried setting the time dimensions to the same value (t=1) and it looks like the climo data is showing but, now I need to subtract the operational model data from the climo data to get the anomalies.<br>
>><br>
>> Here's the revised script:<br>
>><br>
>> 'open 'ctl1''<br>
>> 'open 'ctl2''<br>
>><br>
>> 'set mproj nps'<br>
>> 'set lat 'latmin' 'latmax''<br>
>> 'set lon 'lonmin' 'lonmax''<br>
>> 'set mpvals -130 -64 20 70'<br>
>> 'set display color white'<br>
>> 'clear'<br>
>> 'set map 1 1 10'<br>
>> 'set mpdset hires'<br>
>> 'set grid off'<br>
>><br>
>> 'set grads off'<br>
>> 'set timelab off'<br>
>> 'set gxout shade2b'<br>
>><br>
>> 'define terp2 = TMPprs.2(lev=850)'<br>
>> 'define terp1 = TMP850mb.1'<br>
>><br>
>> 'define newvar = lterp(terp1(t=1),terp2(t=1))'<br>
>><br>
>> 'd newvar'<br>
>><br>
>> Any ideas for getting this to work? It seems simple enough but...<br>
>><br>
>><br>
>><br>
>> Onward and Upward,<br>
>><br>
>> Aaron D. Perry<br>
>> Boston, MA<br>
>> Lyndon State College '16<br>
>> Twitter: @AaronPerryLSC<br>
>> Mobile: 617-780-4312<br>
>> From: gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of Jeff Duda [jeffduda319@gmail.com]<br>
>> Sent: Thursday, January 16, 2014 1:03 PM<br>
>> To: GrADS Users Forum<br>
>> Subject: Re: [gradsusr] Serious Lterp Issue With GFS Data & Climatology Data<br>
>><br>
>> Remove the double quotes in your time strings and see if that makes any difference.<br>
>><br>
>><br>
>> On Thu, Jan 16, 2014 at 7:00 AM, Andrew Friedman <andfried@berkeley.edu> wrote:<br>
>> Hi Aaron,<br>
>><br>
>> I think your time indices need to be the same for lterp to work. Check out examples 4 and 5 in the documentation entry for lterp, which show how to manually override the time values. Hope this helps.<br>
>><br>
>> Andrew<br>
>><br>
>> On Jan 16, 2014, at 3:26 AM, Perry, Aaron @ LSC <Aaron.Perry@lsc.vsc.edu> wrote:<br>
>><br>
>>> Good Evening Everyone,<br>
>>><br>
>>> Currently, I'm having an issue with lterping two variables, one being 850mb temp with a recent GFS .ctl file and a climatological 850mb temp .ctl. The recent GFS .ctl is fine and the climate .ctl is coarse. I have been trying to read up on lterp as much
as possible but, I can't seem to get this work and I've reached the end of possible resources it seems like...<br>
>>><br>
>>> The script in question is below:<br>
>>><br>
>>> ctl1 = '~/grads/data/cfsrclim/t850.cfsr.mean.clim.daily.1982.2010.ctl'<br>
>>> ctl2 = '~/grads/data/gfs.t12z.pgrbf33.ctl'<br>
>>><br>
>>> 'open 'ctl1''<br>
>>> 'set x 1 360'<br>
>>> 'open 'ctl2''<br>
>>><br>
>>> 'set dfile 1'<br>
>>> 'define t850clim = (TMP850mb-273.15)*(9/5)+32)'<br>
>>><br>
>>> 'set dfile 2'<br>
>>> 'define t850curr = (tmpprs(lev=850)-273.15)*(9/5)+32)'<br>
>>><br>
>>> 'define t850a = (t850clim(t="00Z01jan1984"))'<br>
>>><br>
>>> 'define t850b = (t850curr(t="18Z15jan2014"))'<br>
>>><br>
>>> 'define final = lterp(t850b,t850a)'<br>
>>><br>
>>> 'define hi = t850a-final'<br>
>>><br>
>>> 'd hi'<br>
>>><br>
>>><br>
>>> Onward and Upward,<br>
>>><br>
>>> Aaron D. Perry<br>
>>> Boston, MA<br>
>>> Lyndon State College '16<br>
>>> Twitter: @AaronPerryLSC<br>
>>> Mobile: 617-780-4312<br>
>>> <climtest.gs>_______________________________________________<br>
>>> gradsusr mailing list<br>
>>> gradsusr@gradsusr.org<br>
>>> http://gradsusr.org/mailman/listinfo/gradsusr<br>
>><br>
>> _______________________________________________<br>
>> gradsusr mailing list<br>
>> gradsusr@gradsusr.org<br>
>> http://gradsusr.org/mailman/listinfo/gradsusr<br>
>><br>
>><br>
>><br>
>> --<br>
>> Jeff Duda<br>
>> Graduate research assistant<br>
>> University of Oklahoma School of Meteorology<br>
>> Center for Analysis and Prediction of Storms<br>
>> _______________________________________________<br>
>> gradsusr mailing list<br>
>> gradsusr@gradsusr.org<br>
>> http://gradsusr.org/mailman/listinfo/gradsusr<br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> gradsusr@gradsusr.org<br>
> http://gradsusr.org/mailman/listinfo/gradsusr<br>
><br>
> _______________________________________________<br>
> gradsusr mailing list<br>
> gradsusr@gradsusr.org<br>
> http://gradsusr.org/mailman/listinfo/gradsusr<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
gradsusr@gradsusr.org<br>
http://gradsusr.org/mailman/listinfo/gradsusr<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
gradsusr@gradsusr.org<br>
http://gradsusr.org/mailman/listinfo/gradsusr<br>
</div>
</div>
</div>
</div>
</div>
<br>_______________________________________________<br>gradsusr mailing list<br>gradsusr@gradsusr.org<br>http://gradsusr.org/mailman/listinfo/gradsusr<br></div></div></body></html>