
addpath('../')
addpath(genpath('boundedline-pkg/'))
load('./colormaps/distinct_colors_20.mat')

%-----------------------------------------

% ALL PROXY TYPES, LMR DATABASE
%efl='cesm_lme010_r12000_p998_state_t2qjozn_avg_AprMar_prxydtst_5_prxtp_tca_2896_swtchbld100_05-Jan-2018_14:33:07.mat';
%efl='cesm_lme010_r12000_p998_state_t2qjozn_avg_DecFeb_prxydtst_5_prxtp_tca_2896_swtchbld100_05-Jan-2018_14:40:23.mat';
efl='cesm_lme010_r12000_p998_state_t2qjozn_avg_JunAug_prxydtst_5_prxtp_tca_2896_swtchbld100_05-Jan-2018_14:49:45.mat';

disp(' ')
disp(['Analyzing file: ' efl])


%pth='/d2/nsteiger/output-da/hydroclimate/'; load([pth,efl])
load(['./output-recon/',efl],'proxy')

% new path for figures
%prntpth='./figs/';


%------------------------------------------
% HOW WERE THE TREE RINGS MODELED?
j=0;k=0;
for i=1:length(proxy)

   if strcmp(proxy{i}.ptype,'Tree Rings')

      if strcmp(proxy{i}.psm_mthd,'p')
         j=j+1;
      elseif strcmp(proxy{i}.psm_mthd,'t')
         k=k+1;
      end

   end

end

disp(['PDSI PSM = ' num2str(j)])
disp(['Temp PSM = ' num2str(k)])


