function [mean_change] = func_cal_domain_mean(soil, soil_diff)

soil_domain1 = soil(:,:,1,5);     % no routing
soil_domain2 = soil(:,:,1,5);     % no routing
soil_domain3 = soil(:,:,1,5);     % no routing

soil_diff_domain1 = soil_diff(:,:,2,5) ; % routing 100m - routing 250m
soil_diff_domain2 = soil_diff(:,:,3,5) ; % routing 100m - routing 500m
soil_diff_domain3 = soil_diff(:,:,4,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
