Changeset 2178 for palm/trunk/SOURCE/init_pegrid.f90

Timestamp:

Mar 17, 2017 11:07:39 AM (7 years ago)

Author:

hellstea

Message:

Nesting bugfixes

File:

• palm/trunk/SOURCE/init_pegrid.f90 (modified) (1 diff)

palm/trunk/SOURCE/init_pegrid.f90

@@ -261 +261 @@
        numproc_sqr = SQRT( REAL( numprocs, KIND=wp ) )
        pdims(1)    = MAX( numproc_sqr , 1 )
-       DO  WHILE ( MOD( numprocs , pdims(1) ) /= 0 )
+       pdims(2) = numprocs / pdims(1)
+       DO  WHILE ( MOD( numprocs , pdims(1) ) /= 0  .OR.                        &
+              MOD( nx+1, pdims(1) ) /= 0  .OR.                                  &
+              MOD( ny+1, pdims(2) ) /= 0 ) 
           pdims(1) = pdims(1) - 1
+          pdims(2) = numprocs / pdims(1)
        ENDDO
-       pdims(2) = numprocs / pdims(1)
 
     ELSEIF ( npex /= -1  .AND.  npey /= -1 )  THEN