Changeset 13298

Timestamp:

2020-07-13T12:12:33+02:00 (3 years ago)

Author:

dancopsey

Message:

Improve error message when there is a zero river area for one of the rivers.

File:

• branches/UKMO/dev_r5518_GO6_fix_zemp_ice_10681/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_rnf_1d.F90 (modified) (3 diffs)

branches/UKMO/dev_r5518_GO6_fix_zemp_ice_10681/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_rnf_1d.F90

@@ -61 +61 @@
       INTEGER                                   ::   max_river
       REAL(wp), POINTER, DIMENSION(:,:)         ::   river_number        ! 2D array containing the river outflow numbers
+      REAL(wp)                                  ::   min_river_area      ! Mininum river area used for debugging
+      INTEGER                                   ::   min_river_number    ! Number of minimum river area used for debugging
       
       NAMELIST/nam_cpl_rnf_1d/file_riv_number, nn_cpl_river, ln_print_river_info
@@ -117 +119 @@
       ALLOCATE( rivers%river_area_c( nn_cpl_river ) )
       rivers%river_area_c(:) = CMPLX( 0.e0, 0.e0, wp )
-      DO ii = nldi, nlei      
-        DO jj = nldj, nlej
+      DO ii = 1, jpi      
+        DO jj = 1, jpj
           IF ( tmask_i(ii,jj) > 0.5 ) THEN  ! This makes sure we are not at a duplicated point (at north fold or east-west cyclic point)
             IF ( rivers%river_number(ii,jj) > 0 .AND. rivers%river_number(ii,jj) <= nn_cpl_river ) THEN
@@ -139 +141 @@
       IF ( ln_print_river_info ) THEN
         WRITE(numout,*) 'Area of rivers 1 to 10 are ',rivers%river_area(1:10)
+      END IF
+
+      IF ( MINVAL(rivers%river_area(1:nn_cpl_river)) <= 0 ) THEN
+         min_river_area=9E15
+         min_river_number=0
+         DO rr = 1, nn_cpl_river
+            IF ( rivers%river_area(rr) < min_river_area ) THEN
+               min_river_area = rivers%river_area(rr)
+               min_river_number = rr
+            END IF
+         END DO
+         WRITE(numout,*) 'Minimum river area = ',min_river_area
+         WRITE(numout,*) 'First river number with minimum area = ',min_river_number
+         CALL ctl_stop( 'There is at least one river with a river outflow area of zero. Please check your file_riv_number file that all the allocated river numbers are at ocean points as defined by the bathymetry file.')
       END IF