<div dir="ltr"><div>There is no magic solution.<br><br>Daily data <b>[rf</b>]:<br><br>Continuous variables like temperature, pressure, etc have a *high* degree of day-to-day correlation.<br>Typically, there is 4-to-7 days between independent estimates. Hence, interpolating *short*<br>spans of missing daily values may be acceptable.<br><br>Unfortunately, daily precipitation has *low* day-to-day autocorrelation. Therefore, interpolating between <br>even short spans of missing daily precipitation values may not be appropriate. <br><br>Raw daily climatology: <b>xClmDay</b><br><br>Again, you could interpolate and you could get 'numbers' but you will have to decide how to proceed.<br><br>If there are large daily gaps and you do decide to interpolate the <b>xClmDay </b>variable, the <b>smthClmDay </b>result will be<br>affected by the interpolated values. Again, this is your decision.</div><div><br></div><div>You might look at <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/linmsg_n.shtml" target="_blank"><b>linmsg_n</b></a> to fill in missing values and <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/runave_n_Wrap.shtml" target="_blank"><b>runave_n_Wrap</b></a> for smoothing <br></div><div><br></div><div><b></b></div><div><b>Stay Healthy</b><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, May 2, 2020 at 1:36 PM Dr. Suman Maity via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@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>Dear all</div><div>This is Suman. I have 10 years daily data of rainfall. I want to calculate the anomaly using smoothed climatology. Here is my code:</div><div>load "./Extra.ncl"<br><i><b>begin<br> ymdStrt = 19940101 <br> ymdLast = 19981231 <br> nhar = 3 <br> fn = addfile ("./.<a href="http://rain.daily.nc" target="_blank">rain.daily.nc</a>", "r")<br> time = fn->time<br> ymd = cd_calendar(time, -2) <br> iStrt = ind(ymd.eq.ymdStrt) <br> iLast = ind(ymd.eq.ymdLast)<br> delete(time)<br> delete(ymd)<br>;***********************************************************<br>; Read user specified time and create required yyyyddd <br>;***********************************************************<br> time = fn->time(iStrt:iLast)<br> TIME = cd_calendar(time, 0)<br> year = floattointeger( TIME(:,0) )<br> month = floattointeger( TIME(:,1) )<br> day = floattointeger( TIME(:,2) ) <br> ddd = day_of_year(year, month, day) <br> yyyyddd = year*1000 + ddd<br><br> x = fn->rf(iStrt:iLast,:,:)<br> xClmDay = clmDayTLL(x, yyyyddd)<br> xClmDay_sm = smthClmDayTLLS(xClmDay, nhar)<br> printVarSummary(xClmDay_sm)</b></i></div><div><i><b>end</b></i></div><div><br></div><div>Now during the execution of the code, it shows "warning:ezfftf: 12451 input arrays contained missing values. No calculations performed on these arrays".<br></div><div>It seems to mean that "ezfftf" is not running because of "missing values". Yes my data contains missing values.</div><div><br></div><div>Therefore my query:</div><div>Is it not possible to find smoothed anomaly for a data with missing value using the above method? What is other ways in ncl to approach the issue?</div><div>Please answer me. Any sort of help/suggestions are highly appreciated.</div><div>Thanking you,</div><div><br></div><div>Best</div><div>Suman</div><div><div><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div>++++++++++++++++++++++++++++++++++++++++++</div><div><font size="2">Dr. Suman Maity</font></div><div><font size="2">Research Associate</font></div><div><font size="2">School of Atmospheric Science</font></div><div><font size="2">Sun Yat-sen University, Zhuhai campus,Tangjiawan Town<br></font></div><div><font size="2"> Zhuhai, Guangdong, Postcode-<span>519082,</span>China.</font></div><div><font size="2">Phone:+86-13543861045 (China), +91-9732636778 (India)<br></font></div><div><font size="2">skype: reach2suman</font></div><div><font size="2"><a href="https://www.researchgate.net/profile/Suman_Maity" target="_blank">https://www.researchgate.net/profile/Suman_Maity</a></font></div><span>++++++++++++++++++++++++++++++++++++++++++</span></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>