Changeset 14848

Timestamp:

2021-05-12T08:41:26+02:00 (3 years ago)

Author:

smasson

Message:

trunk: reduce minimum domain size, #2671

File:

• NEMO/trunk/src/OCE/LBC/mppini.F90 (modified) (4 diffs)

NEMO/trunk/src/OCE/LBC/mppini.F90

@@ -632 +632 @@
       klci(1:iresti      ,:) = kimax
       klci(iresti+1:knbi ,:) = kimax-1
-      IF( MINVAL(klci) < 2*i2hls ) THEN
-         WRITE(ctmp1,*) '   mpp_basesplit: minimum value of jpi must be >= ', 2*i2hls
+      IF( MINVAL(klci) < 3*khls ) THEN
+         WRITE(ctmp1,*) '   mpp_basesplit: minimum value of jpi must be >= ', 3*khls
          WRITE(ctmp2,*) '   We have ', MINVAL(klci)
-        CALL ctl_stop( 'STOP', ctmp1, ctmp2 )
+         CALL ctl_stop( 'STOP', ctmp1, ctmp2 )
       ENDIF
       IF( l_NFold ) THEN
@@ -650 +650 @@
       ENDIF
       klcj(:,1:irestj) = kjmax
-      IF( MINVAL(klcj) < 2*i2hls ) THEN
-         WRITE(ctmp1,*) '   mpp_basesplit: minimum value of jpj must be >= ', 2*i2hls
+      IF( MINVAL(klcj) < 3*khls ) THEN
+         WRITE(ctmp1,*) '   mpp_basesplit: minimum value of jpj must be >= ', 3*khls
          WRITE(ctmp2,*) '   We have ', MINVAL(klcj)
          CALL ctl_stop( 'STOP', ctmp1, ctmp2 )
@@ -729 +729 @@
       iszjref = jpiglo*jpjglo+1
       !
-      iszimin = 4*nn_hls          ! minimum size of the MPI subdomain so halos are always adressing neighbor inner domain
-      iszjmin = 4*nn_hls
+      iszimin = 3*nn_hls          ! minimum size of the MPI subdomain so halos are always adressing neighbor inner domain
+      iszjmin = 3*nn_hls
       IF( c_NFtype == 'T' )   iszjmin = MAX(iszjmin, 2+3*nn_hls)   ! V and F folding must be outside of southern halos
       IF( c_NFtype == 'F' )   iszjmin = MAX(iszjmin, 1+3*nn_hls)   ! V and F folding must be outside of southern halos
@@ -760 +760 @@
          ENDIF
       END DO
+      IF( inbimax == 0 ) THEN
+         WRITE(ctmp1,'(a,i2,a,i2)') '   mpp_ini bestpartition: Ni0glo (', Ni0glo, ') is too small to be used with nn_hls = ', nn_hls
+         CALL ctl_stop( 'STOP', ctmp1 )
+      ENDIF
+      IF( inbjmax == 0 ) THEN
+         WRITE(ctmp1,'(a,i2,a,i2)') '   mpp_ini bestpartition: Nj0glo (', Nj0glo, ') is too small to be used with nn_hls = ', nn_hls
+         CALL ctl_stop( 'STOP', ctmp1 )
+      ENDIF
 
       ! combine these 2 lists to get all possible knbi*knbj <  inbijmax