    
0148     <span class="comment">% Since there are many missing values before 1989, we will only</span>
0149     <span class="comment">% consider the period from October 1989 to September 2009.</span>
0150     t = datetime(CAMELS_GB_data.Q{i}(:,1),<span class="string">'ConvertFrom'</span>,<span class="string">'datenum'</span>);
0151     Q = CAMELS_GB_data.Q{i}(:,2);
0152     P = CAMELS_GB_data.P{i}(:,2);
0153     PET = CAMELS_GB_data.PET{i}(:,2);
0154     
0155     indices = 1:length(t);
0156     start_i