[gradsusr] Trend MAP help!

Andrew Friedman andfried at gmail.com
Thu Dec 17 15:59:53 EST 2015


Hi Malcolm,

If you are looking for the slope (the rate of change per time step) you need to add another term when you call trend

For example: 
'ltrend prec precl precslope'

>From the ltrend.gs”

ga-> ltrend
  Calculate linear trend over time (based on least-squares fitting).

  USAGE: ltrend <input> [<output> [<slope> [<rmse>]]]
    <input>: input field. Can be any GrADS expression.
    <output>: output field, i.e., fitted trend line. Default=<input>.
    <slope>: slope of fitted trend line, i.e., change of <input> per time step.
    <rmse>: root mean square error.

Andrew

> On Dec 17, 2015, at 1:28 PM, Malcolm Tang <malct32 at gmail.com> wrote:
> 
> Hi Andrew,
> Yes I've set it as
> 
> "set t 1 last"
> 
> For your info, I tried this and a spatial map came out:
> 
> 'set lon 0 360'
> 'set lat -90 -75'
> 'set mproj sps'
> 'set t 1 last'
> 
> 'd prec'
> 
> 'ltrend prec precl'
> 
> 'd precl'
> 
> I'm not sure if it's the spatial trend map that I'm looking for though....
> 
> But since the time is set from 1 to last, I can observe that the map actually changes from 1 to the last time step.
> 
> Any idea or suggestion?
> 
> Thanks.
> Malcolm
> 
> 
> On Thu, Dec 17, 2015 at 5:32 PM, Andrew Friedman <andfried at gmail.com> wrote:
> Hi Malcolm,
> 
> In your text earlier, it looks like your time is set to 1. Can you try setting time to varying?
> 
> Andrew
> 
> > On Dec 17, 2015, at 3:38 AM, Malcolm Tang <malct32 at gmail.com> wrote:
> >
> > Hi Andrew,
> >
> > I removed the bracket and ....... nothing came out?
> >
> > I typed:
> >
> > "ltrend prec fitted trend line" nothing came out.
> >
> > "ltrend prec fitted trend line 0.2 0.025" nothing came out.
> >
> > Should I be expecting a result or something? Because I'm expecting a spatial map to come up.
> >
> > Thanks!
> > Malcolm
> >
> > On Wed, Dec 16, 2015 at 4:05 PM, Andrew Friedman <andfried at gmail.com> wrote:
> > HI Malcolm,
> >
> > You shouldn’t type the brackets - they’re placeholders to indicate optional values.
> >
> > Andrew
> >
> > > On Dec 16, 2015, at 4:55 AM, Malcolm Tang <malct32 at gmail.com> wrote:
> > >
> > > Hi all,
> > >
> > > I want to plot the trend map of precipitation in antarctica (75s-90s) from 1999-2015.
> > >
> > > I have done some searching and saw someone suggested using BinGuan's script ltrend.gs
> > > I have copied the script and the dependencies qdims.gsf.
> > >
> > > So these are what I type in grads:
> > >
> > > 'set lon 0 360'
> > > 'set lat -90 -75'
> > > 'set mproj sps'
> > > 'set t 1'
> > >
> > > 'd prec'
> > > 'ltrend prec [fitted trend line [0.2[0.025]]]'
> > >
> > > and there will be error coming out: DEFINE error: Name too long; missing '='
> > >
> > > Anyone who has experience with ltrend.gs PLEASEEEEEEE help!!!
> > >
> > > Thanks
> > > Malcolm
> > > _______________________________________________
> > > 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
> >
> > _______________________________________________
> > 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
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr




More information about the gradsusr mailing list