-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathplot_climato_puertos.m
More file actions
36 lines (33 loc) · 978 Bytes
/
plot_climato_puertos.m
File metadata and controls
36 lines (33 loc) · 978 Bytes
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
clear all; close all; clc;
%%
cd D:\trabajo\IGP\CLIM_PEZ\Dra_Jaramillos
%% climatologia import
[data,raw,j]=xlsread('promedios_caudales_correcto_requena_yurimaguas.xlsx',...
'Climatologia_mensual');
%%
mes=data(:,1);
chazuta_caudal=data(:,3);
sregis_caudal=data(:,4);
req_caudal=data(:,6);
tamshi_caudal=data(:,7);
%% plots
puertos={'Tamshiyacu','San Regis','Requena','Chazuta'};
labels={'Ene','Feb','Mar','Abr','May','Jun','Jul','Ago','Set','Oct','Nov','Dic'};
figure
% P=get(gcf,'position');
% P(3)=P(3)*4;
% P(4)=P(4)*1;
% set(gcf,'position',P);
% set(gcf,'PaperPositionMode','auto');
plot(mes,tamshi_caudal,'.--');
hold on
plot(mes,sregis_caudal,'^--');
hold on
plot(mes,req_caudal,'d--');
hold on
plot(mes,chazuta_caudal,'.:');
grid on
set(gca,'xtick',[1:1:12],'xticklabel',labels,'xlim',[1 12],'Fontsize',8);
ylabel('Caudal [m^3/s]'); xlabel('Meses');
title('Climatologia mensual Caudal [m^3/s]','Fontsize',12);
legend(puertos);