%####################################################
%#########                                    #######
%######### put infilled into allproxy1227       #######
%######### all back to 1998                   #######
%#########                                    #######
%####################################################

clear all
close all


y=load('d:\scozztemann\trimdata\triminfilled1227');
means=load('d:\scozztemann\trimdata\means');
stdevs=load('d:\scozztemann\trimdata\stdevs');

    [ny mpc]=size(y)
    y=y.*repmat(stdevs,ny,1);
    y=y+repmat(means,ny,1);


dataf = 'd:\scozztemann\trimdata\1227proxynamenolu.txt';

filee = textread(dataf,'%s','delimiter','\n','whitespace','');
[m1,n1] = size(filee)

for i =1:m1 
   cc=lower(filee{i}); dd='d:\scozztemann\dodist\allproxy1\'; aa='.ppd';
   ccc=strcat(dd,cc,aa);
   cf = load(strtok(ccc));
   m=length(cf);
   endyr=cf(m,1);
   if endyr>=1998
       cff=cf;
   else
       mm=m+1998-endyr;
       cff=repmat(NaN,mm,3);
       cff(1:mm-199,:)=cf(1:mm-199,:);
       cff(mm-198:mm,1)=(1800:1:1998)';
       cff(mm-198:mm,2)=y(:,i);
       cff(1:m,3)=cf(1:m,3);
       cff(m+1:mm,3)=cf(m,3);
   end
   save(strcat('d:\scozztemann\trimdata\allproxy1227\',cc,'.ppd'),'cff','-ASCII') 
%   save(strcat('d:\scozztemann\trimdata\test\',cc,'.ppd'),'cff','-ASCII') 
   clear cc ccc cf cff m mm endyr
end

%----------
%####################################################
%####################################################