Improvement which supresses the singularity in the FV case

F90subroutines/routing_reg.f90

@@ -2805,7 +2805,7 @@ SUBROUTINE routing_reg_end_truncate(nbpt, nbasmax, gridarea, contfrac, gridcente
          basin_cg(ib,1,2) = gridcenters(ib,2)
          ! To define in function of the minimum topoindex on the map
          ! (to revise)
-         basin_topoind(ib,1) = topoindexmin
+         basin_topoind(ib,1) = SQRT(hydro_meanlen**3./(hydro_meandz*1.0e6))
          basin_rlen(ib,1) = hydro_meanlen
          basin_rdz(ib,1) = hydro_meandz
          basin_id(ib,1) = 0