function [yeartcprox,tcprox] = MakeProxMCNoPRRecent(randpert,iweight,inorm);

% create surrogate composite paleo Hurricane series from Woodruff data
% by perturbing age model within indicated uncertainties
%
for i=1:2200
    tcprox(i)=0.0;
    wsum(i)=0.0;
end
% 
% determine weighting of regional composites
%
if (iweight==1)
% use return times calculated in Woodruff et al G^3 article to define weights of different records in composite
% take into account variable sensitivity of sites
weight1 = 1/5.37; % Vieques PR 5.37 yrs 
weight2=1/45.66; % Western Lake, FL	45.66 yrs 
weight3=1/10.15; % Alder Island, NY	10.15 yrs
weight4=1/8.52; % Mattapoisett, MA	8.52 yrs
weight5=1/7.84; % Singleton Swash SC 7.84 yrs
end
if (iweight==2)
% use return times calculated in Woodruff et al G^3 article to define weights of different records in composite
% assume uniform sensitivity to cat 3 storms
weight1 = 1/5.37; % Vieques PR 5.37 yrs 
weight2=1/8.78; % Western Lake, FL	8.78 yrs 
weight3=1/25.71; % Alder Island, NY	25.71 yrs
weight4=1/23.81; % Mattapoisett, MA	23.81 yrs
weight5=1/7.84; % Singleton Swash SC 7.84 yrs
end
if (iweight==3)
   weight1=1;
   weight2=1;
   weight3=1;
   weight4=1;
   weight5=1;
end
%
weights=[weight1 weight2 weight3 weight4 weight5];
weightmin=min(weights);
%
weight0 = 0.0*(11.0/32.0)*weight1; % don't use post 1700 Vieques
weight0=weight0/weightmin;
weight1=weight1/weightmin;
weight2=weight2/weightmin;
weight3=weight3/weightmin;
weight4=weight4/weightmin;
weight5=weight5/weightmin;
%
sigma0=1; % perturb age models within 1 sigma age model uncertainties
%
load Data/PR000.dat
nv=size(PR000,1);
yearsv=PR000(nv:-1:1,1);
yearsvl=PR000(nv:-1:1,2);
yearsvh=PR000(nv:-1:1,3);
sigmaplus = yearsvh-yearsv;
sigmaminus = yearsv-yearsvl;
yearsvhnew = yearsv+sigma0*sigmaplus;
yearsvlnew = yearsv-sigma0*sigmaminus;
anorm=1;
for i=1:nv
    if (inorm==1) 
        anorm=nv;
    end
    if (yearsv(i)>0) 
    jj=max(1,round(yearsvlnew(i)+randpert(i,1)*(yearsvhnew(i)-yearsvlnew(i))));
    if (jj<1700)
    weight(1)=weight1/anorm;
    else
    weight(1)=weight0/anorm;
    end
    wsum(jj)=wsum(jj)+weight(1);  
    end
end  
%
% use (Western Lake) Florida record from AD 0-2000
%
weight(2)=weight2;
load Data/FL000.dat
nw=size(FL000,1);
if (inorm==1) 
    anorm=nw;
end
weight(2)=weight(2)/anorm;
yearsw=FL000(nw:-1:1,1);
yearswl=FL000(nw:-1:1,2);
yearswh=FL000(nw:-1:1,3);
sigmaplus = yearswh-yearsw;
sigmaminus = yearsw-yearswl;
yearswhnew = yearsw+sigma0*sigmaplus;
yearswlnew = yearsw-sigma0*sigmaminus;
% keep all events after year 0
for i=1:nw
    if (yearsw(i)>0) 
    jj=max(1,round(yearswlnew(i)+randpert(i,2)*(yearswhnew(i)-yearswlnew(i))));
    wsum(jj)=wsum(jj)+weight(2);
    end
end 
%
% use New Jersey/New York (3 record) composite AD 0-date
% (Brigantine NJ AD 600-date; Whale Lake NJ 1300-date; 
% + Alder Lake NY AD 0-date)
%
weight(3) = weight3;
load Data/NYNJ000.dat
nb=size(NYNJ000,1);
if (inorm==1) 
    anorm=nb;
end
weight(3)=weight(3)/anorm;
yearsb=NYNJ000(nb:-1:1,1);
yearsbl=NYNJ000(nb:-1:1,2);
yearsbh=NYNJ000(nb:-1:1,3);
sigmaplus = yearsbh-yearsb;
sigmaminus = yearsb-yearsbl;
yearsbhnew = yearsb+sigma0*sigmaplus;
yearsblnew = yearsb-sigma0*sigmaminus;
% keep all events after year 0
for i=1:nb
    if (yearsb(i)>0)
    jj=max(1,round(yearsblnew(i)+randpert(i,3)*(yearsbhnew(i)-yearsblnew(i))));
    wsum(jj)=wsum(jj)+weight(3);
    end
end 
%
%
% use Rhode Island/MA (2 record) composite AD 0-date
% (Succotash Marsh RI AD 1300-date; Mattapoisett Marsh MA O-date; 
%
weight(4) = weight4;
load Data/RIMA000.dat
nn=size(RIMA000,1);
if (inorm==1) 
    anorm=nn;
end
weight(4)=weight(4)/anorm;
yearsn=RIMA000(nn:-1:1,1);
yearsnl=RIMA000(nn:-1:1,2);
yearsnh=RIMA000(nn:-1:1,3);
sigmaplus = yearsnh-yearsn;
sigmaminus = yearsn-yearsnl;
yearsnhnew = yearsn+sigma0*sigmaplus;
yearsnlnew = yearsn-sigma0*sigmaminus;
% keep all events after AD 0
for i=1:nn
    if (yearsn(i)>0)
    jj=max(1,round(yearsnlnew(i)+randpert(i,4)*(yearsnhnew(i)-yearsnlnew(i))));
    wsum(jj)=wsum(jj)+weight(4);
    end
end 
%
% use South Carolina site (1 record) composite AD 0-date
% (Singleton Swash SC AD O-date)
%
weight(5) = weight5;
load Data/SC000.dat
nc=size(SC000,1);
if (inorm==1) 
    anorm=nc;
end
weight(5)=weight(5)/anorm;
yearsc=SC000(nc:-1:1,1);
yearscl=SC000(nc:-1:1,2);
yearsch=SC000(nc:-1:1,3);
sigmaplus = yearsch-yearsc;
sigmaminus = yearsc-yearscl;
yearschnew = yearsc+sigma0*sigmaplus;
yearsclnew = yearsc-sigma0*sigmaminus;
% keep all events after AD 0
for i=1:nc
    if (yearsc(i)>0)
    jj=max(1,round(yearsclnew(i)+randpert(i,5)*(yearschnew(i)-yearsclnew(i))));
    wsum(jj)=wsum(jj)+weight(5);
    end
end 
%
% adjust for fewer series (i.e. no Vieques) after AD 1700  
%
for i=1700:2200
   wsum(i)=wsum(i)*(weight1+weight2+weight3+weight4+weight5)/(weight2+weight3+weight4+weight5);
end
yeartcprox=[500:1991];
%
% smooth series on 40 year timescale to insure roughly uniform resolution
wsum=wsum';
tcfull=lowpassmin(wsum,0.025);
tcprox=tcfull(500:1991);
end



