<div dir="ltr"><div><b>This is pure dog-work programming. </b><br></div><div><br></div><div>NTIM - maximum number of hourly time steps**IF** all files were available. <br></div><div><br></div><div>So if you had January and February (leap year) data files.</div><div><br></div><div>NTIM = 31*24 + 29*24 = 1440<br></div><div><br></div><div>*IF* the yyyy_mm_ddhh files are (say) (lat,lon) => (nlat,mlon), create a 'super array'<br></div><div><br></div><div> NTIM = 1440<br></div><div> DATA = new ( (/NTIM,nlat,mlon/), typeof(data), getVarFillValue(data)) ; super array<br></div><div> DATA!0 = "YMDH"</div><div><br></div><div> fili = systemfunc(...)</div><div> nfili = dimsizes(fili)</div><div><br></div><div> do nf=0,nfili-1<br></div><div> f = addfile(...fili(nf),"r")</div><div> x = f->X</div><div> ; nt is the appropriate index into super array<br></div><div> DATA(nt,:,:) = (/ x /)</div><div> end do</div><div><br></div><div> DATA = linmsg_n(DATA,0,0)</div> <br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 1, 2019 at 9:49 AM Ipsita Putatunda <<a href="mailto:ipsita.putatunda@gmail.com" target="_blank">ipsita.putatunda@gmail.com</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">Hi Dennis,<div>The link opened after removing the "test". But my issue is reading the hourly files, which are missing; not the variables.</div><div>Suppose I have the files as below-</div><div><a href="http://1980_01_0100.nc" target="_blank">1980_01_0100.nc</a></div><div><a href="http://1980_01_0102.nc" target="_blank">1980_01_0102.nc</a><br></div><div><a href="http://1980_01_0103.nc" target="_blank">1980_01_0103.nc</a><br></div><div><a href="http://1980_01_0105.nc" target="_blank">1980_01_0105.nc</a><br></div><div><a href="http://1980_01_0107.nc" target="_blank">1980_01_0107.nc</a> etc......<br></div><div>In that case how the interpolation will work to read the files and making 3 hourly files? Hope I make you understand my problem.</div><div><br></div><div>Thanks,</div><div>Ipsita</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 1, 2019 at 8:21 PM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@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>Hi Ipsita,</div><div><br></div><div>Try removing the "test." from the first part of that URL -- that's an internal testing site for preparing documentation.</div><div><br></div><div>RIck<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 1, 2019 at 8:15 AM Ipsita Putatunda 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">Hi Dennis,<div>I tried several times but the link you have shared is not opening. I have hourly data files for 30 years. In some cases there are 3 successive hourly files, and in some cases it is missing. In that case how can I put it in a loop for 30 years?</div><div><br></div><div>Thanks,</div><div>Ipsita</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 1, 2019 at 5:24 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@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>There is no magic ...</div><div><br></div><div>If you do not have 3 successive hourly files then you should perform linear interpolation"</div><div><a href="http://test.www.ncl.ucar.edu/Document/Functions/Built-in/linmsg_n.shtml" target="_blank"><b>http://test.www.ncl.ucar.edu/Document/Functions/Built-in/linmsg_n.shtml</b></a><br></div><div><br></div><div>If in some cases you have<b> 0 </b>files then the derived values are just place holders.<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Aug 1, 2019 at 12:57 AM Ipsita Putatunda 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">Dear NCl users,<div>I have hourly files, which I want to convert in 3 hourly files by averaging. But for everyday the hourly files are not continuous, and there are some missing hourly files. That is why the number of files present in each three hours is varying from 0 to 3. Hence using systemfunc function is giving error, where there are missing files.</div><div>Can anyone please suggest me how can I get a 3 hourly average in this case.</div><div><br></div><div>Thanks in advance,</div><div>Ipsita</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>
</blockquote></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>
</blockquote></div>
</blockquote></div>