Correction in topoindex_stream

F90subroutines/routing_reg.f90

@@ -2062,16 +2062,13 @@ SUBROUTINE routing_reg_divbas_divide(nb_htu, nbv,nbi,nbj, ijdimmax, tbname, tsz,
    IF (cnt .GT. 1) THEN
       dorog = MAXVAL(main_orog(:cnt)) - MINVAL(main_orog(:cnt))
       dorog = MAX(0.1, dorog)
+      distance = SUM(main_len(:cnt))
    ELSE
       ! if the main river has only one pixel we use directly rdz 
       dorog = MAX(0.1, rdz(il,jl))
+      distance = main_len(1)
    END IF
   
-   ! For the distance : sum rlen along the main river 
-   DO l = 1,cnt
-      distance = SUM(main_len(:cnt))
-   END DO 
-
    topoindex_stream = SQRT(distance **3 / dorog)
 !
 END SUBROUTINE mainstream_topoindex