function [mean_change] = func_cal_domain_mean(soil, soil_diff)

soil_domain1 = soil(:,:,1,5);     % routing 250 m
soil_domain2 = soil(:,:,1,5);     % routing 500 m
soil_domain3 = soil(:,:,1,5);     % routing 1000 m

soil_diff_domain1 = soil_diff(:,:,6,5) ; % routing 100m - routing 250m
soil_diff_domain2 = soil_diff(:,:,7,5) ; % routing 100m - routing 500m
soil_diff_domain3 = soil_diff(:,:,8,5) ; % routing 100m - routing 1000m


mean_change(:,:,1) = (soil_diff_domain1)./(soil_domain1)*100;
mean_change(:,:,2) = (soil_diff_domain2)./(soil_domain2)*100;
mean_change(:,:,3) = (soil_diff_domain3)./(soil_domain3)*100;

end