[gradsusr] Serious Lterp Issue With GFS Data & Climatology Data

[Perry, Aaron @ LSC]

Hi Everyone,

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.

Here's the revised script:

'open 'ctl1''
'open 'ctl2''

'set mproj nps'
'set lat 'latmin' 'latmax''
'set lon 'lonmin' 'lonmax''
'set mpvals -130 -64 20 70'
'set display color white'
'clear'
'set map 1 1 10'
'set mpdset hires'
'set grid off'

'set grads off'
'set timelab off'
'set gxout shade2b'

'define terp2 = TMPprs.2(lev=850)'
'define terp1 = TMP850mb.1'

'define newvar = lterp(terp1(t=1),terp2(t=1))'

'd newvar'

Any ideas for getting this to work? It seems simple enough but...



Onward and Upward,

Aaron D. Perry
Boston, MA
Lyndon State College '16
Twitter: @AaronPerryLSC<https://twitter.com/AaronPerryLSC>
Mobile: 617-780-4312
________________________________
From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] on behalf of Jeff Duda [jeffduda319 at gmail.com]
Sent: Thursday, January 16, 2014 1:03 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] Serious Lterp Issue With GFS Data & Climatology Data

Remove the double quotes in your time strings and see if that makes any difference.


On Thu, Jan 16, 2014 at 7:00 AM, Andrew Friedman <andfried at berkeley.edu<mailto:andfried at berkeley.edu>> wrote:
Hi Aaron,

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.

Andrew

On Jan 16, 2014, at 3:26 AM, Perry, Aaron @ LSC <Aaron.Perry at lsc.vsc.edu<mailto:Aaron.Perry at lsc.vsc.edu>> wrote:

> Good Evening Everyone,
>
> 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...
>
> The script in question is below:
>
> ctl1 = '~/grads/data/cfsrclim/t850.cfsr.mean.clim.daily.1982.2010.ctl'
> ctl2 = '~/grads/data/gfs.t12z.pgrbf33.ctl'
>
> 'open 'ctl1''
> 'set x 1 360'
> 'open 'ctl2''
>
> 'set dfile 1'
> 'define t850clim = (TMP850mb-273.15)*(9/5)+32)'
>
> 'set dfile 2'
> 'define t850curr = (tmpprs(lev=850)-273.15)*(9/5)+32)'
>
> 'define t850a = (t850clim(t="00Z01jan1984"))'
>
> 'define t850b = (t850curr(t="18Z15jan2014"))'
>
> 'define final = lterp(t850b,t850a)'
>
> 'define hi = t850a-final'
>
> 'd hi'
>
>
> Onward and Upward,
>
> Aaron D. Perry
> Boston, MA
> Lyndon State College '16
> Twitter: @AaronPerryLSC
> Mobile: 617-780-4312
> <climtest.gs<http://climtest.gs>>_______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
> http://gradsusr.org/mailman/listinfo/gradsusr


_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr



--
Jeff Duda
Graduate research assistant
University of Oklahoma School of Meteorology
Center for Analysis and Prediction of Storms
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140116/d488f9f9/attachment.html