Fix the bug when calculating topographic index and put comment to explain the idea of routing_reg_divbas