% 
% Analiza datos de ECCO en region 2
% Orpa Nov 2020

load Salt_Reg_2.mat
% Retorna
%   fec        312x3                  7488  double               
%   latr        80x1                   640  double               
%   lonr        44x1                   320  double               
%   Salt        4-D             399360000  double  [44    80    50   312]
%   z           50x1                   400  double               

% El arreglo Salt contiene 4-D:   [lon, lat, z, tiempo]

% % Se deben remover valore de datos no válidos y convertir a NaN
 ind = (salt > 100 | salt < -10);
 salt(ind) = NaN;
 
 % Calcula ciclo anual y anomalías
 SAsalt = NaN(44,80,50,12);
 for i = 1:12
    SAsalt(:,:,:,i) = mean(salt(:,:,:,i:12:end),4);
 end
 % genera arreglo del ciclo anual con tamaño del arreglo mayor 26 años
 SAsaltL = cat(4,SAsalt,SAsalt,SAsalt,SAsalt, SAsalt, SAsalt, ...
                  SAsalt,SAsalt,SAsalt,SAsalt, SAsalt, SAsalt, ...
                  SAsalt,SAsalt,SAsalt,SAsalt, SAsalt, SAsalt, ...
                  SAsalt,SAsalt,SAsalt,SAsalt, SAsalt, SAsalt, ...
                  SAsalt,SAsalt);
 % Anomalías
 Asalt = salt - SAsaltL;
 
 
% % por ejemplo grafica un grafico temporal en una posición (ditintas prof)
lat0 = 15;  % 15=>36.4°S    40 => -30.1250
lon0 = 20;  % 25=>73.875    20 =>-75.1250
figure
pcolor(datenum(fec), -z(1:30), squeeze(salt(lon0,lat0,1:30,:)))
hold on; contour(datenum(fec), -z(1:30), squeeze(salt(lon0,lat0,1:30,:)),[34.5 34.5],'k')
shading interp; colormap(jet); datetick; set(gca,'tickdir', 'out')
colorbar

figure
pcolor(datenum(fec), -z(1:30), squeeze(Asalt(lon0,lat0,1:30,:)))
shading interp; colormap(jet); datetick; set(gca,'tickdir', 'out')
colorbar


% % Ciclo anual en posición lon0 lat0
% pcolor((1:12), -z(1:15), squeeze(SAsalt(lon0,lat0,1:15,:)))
% shading interp; colormap(jet); 

% Hovermole lon tiempo a un z dado
 lat0 = 40; % 76;  % 21°S
 figure
 pcolor(lonr, datenum(fec), squeeze(Asalt(:,lat0,17,:))')
 shading interp; colormap(jet); datetick('y'); colorbar




