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

[Andrew Friedman]

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