<div dir="auto">One can use splithour in CDO.<div dir="auto"><br></div><div dir="auto">Soma</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Aug 12, 2019, 19:54 Dennis Shea via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><b>NCL: <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml" target="_blank" rel="noreferrer">calculate_monthly_values</a><br></b></div><div>----<br></div><div><b>T</b>he CDO has a <b>monmean </b>operator but I am not sure how to extract 6-hrly subsets.<br></div><div><b></b></div><div><b>CDO:</b> <a href="https://code.mpimet.mpg.de/projects/cdo" target="_blank" rel="noreferrer">https://code.mpimet.mpg.de/projects/cdo</a></div><div>----</div><div>NCL:</div><div><br></div><div>   f    = addfile("...nc", "r")<br>   x   = f>X<br>   printVarSummary(x)                        ; [time | ...] x [lat | ...] x [lon | ...]</div><div>                                                           ; [time | ...] x [lev |...] x [lat | ...] x [lon | ...]</div><div><br></div><div>   time = x&time                                  ; "seconds/minutes/days since ..."</div><div>   TIME   = <b>cd_calendar</b>(time, 0)             ; type float <br>   year    = toint( TIME(:,0) )               ; toint strips meta data<br>   month = toint( TIME(:,1) )<br>   day     = toint( TIME(:,2) )               ; day of month <br>   hour   = toint( TIME(:,3) )               ; hour of day</div><div><br></div><div>;  x(time,lat,lon)</div><div>                             <br></div><div>   do hh=0,18,6</div><div>        ihr  := ind(hour.eq.hh)      ; indices</div><div>        if (hh.eq.0) then</div><div>           x00 = <b>calculate_monthly_values</b>(x(ihr,:,:), "avg", 0, False)    <br>           x00@long_name = "00Z: Monthly Mean: "+x00@long_name<br>           printVarSummary(x00)<br></div><div><br> </div><div>        elseif (hh.eq. 6) then<br>          x06 = calculate_monthly_values(x(ihr,:,:), "avg", 0, False)<br>          x06@long_name = "06Z: Monthly Mean: "+x06@long_name<br>          printVarSummary(x06)<br><br>      elseif (hh.eq.12) then<br>          x12 = calculate_monthly_values(x(ihr,:,:), "avg", 0, False)<br>          x12@long_name = "12Z: Monthly Mean: "+x12@long_name<br>          printVarSummary(x12)<br><br>      elseif (hh.eq.18) then<br>          x18 = calculate_monthly_values(x(ihr,:,:), "avg", 0, False)<br>          x18@long_name = "18Z: Monthly Mean: "+x18@long_name<br>          printVarSummary(x18)<br>      end if<br>      print("------------------------------------------------")<br>   end do</div><div><br></div><div>Make appropriate changes for X being other dimension ranks<br></div><div><br></div><div>X(time) => x(ihr)</div><div>X(time,pts) => x(ihr,:)</div><div>X(time,lat,lon) => x(ihr,:,:)</div><div>X(time,lev,lat,lon) => x(ihr,:,:,:)</div><div><br></div><div>Good luck<br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Aug 11, 2019 at 1:52 PM Setareh Rahimi via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank" rel="noreferrer">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">Dear all,<div><br></div><div>I have a nc file for my study area, containing hourly values for 6 months at <span style="background-color:rgb(217,210,233)">00, 06, 12, 18.</span> What I am looking for is to have an average for each month over each mentioned time, like :</div><div> April average at 00, 06, 12, 18.</div><div><br></div><div>How can I do this? Is there any function doing this, please?</div><div>Thanks for your advice in advance,</div><div>Best wishes,<br clear="all"><div><br></div>-- <br><div dir="ltr" class="m_-7226652315323787653gmail-m_180681067769693652gmail-m_-4264337966189119083gmail_signature">S.Rahimi<br><br></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank" rel="noreferrer">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank" rel="noreferrer">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>