<div dir="ltr"><div>Not tested</div><div><br></div><div>Explicit approach<br></div><div><br></div><div>   prc1d(time)</div><div>   ntim = dimsizes(prc1d)</div><div>   nlat = 1</div><div>   mlon = 1</div><div><br></div><div>; expand 1D to 3D<br></div><div><br></div><div>   prc3d  = <b>conform_dims</b>((/ntim,nlat,mlon/), prc1d, 0)</div><div>   printVarSummary(prc3d)</div><div><br></div><div>    yyyyddd =</div><div><br></div><div>   prc3d_ClmDay    = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/clmDayTLL.shtml" target="_blank"><b>clmDayTLL</b></a>(prc3d, yyyyddd)</div><div>   printVarSummary(prc_3d)<br></div><div><br></div><div>   <br></div><div>   prc3d_ClmDay    = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/clmDayTLL.shtml" target="_blank"><b>clmDayTLL</b></a>(prc3d, yyyyddd)</div><div>   prc3d_Anom   = <b>calcDayAnomTLL</b> (prc3d, yyyyddd, prc3d_ClmDay)</div><div>    :</div><div>    :</div><div>======================</div><div>2nd approach, cleaner and easier</div><div><br></div><div>there are <i>undocumented functions</i>  <b>clmDayT</b>, <b>calcDayAnomT</b></div><div><br></div><div>prc1d_ClmDay    = <b>clmDayT</b>(prc1d, yyyyddd)</div><div>printVarSummary(prc1d_ClmDay )<br></div><div>prc1d_AnomDay =<b> calcDayAnomT</b>(prc1d, yyyyddd, prc1d_ClmDay)<b><br></b></div><div>printVarSummary(prc1d_AnomDay)</div><div><br></div><div>Look at <b>contributed.ncl</b></div><div> search for the above functions if you want<br></div><div><b></b></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Aug 6, 2021 at 12:20 PM Lyndz via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear NCL-Experts,</div><div><br></div><div>I have daily rainfall data from 1979-2017.</div><div><br></div><div>I would like to:</div><div><br></div><div>[a] Remove the annual cycle</div><div>[b] Then remove the seasonal cycle (3rd harmonics).</div><div><br></div><div>I saw this function in NCL which computes the daily climatology first before getting the anomalies (including smoothed anomalies). </div><div><br></div><div><b><a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/calcDayAnomTLL.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/calcDayAnomTLL.shtml</a></b></div><div><b><br></b></div><div>But this can only be applied for <b>time x lat x lon data</b>. </div><div><br></div><div>Is there an equivalent function of one-dimationsal (time series) data like this?</div><div>If none, do you have any suggestions on how I can do this in NCL?</div><div><br></div><div><br></div><div>I'll appreciate any help.</div><div><br></div><div>-Lyndz</div><div><br></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>