<div dir="ltr"><div class="gmail_extra"><div class="gmail_default" style="font-family:verdana,sans-serif">The easiest way to do this is using the time functions and indexing. </div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">ymdhms = cd_calendar(x, 0) ; </div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">dec_ind = ind(ymdhms(:,1) .eq. 12) ;; find locations where month .eq. 12. </div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">x_dec = x(dec_ind,:,:)</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Adding another 2 dimensions, doesn&#39;t seem necessary. You can stack the ind() function to be as specific or as broad as you like very easily. Will be more flexible than multiple time dimension. </div><br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><div class="gmail_default" style="font-family:verdana,sans-serif">​Alan. ​</div><br></div><div class="gmail_extra"><br></div><div class="gmail_extra"><br><div class="gmail_quote">On 10 June 2015 at 11:58, Zachary Suriano <span dir="ltr">&lt;<a href="mailto:zsuriano@udel.edu" target="_blank">zsuriano@udel.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><br></div>Ultimately I would like to pull out data from each month individually with a statement such as: dec=x( : , 12 , : , : , : ). With leap days present within the dataset, I cannot figure out a way to achieve this without having year, month, day being their own dimensions.<br></blockquote></div><br><br></div></div>