Changeset 4264 for palm


Ignore:
Timestamp:
Oct 15, 2019 4:00:23 PM (4 years ago)
Author:
scharf
Message:

corrected error message string for incompatible combinations of nx, ny or nz with pdims(1) or pdims(2)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_pegrid.f90

    r4241 r4264  
    2525! -----------------
    2626! $Id$
     27! corrected error message string
     28!
     29! 4241 2019-09-27 06:32:47Z raasch
    2730! Check added to ensure that subdomain grid has at least the size as given by the number
    2831! of ghost points
     
    300303    IF ( MOD( nx+1 , pdims(1) ) /= 0 )  THEN
    301304       WRITE( message_string, * ) 'x-direction: gridpoint number (',nx+1,') ', &
    302                                'is not an& integral divisor of the number ',    &
    303                                'of processors (', pdims(1),')'
     305                                  'is not an& integral multiple of the number',&
     306                                  ' of processors (',pdims(1),')'
    304307       CALL message( 'init_pegrid', 'PA0225', 1, 2, 0, 6, 0 )
    305308    ELSE
    306309       nnx  = ( nx + 1 ) / pdims(1)
    307     ENDIF   
     310    ENDIF
    308311
    309312!
     
    318321    IF ( MOD( ny+1 , pdims(2) ) /= 0 )  THEN
    319322       WRITE( message_string, * ) 'y-direction: gridpoint number (',ny+1,') ', &
    320                            'is not an& integral divisor of the number of',      &
    321                            'processors (', pdims(2),')'
     323                                  'is not an& integral multiple of the number',&
     324                                  ' of processors (',pdims(2),')'
    322325       CALL message( 'init_pegrid', 'PA0227', 1, 2, 0, 6, 0 )
    323326    ELSE
    324327       nny  = ( ny + 1 ) / pdims(2)
    325     ENDIF   
     328    ENDIF
    326329
    327330!
     
    363366       IF ( pdims(2) /= 1 )  THEN
    364367          IF ( MOD( nz , pdims(1) ) /= 0 )  THEN
    365              WRITE( message_string, * ) 'transposition z --> x:',              &
    366                        '& nz=',nz,' is not an integral divisior of pdims(1)=', &
    367                                                                    pdims(1)
     368             WRITE( message_string, * ) 'transposition z --> x:& ',              &
     369                                        'nz=',nz,' is not an integral multiple ',&
     370                                        'of pdims(1)=',pdims(1)
    368371             CALL message( 'init_pegrid', 'PA0230', 1, 2, 0, 6, 0 )
    369372          ENDIF
     
    381384
    382385
    383     IF ( psolver == 'poisfft' )  THEN 
     386    IF ( psolver == 'poisfft' )  THEN
    384387!
    385388!--    2. transposition  x --> y
    386389       IF ( MOD( nx+1 , pdims(2) ) /= 0 )  THEN
    387           WRITE( message_string, * ) 'transposition x --> y:',                 &
    388                             '& nx+1=',nx+1,' is not an integral divisor of ',  &
    389                             'pdims(2)=',pdims(2)
     390          WRITE( message_string, * ) 'transposition x --> y:& ',              &
     391                                     'nx+1=',nx+1,' is not an integral ',     &
     392                                     'multiple of pdims(2)=',pdims(2)
    390393          CALL message( 'init_pegrid', 'PA0231', 1, 2, 0, 6, 0 )
    391394       ENDIF
     
    399402       sendrecvcount_xy = nnx_y * nny_x * nnz_y
    400403!
    401 !--    3. transposition  y --> z 
     404!--    3. transposition  y --> z
    402405!--    (ELSE:  x --> y  in case of 1D-decomposition along x)
    403406       nxl_z = nxl_y
     
    414417!--       along x, except that the uptream-spline method is switched on
    415418          IF ( MOD( ny+1 , pdims(1) ) /= 0 )  THEN
    416              WRITE( message_string, * ) 'transposition y --> z:',              &
    417                                '& ny+1=',ny+1,' is not an integral divisor of',&
    418                                ' pdims(1)=',pdims(1)
     419             WRITE( message_string, * ) 'transposition y --> z:& ',            &
     420                                        'ny+1=',ny+1,' is not an integral ',   &
     421                                        'multiple of pdims(1)=',pdims(1)
    419422             CALL message( 'init_pegrid', 'PA0232', 1, 2, 0, 6, 0 )
    420423          ENDIF
     
    425428!--       This condition must be fulfilled for a 1D-decomposition along x
    426429          IF ( MOD( ny+1 , pdims(1) ) /= 0 )  THEN
    427              WRITE( message_string, * ) 'transposition x --> y:',              &
    428                                '& ny+1=',ny+1,' is not an integral divisor of',&
    429                                ' pdims(1)=',pdims(1)
     430             WRITE( message_string, * ) 'transposition x --> y:& ',            &
     431                                        'ny+1=',ny+1,' is not an integral ',   &
     432                                        'multiple of pdims(1)=',pdims(1)
    430433             CALL message( 'init_pegrid', 'PA0233', 1, 2, 0, 6, 0 )
    431434          ENDIF
     
    440443       IF ( MOD( nz, pdims(2) ) /= 0 )  THEN
    441444          WRITE( message_string, * ) 'direct transposition z --> y (needed ',  &
    442                     'for spectra): nz=',nz,' is not an integral divisor of ',  &
    443                     'pdims(2)=',pdims(2)
     445                                     'for spectra):& nz=',nz,' is not an ',    &
     446                                     'integral multiple of pdims(2)=',pdims(2)
    444447          CALL message( 'init_pegrid', 'PA0234', 1, 2, 0, 6, 0 )
    445448       ELSE
Note: See TracChangeset for help on using the changeset viewer.