[ncl-talk] question about "dtrend" function
Dennis Shea
shea at ucar.edu
Fri Sep 5 11:33:36 MDT 2014
This question is a bit ambiguous.
Assuming y(time,lat,lon) ; (ntim,nlat,mlon)
yNew = dtrend_n(y, True, 0) ;
yNew=dtrend(y(lat|:,lon|:,time|:), True)
slope = *onedtond*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/onedtond.shtml>
(yNew at slope,(/nlat,nlon/)) ; slope(nlat,mlon)
copy_VarCoords(y(0,:,:), slope)
slope at long_name = "..."
slope at units = "..."
printVarSummary(slope)
; arithmetic stddev
slope_std = stddev(slope)
The weighted stddev could also be determined by weight the grid by (say)
cos(lat*rad)
On Sat, Aug 30, 2014 at 2:51 PM, Verena Lili <verena.prick at gmail.com> wrote:
> Hallo,
>
> Anybody try to have a standard deviation of the trend using "dtrend"
> function? I'd like to plot the trend values together with the
> their stdev. Thank you..,
>
> Cheers
>
>
> --
> Regards,
> Dr. Verena.
> School of Marine and Atmospheric Sciences
> Stony Brook University
> homepage: http://www.somas.stonybrook.edu/
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140905/d1425409/attachment.html
More information about the ncl-talk
mailing list