%% epoch 1000
n_c=[30	40	50	60	70	80];
rmse1=[0.0125	0.0285
0.0089	0.0203
0.0096	0.0212
0.0076	0.0162
0.0061	0.0134
0.0074	0.0169
];
rmse2=[0.0065	0.0145
0.0073	0.0172
0.0129	0.0265
0.0073	0.017
0.0075	0.0156
0.0081	0.018
];
rmse3=[0.0071	0.0166
0.0073	0.0169
0.0065	0.0152
0.008	0.0174
0.0127	0.0275
0.0066	0.0146
];
rmse4=[0.0078	0.0178
0.0079	0.0172
0.007	0.0156
0.0061	0.0137
0.0078	0.0175
0.0072	0.0163
];

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);
%ylim([0, 0.3])
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);
%ylim([0, 0.25])
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)
%% epoch 2000
n_c=[30	40	50	60	70	80];
rmse1=[0.0118	0.028
0.0067	0.013
0.0071	0.0151
0.007	0.0149
0.0063	0.0143
0.0068	0.0137
];
rmse2=[0.0081	0.0182
0.0059	0.0132
0.0119	0.0277
0.0073	0.0154
0.0068	0.0146
0.0071	0.0151
];
rmse3=[0.0064	0.0139
0.0059	0.0127
0.0058	0.013
0.0069	0.0154
0.0123	0.0286
0.0079	0.0186
];
rmse4=[0.0069	0.0151
0.0067	0.014
0.0085	0.018
0.0055	0.0111
0.0075	0.0167
0.0058	0.0139
];

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);
%ylim([0, 0.3])
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);
%ylim([0, 0.25])
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)
%% epoch 3000
n_c=[30	40	50	60	70	80];
rmse1=[0.0119	0.0270
0.0073	0.0177
0.0089	0.0218
0.0059	0.012
0.0072	0.0153
0.0086	0.0206
];
rmse2=[0.0066	0.0138
0.0063	0.0136
0.0116	0.0283
0.0072	0.0163
0.0064	0.0141
0.0063	0.012
];
rmse3=[0.0065	0.0137
0.0071	0.0169
0.006	0.0121
0.006	0.0137
0.0136	0.0287
0.0062	0.0119
];
rmse4=[0.0065	0.0159
0.0078	0.0193
0.0069	0.0141
0.0067	0.0145
0.0062	0.0133
0.0068	0.0155
];

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);
%ylim([0, 0.3])
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);
%ylim([0, 0.25])
set(gca, 'FontSize',25)
title('Water Discharge','Interpreter','latex','FontSize',25)