They are different because the first set is 5 files of 10 years each, and the second set is 50 files of 1 year each. The code is in a loop d = 0,364 to get a single day from all 50 years. So the first set gets 10 records (d::365) per file, while the second gets only 1 record (d) per file and so doesn't need the stride. The first concatenates the five groups of 10 records, as shown in printVarSummary, while the second does not, resulting in only 1 timestep instead of 50. I can add the stride to the second set to see if it then concatenates, but shouldn't it work as written?<br>
<br><div class="gmail_quote">On Sat, Jan 23, 2010 at 7:04 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">concatenation is the default. do you do not need<br>
ListSetType(...,"cat")<br>
<br>
the 1st set of files you did:<br></div><div class="im">
day_data_20c = fall20c[:]->$var$(d::365,:,:)<br>
<br>
<br></div><div class="im">
the 2nd set of files you did<br>
sortobs = fallobs[:]->$var$(d,:,:)<br>
<br>
I think the (d,:,:) should be (d::365,:,:)<br>
<br>
sortobs = fallobs[:]->$var$(d::365,:,:)<br>
<br>
Good luck<br>
<br>
Bridget Thrasher wrote:<br>
</div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div class="h5">
I am finding some peculiar behavior when getting data from a set resulting from a call to addfiles. Both arrays should be the same size (50x180x360), but the second is not concatenating. Am I doing something wrong?<br>
<br>
Code snippets:<br>
<br>
fall20c = addfiles(files20c,"r") ; files20c contains 5 files, each with 10 years of daily data<br>
ListSetType(fall20c,"cat")<br>
<br>
fallobs = addfiles(filesobs,"r") ; filesobs contains 50 files, each with 1 year of daily data<br>
ListSetType(fallobs,"cat")<br>
<br>
day_data_20c = fall20c[:]->$var$(d::365,:,:)<br>
printVarSummary(day_data_20c)<br>
<br>
sortobs = fallobs[:]->$var$(d,:,:)<br>
printVarSummary(sortobs)<br>
<br>
<br>
Output to screen:<br>
<br>
Variable: day_data_20c<br>
Type: float<br>
Total Size: 12960000 bytes<br>
3240000 values<br>
Number of Dimensions: 3<br>
Dimensions and sizes: [time | 50] x [latitude | 180] x [longitude | 360]<br>
Coordinates:<br>
time: [711751.5..729636.5]<br>
latitude: [-89.5..89.5]<br>
longitude: [ 0.5..359.5]<br>
<br>
<br>
Variable: sortobs<br>
Type: float<br>
Total Size: 259200 bytes<br>
64800 values<br>
Number of Dimensions: 2<br>
Dimensions and sizes: [lat | 180] x [lon | 360]<br>
Coordinates:<br>
lat: [-89.5..89.5]<br>
lon: [ 0.5..359.5]<br>
<br>
<br>
<br>
-- <br>
Bridget Thrasher, PhD<br>
Postdoctoral Researcher<br>
Climate Central<br>
</div></div><div class="im"><a href="http://www.climatecentral.org" target="_blank">www.climatecentral.org</a> <<a href="http://www.climatecentral.org" target="_blank">http://www.climatecentral.org</a>><br>
<br>
<br>
<br>
------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</div></blockquote>
</blockquote></div><br><br clear="all"><br>-- <br>Bridget Thrasher, PhD<br>Postdoctoral Researcher<br>Climate Central<br><a href="http://www.climatecentral.org">www.climatecentral.org</a><br><br><br>