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

[Andrew Friedman]

Hi Aaron,
Have you tried the ‘modify’ command for making the climatology apply to the whole record?
Andrew

On Jan 16, 2014, at 7:57 PM, Perry, Aaron @ LSC <Aaron.Perry at lsc.vsc.edu> wrote:

> 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
> 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> 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> 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>_______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> 
> 
> _______________________________________________
> gradsusr mailing list
> 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
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr