%
% Constanza Amaya - Universidad de Concepción  01/2021
%

close all
clear all


%% VEO LOS DATOS DE PATHFINDER ANUAL.

ncdisp('climato_pathfinder.nc')


%% Abro el archivo y separo por cada variable.

Longitud = ncread('climato_pathfinder.nc','X');
Latitud = ncread('climato_pathfinder.nc','Y');
Meses  = ncread('climato_pathfinder.nc','T');
SST = ncread('climato_pathfinder.nc','SST');

%% Selecciono mi area de estudio (Juan Fernandez)
Lon=Longitud(1104:1195);
Lat=Latitud(603:683);
SST_JF=SST(1104:1195,603:683,:); % Temperatura del mar en JF

%Calculo promedio anual de SST
SST_JF_ANUAL=nanmean(SST_JF,3);

%Calculo promedios estacionales de SST
SST_JF_otono=nanmean(SST_JF(:,:,3:5),3);%Marzo, Abril y Mayo 
SST_JF_invierno=nanmean(SST_JF(:,:,6:8),3);%Junio, Julio y Agosto
SST_JF_primavera=nanmean(SST_JF(:,:,9:11),3);%Septiembre, Octubre y Noviembre
SST_JF_V=cat(3,SST_JF(:,:,1:2),SST_JF(:,:,12));% Diciembre, Enero, Febrero
SST_JF_verano=nanmean(SST_JF_V,3);


%% PLOTEO 
load coastline_l.mat


figure(1)

pcolor(Lon,Lat,SST_JF_ANUAL')
xlabel('Longitud')
ylabel('Latitud')
title('TSM PATHFINDER: promedio anual ')
s.FaceColor = 'interp';
colorbar



subplot(2,2,1)
pcolor(Lon,Lat,SST_JF_otono')
xlabel('Longitud')
ylabel('Latitud')
title('TSM PATHFINDER: Otoño ')
s.FaceColor = 'interp';
colorbar


subplot(2,2,2)
pcolor(Lon,Lat,SST_JF_invierno')
xlabel('Longitud')
ylabel('Latitud')
title('TSM PATHFINDER: Invierno ')
s.FaceColor = 'interp';
colorbar


subplot(2,2,3)
pcolor(Lon,Lat,SST_JF_primavera')
xlabel('Longitud')
ylabel('Latitud')
title('TSM PATHFINDER: Primavera ')
s.FaceColor = 'interp';
colorbar

subplot(2,2,4)
pcolor(Lon,Lat,SST_JF_verano')
xlabel('Longitud')
ylabel('Latitud')
title('TSM PATHFINDER: Verano ')
s.FaceColor = 'interp';
colorbar %jet

%% CARGO LOS PROMEDIOS ESTACIONALES (CROCO)

Lon1 = ncread('croco_OTONO.nc','lon_rho');
Lon1=Lon1(:,1);
Lat1 = ncread('croco_OTONO.nc','lat_rho');
Lat1=Lat1(1,:)';


Temp_otono = ncread('croco_OTONO.nc','temp');
Temp_invierno = ncread('croco_INVIERNO.nc','temp');
Temp_primavera = ncread('croco_PRIMAVERA.nc','temp');
Temp_verano = ncread('croco_VERANO.nc','temp');

%% CARGO LOS PROMEDIO ANUAL (CROCO)
Temp_anual = ncread('croco_PromedioAnual.nc','temp');

%% Selecciono la capa superficial
SST_otono=Temp_otono(:,:,32);
SST_invierno=Temp_invierno(:,:,32);
SST_primavera=Temp_primavera(:,:,32);
SST_verano=Temp_verano(:,:,32);
SST_anual=Temp_anual(:,:,32);

%% INTERPOLO PARA DEJAR LOS DATOS DE PATHFINDER IGUAL A LOS DE CROCO
SST_JF_otonoo=[];
SST_JF_inviernoo=[];
SST_JF_primaveraa=[];
SST_JF_veranoo=[];
SST_JF_anuall=[];

c=0;
for i=1:81
    c=c+1;
SST_JF_otonoo(:,c) = interp1(Lon,SST_JF_otono(:,i),Lon1,'spline');
SST_JF_inviernoo(:,c) = interp1(Lon,SST_JF_invierno(:,i),Lon1,'spline');
SST_JF_primaveraa(:,c) = interp1(Lon,SST_JF_primavera(:,i),Lon1,'spline');
SST_JF_veranoo(:,c) = interp1(Lon,SST_JF_verano(:,i),Lon1,'spline');
SST_JF_anuall(:,c) = interp1(Lon,SST_JF_ANUAL(:,i),Lon1,'spline');
end



SST_JF_otono2=[];
SST_JF_invierno2=[];
SST_JF_primavera2=[];
SST_JF_verano2=[];
SST_JF_anual2=[];

c=0;
for i=1:81
    c=c+1;
SST_JF_otono2(c,:) = interp1(Lat,SST_JF_otonoo(i,:),Lat1,'spline');
SST_JF_invierno2(c,:) = interp1(Lat,SST_JF_inviernoo(i,:),Lat1,'spline');
SST_JF_primavera2(c,:) = interp1(Lat,SST_JF_primaveraa(i,:),Lat1,'spline');
SST_JF_verano2(c,:) = interp1(Lat,SST_JF_veranoo(i,:),Lat1,'spline');
SST_JF_anual2(c,:) = interp1(Lat,SST_JF_anuall(i,:),Lat1,'spline');

end

%% DIFERENCIA ENTRE PATHFINDER Y RESULTADOS DE CROCO

sst_dif_otono=SST_JF_otono2-SST_otono;
sst_dif_invierno=SST_JF_invierno2-SST_invierno;
sst_dif_primavera=SST_JF_primavera2-SST_primavera;
sst_dif_verano=SST_JF_verano2-SST_verano;
sst_dif_anual=SST_JF_anual2-SST_anual;


%% PLOTEO
load coastline_l_mask.mat


pcolor(Lon1,Lat1,sst_dif_otono')
geoshow(lat,lon,'color','k')
xlabel('Longitud')
ylabel('Latitud')
title('Diferencia de TSM entre PATHFINDER y CROCO : Otoño ')
s.FaceColor = 'interp';
colormap(whitejet);
c=colorbar;
caxis([-1.5 1.5]); 

pcolor(Lon1,Lat1,sst_dif_invierno')
geoshow(lat,lon,'color','k')
xlabel('Longitud')
ylabel('Latitud')
title('Diferencia de TSM entre PATHFINDER y CROCO : Invierno ')
s.FaceColor = 'interp';
colormap(whitejet);
c=colorbar %jet
caxis([-1.5 1.5]);


pcolor(Lon1,Lat1,sst_dif_primavera')
geoshow(lat,lon,'color','k')
xlabel('Longitud')
ylabel('Latitud')
title('Diferencia de TSM entre PATHFINDER y CROCO : Primavera ')
s.FaceColor = 'interp';
colormap(whitejet);
c=colorbar %jet
caxis([-1.5 1.5]);


pcolor(Lon1,Lat1,sst_dif_verano')
geoshow(lat,lon,'color','k')
xlabel('Longitud')
ylabel('Latitud')
title('Diferencia de TSM entre PATHFINDER y CROCO : Verano ')
s.FaceColor = 'interp';
colormap(whitejet);
c=colorbar;
caxis([-1.5 1.5]); 

figure(3)
pcolor(Lon1,Lat1,sst_dif_anual')
geoshow(lat,lon,'color','k')
xlabel('Longitud')
ylabel('Latitud')
title('Diferencia de TSM entre PATHFINDER y CROCO : Anual ')
s.FaceColor = 'interp';
colormap(whitejet);
c=colorbar;
caxis([-1.2 1.2]);

%% Histograma de las diferencias calculadas

figure(4)
subplot(3,2,1)
hist(reshape(sst_dif_otono,[81*86,1]))
ylabel('Frecuencia')
xlabel('[°C]')
xlim ([-1.5 1.5])
subplot(3,2,2)
hist(reshape(sst_dif_invierno,[81*86,1]))
xlim ([-1.5 1.5])
ylabel('Frecuencia')
xlabel('[°C]')
subplot(3,2,3)
hist(reshape(sst_dif_primavera,[81*86,1]))
xlim ([-1.5 1.5])
ylabel('Frecuencia')
xlabel('[°C]')
subplot(3,2,4)
hist(reshape(sst_dif_verano,[81*86,1]))
xlim ([-1.5 1.5])
ylabel('Frecuencia')
xlabel('[°C]')
subplot(3,2,5)
hist(reshape(sst_dif_anual,[81*86,1]))
xlim ([-1.5 1.5])
ylabel('Frecuencia')
xlabel('[°C]')

