Anticipation of error related to hydrological input which are not well masked over ocean.

F90subroutines/routing_reg.f90

@@ -261,7 +261,8 @@ CONTAINS
           !
           IF ( topoindex(ib, ip) < 0 .OR. topoindex(ib, ip) > 10000 ) THEN
              WRITE(numout,*) "There is an issue with Topoindex : ", topoindex(ib, ip), min_topoind
-             WRITE(numout,*) "Check other variables : ", trip(ib, ip), basins(ib, ip), sub_area(ib, ip)
+             WRITE(numout,*) "Check other variables : ", trip(ib, ip), basins(ib, ip), sub_area(ib, ip) 
+             WRITE(numout,*) "rdz,rlen:", rdz(ib, ip), rlen(ib, ip)
           ENDIF
           trip_bx(iloc, jloc) = NINT(trip(ib, ip))
           basin_bx(iloc, jloc) = NINT(basins(ib, ip))
@@ -314,7 +315,7 @@ CONTAINS
           ! Compute the destination of the flow on the high resolution grid (if possible !) to see if it still within the points
           ! belonging to the grid we are working on.
           !
-          IF ( trip_bx(ip,jp) < 97 ) THEN
+          IF (( trip_bx(ip,jp) < 97 ) .AND. ( trip_bx(ip,jp) > 0 )) THEN
              !
              ipp = ip+inc(trip_bx(ip,jp),1)
              jpp = jp+inc(trip_bx(ip,jp),2)
@@ -339,6 +340,11 @@ CONTAINS
                    ENDIF
                 ENDIF
              ENDIF
+          ELSE IF ( trip_bx(ip,jp) .LE. 0 ) THEN
+             WRITE(numout,*) "We have a value of trip:", trip_bx(ip,jp)
+             WRITE(numout,*) "The variables in the hydrological input may not be correctly masked. Please check it."
+             CALL ipslerr_p(3,'routing_reg_getgrid','Negative value of trip->mask the original variables over ocean','','')
+
           ELSE IF ( trip_bx(ip,jp) > 100 .AND. trip_bx(ip,jp) < 109 ) THEN
              WRITE(numout,*) 'WARNING : Point flows our of routing.nc file '
              WRITE(numout,*) 'WARNING : Point : basin = ', basin_bx(ip,jp)