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

load Theta_Reg_2.mat
% Retorna
%   fec        312x3                  7488  double               
%   latr        80x1                   640  double               
%   lonr        40x1                   320  double               
%   theta        4-D             399360000  double               
%   z           50x1                   400  double               

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

% % Se deben remover valore de datos no válidos y convertir a NaN
 ind = (theta > 100 | theta < -10);
 theta(ind) = NaN;
 
 % Calcula ciclo anual y anomalías
 SAtheta = NaN(40,80,50,12);
 for i = 1:12
    SAtheta(:,:,:,i) = mean(theta(:,:,:,i:12:end),4);
 end
 % genera arreglo del ciclo anual con tamaño del arreglo mayor 26 años
 SAthetaL = cat(4,SAtheta,SAtheta,SAtheta,SAtheta, SAtheta, SAtheta, ...
                  SAtheta,SAtheta,SAtheta,SAtheta, SAtheta, SAtheta, ...
                  SAtheta,SAtheta,SAtheta,SAtheta, SAtheta, SAtheta, ...
                  SAtheta,SAtheta,SAtheta,SAtheta, SAtheta, SAtheta, ...
                  SAtheta,SAtheta);
 % Anomalías
 Atheta = theta - SAthetaL;
 
 
% % por ejemplo grafica un grafico temporal en una posición (ditintas prof)
lat0 = 40; % corresponde a -30.1250
lon0 = 20; % -75.1250
figure
pcolor(datenum(fec), -z(1:30), squeeze(Atheta(lon0,lat0,1:30,:)))
shading interp; colormap(jet); datetick

% % Ciclo anual en posición lon0 lat0
% pcolor((1:12), -z(1:15), squeeze(SAtheta(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(Atheta(:,lat0,17,:))')
 shading interp; colormap(jet); datetick('y'); colorbar




