<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">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>
<div>Any ideas for getting this to work? It seems simple enough but...<br>
<br>
<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">@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 tabindex="-1">
<div style="direction: ltr;" id="divRpF118301"><font size="2" color="#000000" face="Tahoma"><b>From:</b> gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] on behalf of Jeff Duda [jeffduda319@gmail.com]<br>
<b>Sent:</b> Thursday, January 16, 2014 1:03 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 dir="ltr">Remove the double quotes in your time strings and see if that makes any difference.<br>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Thu, Jan 16, 2014 at 7:00 AM, Andrew Friedman <span dir="ltr">
<<a href="mailto:andfried@berkeley.edu" target="_blank">andfried@berkeley.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex; border-left:1px #ccc solid; padding-left:1ex">
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>
<div>
<div class="h5"><br>
On Jan 16, 2014, at 3:26 AM, Perry, Aaron @ LSC <<a href="mailto:Aaron.Perry@lsc.vsc.edu" target="_blank">Aaron.Perry@lsc.vsc.edu</a>> 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>
</div>
</div>
> <<a href="http://climtest.gs" target="_blank">climtest.gs</a>>_______________________________________________<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>
<br>
<br>
_______________________________________________<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>
</blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
Jeff Duda<br>
Graduate research assistant<br>
University of Oklahoma School of Meteorology<br>
Center for Analysis and Prediction of Storms<br>
</div>
</div>
</div>
</div>
</body>
</html>