<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Thank you very much Dennis, your suggestion works.<br><br>Hyacinth<div><br><div><div id="SkyDrivePlaceholder"></div>> Date: Mon, 22 Oct 2012 23:09:21 -0600<br>> From: shea@ucar.edu<br>> To: hyacinth.1@hotmail.com<br>> CC: ncl-talk@ucar.edu<br>> Subject: Re: [ncl-talk] Looping through dtrend function<br>> <br>> [1] dtrend (dtrend_n) returns anomalies. It removes the mean.<br>> <br>> [2] Is 'x' one dimensional<br>> x(0)=Jan, x(1)=Feb, ... x(11)=Dec, x(12)=Jan_next year , ...<br>> <br>> x_dtrnd = x ; create array including meta data<br>> do nmo=0,nmos-1<br>> work = x(nmo::12) ; every 12-th value<br>> wavg = avg(work) ; average<br>> work = dtrnd(work, False) ; anomalies<br>> work = wavg + work ; add average back in<br>> x_dtrend(nmo::nmos) = work<br>> end do<br>> <br>> <br>> <br>> On 10/22/12 9:20 PM, Hyacinth Nnamchi
wrote:<br>> > Hi;<br>> ><br>> > I am trying to detrend a monthly time series, month by month:<br>> ><br>> > nmos = 12<br>> > ntim = dimsizes(x)<br>> ><br>> > do n = 0,ntim-1,nmos<br>> > x_dtrnd = dtrend(x,False)<br>> > end do<br>> ><br>> > print(x_dtrnd(0:11))<br>> ><br>> > But the result doesn't look correct, looks like deviations from annual<br>> > mean. What could be wrong?<br>> ><br>> > Hyacinth<br>> ><br>> ><br>> ><br>> ><br>> ><br>> ><br>> > _______________________________________________<br>> > ncl-talk mailing list<br>> > List instructions, subscriber options, unsubscribe:<br>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk<br>> ><br></div><embed type="application/x-dap-background" id="DAPPlugin" style="visibility: collapse"></div> </div></body>
</html>