-
Notifications
You must be signed in to change notification settings - Fork 0
/
masas_agua.m
49 lines (39 loc) · 1.46 KB
/
masas_agua.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
%Pasos
%% PASO 1
%Primero descargar los datos de https://resources.marine.copernicus.eu/products
%Seleccionar los modelos resolucion de 8km o 1/12º :
%GLOBAL_REANALYSIS_PHY_001_030/ periodo: 2014-2019 &&
%GLOBAL_ANALYSIS_FORECAST_PHY_001_024/ periodo: 2019-2021
%Descargar los datos diarios de Sal y Temp
%% Paso 2
%Aqui vamos a seleccionar solo los dias 15 de los archivos nc que hemos
%descargado de ambos modelos y los procesamos. Fin. xd
path01='D:\CIO\seminario_CIO\Sesion_2';
hdir=dir(fullfile(path01,'20*sal.nc'));
load('WM_peru');
%aviobj = QTWriter('Masas_awa.mov','FrameRate',1);%aviobj.Quality = 100;
for icmems=1:1:size(hdir,1)
fn=hdir(icmems).name;
time=double(ncread(fn,'time'))./24;
time=time+datenum(1950,1,1,0,0,0);
[yr,mo,da,hr,mi,se]=datevec(time);
lon=double(ncread(fn,'longitude'));
lat=double(ncread(fn,'latitude'));
indx_day=find(da==15 | da==30);%2 quincenas
% indx_day=find(da==15);
sal=double(ncread(fn,'so')); sal=permute(sal,[1 2 4 3]);
sal_sel=sal(:,:,indx_day); %selecciona el dia de la salinidad
t = datetime(yr,mo,da);
fecha=t(indx_day);
for ii=1:1:size(indx_day)
pcolor(lon,lat,sal_sel(:,:,ii)'); shading flat; colormap(peru_wm);
colorbar; caxis([33 35.5]);
title(datestr(fecha(ii)));
% M1=getframe(gcf);
% writeMovie(aviobj, M1);
pause(1)
clf
end
end
% close(aviobj);
%% Superficie