%% epoch 1000
n_c=[5	10	15	20	25	30	35	40];
rmse1=[0.0065	0.0018
0.0254	0.0058
0.0061	0.0056
0.0062	0.0028
0.0071	0.0022
0.0104	0.0023
0.0133	0.0033
0.0148	0.0045
];
rmse2=[0.0107	0.0032
0.0169	0.0051
0.0178	0.0052
0.0205	0.0062
0.0187	0.0067
0.0204	0.0058
0.0396	0.0136
0.0114	0.0032
];
rmse3=[0.0213	0.0061
0.0114	0.0023
0.0188	0.0055
0.0205	0.0061
0.0206	0.0054
0.0115	0.0032
0.0246	0.0091
0.0148	0.0052
];
rmse4=[0.0302	0.0108
0.0145	0.0045
0.0240	0.0088
0.0085	0.0024
0.0167	0.0033
0.0193	0.0058
0.0126	0.0041
0.0759	0.0219
];

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);

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);

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=[5	10	15	20	25	30	35	40];
rmse1=[0.0049	0.0014
0.0033	0.0010
0.0080	0.0023
0.0022	0.0006
0.0013	0.0004
0.0044	0.0013
0.0016	0.0005
0.0021	0.0006
];
rmse2=[0.0049	0.0014
0.0033	0.0010
0.0080	0.0023
0.0022	0.0006
0.0013	0.0004
0.0044	0.0013
0.0016	0.0005
0.0021	0.0006
];
rmse3=[0.0065	0.0019
0.0038	0.0011
0.0063	0.0018
0.0077	0.0022
0.0037	0.0011
0.0036	0.0010
0.0013	0.0004
0.0089	0.0026
];
rmse4=[0.0067	0.0019
0.0116	0.0034
0.0043	0.0013
0.0016	0.0005
0.0047	0.0014
0.0018	0.0005
0.0017	0.0005
0.0017	0.0005
];

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);

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);

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=[5	10	15	20	25	30	35	40];
rmse1=[0.0041	0.0012
0.0060	0.0018
0.0015	0.0004
0.0022	0.0006
0.0068	0.0020
0.0009	0.0003
0.0016	0.0005
0.0012	0.0003
];
rmse2=[0.0043	0.0012
0.0033	0.0010
0.0048	0.0014
0.0015	0.0004
0.0024	0.0007
0.0020	0.0006
0.0015	0.0004
0.0071	0.0021
];
rmse3=[0.0073	0.0021
0.0043	0.0012
0.0064	0.0018
0.0026	0.0007
0.0037	0.0011
0.0099	0.0029
0.0031	0.0009
0.0027	0.0008
];
rmse4=[0.0195	0.0219
0.0057	0.0017
0.0107	0.0031
0.0042	0.0012
0.0030	0.0009
0.0023	0.0007
0.0025	0.0007
0.0027	0.0008
];

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);

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);

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)