an class="comment">% Construct the data filename</span>
0094     dfile = [DataFile,num2str(years(p))];
0095     
0096     <span class="comment">% Open the surge model data file</span>
0097     fid = fopen(dfile,<span class="string">'r'</span>,<span class="string">'b'</span>);
0098     
0099     PASTIT = false;
0100     
0101     <span class="keyword">while</span> PASTIT == false
0102         
0103         <span class="comment">% Read the data file</span>
0104         dump=fread(fid,1,<span class="string">'*