%% epoch 1000
n_c=[15	20	25	30	35	40];
rmse1=[0.0140	0.0049
0.0120	0.0036
0.0081	0.0017
0.0224	0.0070
0.0163	0.0050
0.1250	0.0371
];
rmse2=[0.0202	0.0059
0.0105	0.0031
0.0138	0.0056
0.0169	0.0033
0.0643	0.0148
0.0391	0.0121
];
rmse3=[0.0261	0.0086
0.0303	0.0087
0.0066	0.0023
0.0308	0.0105
0.0652	0.0163
0.0566	0.0160
];
rmse4=[0.0042	0.0034
0.0250	0.0075
0.0277	0.0070
0.0085	0.0053
0.0129	0.0059
0.0149	0.0043
];

figure()
subplot(1,2,1)
plot(n_c,rmse1(:,1),'sk','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','k')
hold on
plot(n_c,rmse2(:,1),'sb','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','b')
hold on
plot(n_c,rmse3(:,1),'sr','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','r')
hold on
plot(n_c,rmse4(:,1),'sg','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','g')
xlabel('numer of channels in 2$^{nd}$ c-layer','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_h$ (m)','Interpreter','latex','FontSize',25);
ylim([0 0.1])
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Depth','Interpreter','latex','FontSize',25)

subplot(1,2,2)
plot(n_c,rmse1(:,2),'ok','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','k')
hold on
plot(n_c,rmse2(:,2),'ob','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','b')
hold on
plot(n_c,rmse3(:,2),'or','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','r')
hold on
plot(n_c,rmse4(:,2),'og','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','g')
xlabel('numer of channels in 2$^{nd}$ c-layer','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_{hu}$ (m$^2$/s)','Interpreter','latex','FontSize',25);
ylim([0 0.02])
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)

figure()
subplot(1,2,1)
bar(n_c,[rmse1(:,1)';rmse2(:,1)';rmse3(:,1)';rmse4(:,1)'],1)
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
xlabel('numer of channels','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_h$ (m)','Interpreter','latex','FontSize',25);
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Depth','Interpreter','latex','FontSize',25)

subplot(1,2,2)
bar(n_c,[rmse1(:,2)';rmse2(:,2)';rmse3(:,2)';rmse4(:,2)'],1)
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
xlabel('numer of channels','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_{hu}$ (m$^2$/s)','Interpreter','latex','FontSize',25);
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)


%% epoch 2000
n_c=[15	20	25	30	35	40];
rmse1=[0.0014	0.0004
0.0014	0.0004
0.0030	0.0009
0.0016	0.0005
0.0031	0.0009
0.0013	0.0004
];
rmse2=[0.0029	0.0008
0.0030	0.0009
0.0052	0.0015
0.0037	0.0011
0.0015	0.0004
0.0014	0.0004
];
rmse3=[0.0060	0.0017
0.0216	0.0063
0.1124	0.0326
0.0166	0.0043
0.0018	0.0005
0.0789	0.0122
];
rmse4=[0.0043	0.0013
0.0030	0.0009
0.0619	0.0179
0.0065	0.0019
0.0149	0.0043
0.0037	0.0011
];

figure()
subplot(1,2,1)
plot(n_c,rmse1(:,1),'sk','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','k')
hold on
plot(n_c,rmse2(:,1),'sb','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','b')
hold on
plot(n_c,rmse3(:,1),'sr','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','r')
hold on
plot(n_c,rmse4(:,1),'sg','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','g')
xlabel('numer of channels in 2$^{nd}$ c-layer','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_h$ (m)','Interpreter','latex','FontSize',25);
ylim([0 0.02])
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Depth','Interpreter','latex','FontSize',25)

subplot(1,2,2)
plot(n_c,rmse1(:,2),'ok','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','k')
hold on
plot(n_c,rmse2(:,2),'ob','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','b')
hold on
plot(n_c,rmse3(:,2),'or','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','r')
hold on
plot(n_c,rmse4(:,2),'og','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','g')
xlabel('numer of channels in 2$^{nd}$ c-layer','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_{hu}$ (m$^2$/s)','Interpreter','latex','FontSize',25);
ylim([0 0.01])
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)

figure()
subplot(1,2,1)
bar(n_c,[rmse1(:,1)';rmse2(:,1)';rmse3(:,1)';rmse4(:,1)'],1)
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
xlabel('numer of channels','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_h$ (m)','Interpreter','latex','FontSize',25);
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Depth','Interpreter','latex','FontSize',25)

subplot(1,2,2)
bar(n_c,[rmse1(:,2)';rmse2(:,2)';rmse3(:,2)';rmse4(:,2)'],1)
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
xlabel('numer of channels','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_{hu}$ (m$^2$/s)','Interpreter','latex','FontSize',25);
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)


%% epoch 3000
n_c=[15	20	25	30	35	40];
rmse1=[0.0033	0.0009
0.0042	0.0012
0.0022	0.0006
0.0034	0.0010
0.0019	0.0005
0.0014	0.0004
];
rmse2=[0.0065	0.0019
0.0057	0.0016
0.2001	0.0580
0.1600	0.0464
0.0035	0.0010
0.2008	0.0582
];
rmse3=[0.0044	0.0013
0.0029	0.0009
0.0119	0.0034
0.0036	0.0010
0.1502	0.0402
0.0011	0.0003
];
rmse4=[0.2001	0.0580
0.0055	0.0016
0.0145	0.0042
0.2001	0.0580
0.0022	0.0006
0.0019	0.0006
];

figure()
subplot(1,2,1)
plot(n_c,rmse1(:,1),'sk','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','k')
hold on
plot(n_c,rmse2(:,1),'sb','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','b')
hold on
plot(n_c,rmse3(:,1),'sr','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','r')
hold on
plot(n_c,rmse4(:,1),'sg','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','g')
xlabel('numer of channels in 2$^{nd}$ c-layer','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_h$ (m)','Interpreter','latex','FontSize',25);
ylim([0 0.02])
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Depth','Interpreter','latex','FontSize',25)

subplot(1,2,2)
plot(n_c,rmse1(:,2),'ok','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','k')
hold on
plot(n_c,rmse2(:,2),'ob','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','b')
hold on
plot(n_c,rmse3(:,2),'or','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','r')
hold on
plot(n_c,rmse4(:,2),'og','LineWidth',1,'MarkerSize',12,'MarkerFaceColor','g')
xlabel('numer of channels in 2$^{nd}$ c-layer','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_{hu}$ (m$^2$/s)','Interpreter','latex','FontSize',25);
ylim([0 0.005])
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)

figure()
subplot(1,2,1)
bar(n_c,[rmse1(:,1)';rmse2(:,1)';rmse3(:,1)';rmse4(:,1)'],1)
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
xlabel('numer of channels','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_h$ (m)','Interpreter','latex','FontSize',25);
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Depth','Interpreter','latex','FontSize',25)

subplot(1,2,2)
bar(n_c,[rmse1(:,2)';rmse2(:,2)';rmse3(:,2)';rmse4(:,2)'],1)
legend('f-neuron-25','f-neuron-50','f-neuron-75','f-neuron-100')
xlabel('numer of channels','Interpreter','latex','FontSize',25);
ylabel('$\mathcal{R}_{hu}$ (m$^2$/s)','Interpreter','latex','FontSize',25);
set(gca, 'FontSize',25)
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)