
% Load all pages2k proxies to get names and right lat lons
load(['../../proxy-data/PAGES2k_v2/v2_0_0/proxy_ama_2.0.0_highres.mat'])
pn=proxy_names;
lt=p_lat;ln=p_lon;
ar=archive;
clearvars -except ar pn lt ln


fnm='AR-Noise_100_proxy_ama_2.0.0_calib-selection_1881_1916_1995_0.67_infilled_DINEOF_PAGES-crit-regionalFDR.mat';
load p1.mat
[~,ip]=intersect(pn,proxy_names);
p_lat=lt(ip);
p_lon=ln(ip);
archive=ar(ip);
save(['./' fnm],'archive','p_lat','p_lon','proxy_ama','proxy_names','year')


clearvars -except ar pn lt ln
fnm='AR-Noise_100_proxy_ama_2.0.0_highres_regional-FDR_screen_infilled.mat';
load p2.mat
[~,ip]=intersect(pn,proxy_names);
p_lat=lt(ip);
p_lon=ln(ip);
archive=ar(ip);
save(['./' fnm],'archive','p_lat','p_lon','proxy_ama','proxy_names','year')


clearvars -except ar pn lt ln
fnm='AR-Noise_100_proxy_ama_2.0.0_highres_regional_screen_infilled.mat';
load p3.mat
[~,ip]=intersect(pn,proxy_names);
p_lat=lt(ip);
p_lon=ln(ip);
archive=ar(ip);
save(['./' fnm],'archive','p_lat','p_lon','proxy_ama','proxy_names','year')





