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

Timestamp:

May 28, 2018 7:55:41 AM (6 years ago)

Author:

Giersch

Message:

Code adjusted according to coding standards, renamed namelists, error messages revised until PA0347, output CASE 108 disabled

File:

• palm/trunk/SOURCE/init_pegrid.f90 (modified) (13 diffs)

palm/trunk/SOURCE/init_pegrid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Error messages revised
+! 
+! 2938 2018-03-27 15:52:42Z suehring
 ! - No checks for domain decomposition in case of turbulence generator 
 !  (is done in stg module)
@@ -341 +344 @@
 !--    must be equal to the number of PEs available to the job
        IF ( ( npex * npey ) /= numprocs )  THEN
-          WRITE( message_string, * ) 'number of PEs of the prescribed ',   &
-              'topology (', npex*npey,') does not match & the number of ', &
+          WRITE( message_string, * ) 'number of PEs of the prescribed ',       &
+              'topology (', npex*npey,') does not match the number of ',       &
               'PEs available to the job (', numprocs, ')'
           CALL message( 'init_pegrid', 'PA0221', 1, 2, 0, 6, 0 )
@@ -353 +356 @@
 !--    If the processor topology is prescribed by the user, the number of
 !--    PEs must be given in both directions
-       message_string = 'if the processor topology is prescribed by th' //  &
-                'e user& both values of "npex" and "npey" must be given' // &
+       message_string = 'if the processor topology is prescribed by th' //     &
+                'e user both values of "npex" and "npey" must be given' //     &
                 ' in the &NAMELIST-parameter file'
        CALL message( 'init_pegrid', 'PA0222', 1, 2, 0, 6, 0 )
@@ -450 +453 @@
 !
 !-- Calculate array bounds along x-direction for every PE.
-    ALLOCATE( nxlf(0:pdims(1)-1), nxrf(0:pdims(1)-1), nynf(0:pdims(2)-1), &
+    ALLOCATE( nxlf(0:pdims(1)-1), nxrf(0:pdims(1)-1), nynf(0:pdims(2)-1),      &
               nysf(0:pdims(2)-1) )
 
     IF ( MOD( nx+1 , pdims(1) ) /= 0 )  THEN
-       WRITE( message_string, * ) 'x-direction: gridpoint number (',nx+1,') ',&
-                               'is not an& integral divisor of the number ',  &
-                               'processors (', pdims(1),')'
+       WRITE( message_string, * ) 'x-direction: gridpoint number (',nx+1,') ', &
+                               'is not an integral divisor of the number ',    &
+                               'of processors (', pdims(1),')'
        CALL message( 'init_pegrid', 'PA0225', 1, 2, 0, 6, 0 )
     ELSE
        nnx  = ( nx + 1 ) / pdims(1)
-       IF ( nnx*pdims(1) - ( nx + 1) > nnx )  THEN
-          WRITE( message_string, * ) 'x-direction: nx does not match the',    & 
-                       'requirements given by the number of PEs &used',       &
-                       '& please use nx = ', nx - ( pdims(1) - ( nnx*pdims(1) &
-                                      - ( nx + 1 ) ) ), ' instead of nx =', nx
-          CALL message( 'init_pegrid', 'PA0226', 1, 2, 0, 6, 0 )
-       ENDIF
     ENDIF    
 
@@ -480 +476 @@
     IF ( MOD( ny+1 , pdims(2) ) /= 0 )  THEN
        WRITE( message_string, * ) 'y-direction: gridpoint number (',ny+1,') ', &
-                           'is not an& integral divisor of the number of',     &
+                           'is not an integral divisor of the number of',      &
                            'processors (', pdims(2),')'
        CALL message( 'init_pegrid', 'PA0227', 1, 2, 0, 6, 0 )
     ELSE
        nny  = ( ny + 1 ) / pdims(2)
-       IF ( nny*pdims(2) - ( ny + 1) > nny )  THEN
-          WRITE( message_string, * ) 'y-direction: ny does not match the',    &
-                       'requirements given by the number of PEs &used ',      &
-                       '& please use ny = ', ny - ( pdims(2) - ( nnx*pdims(2) &
-                                     - ( ny + 1 ) ) ), ' instead of ny =', ny
-          CALL message( 'init_pegrid', 'PA0228', 1, 2, 0, 6, 0 )
-       ENDIF
     ENDIF    
 
@@ -533 +522 @@
           IF ( MOD( nz , pdims(1) ) /= 0 )  THEN
              WRITE( message_string, * ) 'transposition z --> x:',              &
-                       '&nz=',nz,' is not an integral divisior of pdims(1)=',  &
+                       ' nz=',nz,' is not an integral divisior of pdims(1)=',  &
                                                                    pdims(1)
              CALL message( 'init_pegrid', 'PA0230', 1, 2, 0, 6, 0 )
@@ -555 +544 @@
        IF ( MOD( nx+1 , pdims(2) ) /= 0 )  THEN
           WRITE( message_string, * ) 'transposition x --> y:',                 &
-                            '&nx+1=',nx+1,' is not an integral divisor of ',   &
+                            ' nx+1=',nx+1,' is not an integral divisor of ',   &
                             'pdims(2)=',pdims(2)
           CALL message( 'init_pegrid', 'PA0231', 1, 2, 0, 6, 0 )
@@ -584 +573 @@
           IF ( MOD( ny+1 , pdims(1) ) /= 0 )  THEN
              WRITE( message_string, * ) 'transposition y --> z:',              &
-                               '& ny+1=',ny+1,' is not an integral divisor of',&
+                               ' ny+1=',ny+1,' is not an integral divisor of', &
                                ' pdims(1)=',pdims(1)
              CALL message( 'init_pegrid', 'PA0232', 1, 2, 0, 6, 0 )
@@ -595 +584 @@
           IF ( MOD( ny+1 , pdims(1) ) /= 0 )  THEN
              WRITE( message_string, * ) 'transposition x --> y:',              &
-                               '& ny+1=',ny+1,' is not an integral divisor of',&
+                               ' ny+1=',ny+1,' is not an integral divisor of', &
                                ' pdims(1)=',pdims(1)
              CALL message( 'init_pegrid', 'PA0233', 1, 2, 0, 6, 0 )
@@ -609 +598 @@
        IF ( MOD( nz, pdims(2) ) /= 0 )  THEN
           WRITE( message_string, * ) 'direct transposition z --> y (needed ',  &
-                    'for spectra):& nz=',nz,' is not an integral divisor of ', &
+                    'for spectra): nz=',nz,' is not an integral divisor of ',  &
                     'pdims(2)=',pdims(2)
           CALL message( 'init_pegrid', 'PA0234', 1, 2, 0, 6, 0 )
@@ -974 +963 @@
 !
 !--          Check pre-defined value and reset to default, if neccessary
-             IF ( mg_switch_to_pe0_level < mg_switch_to_pe0_level_l  .OR.  &
+             IF ( mg_switch_to_pe0_level < mg_switch_to_pe0_level_l  .OR.      &
                   mg_switch_to_pe0_level >= maximum_grid_level_l )  THEN
-                message_string = 'mg_switch_to_pe0_level ' // &
+                message_string = 'mg_switch_to_pe0_level ' //                  &
                                  'out of range and reset to 0'
                 CALL message( 'init_pegrid', 'PA0235', 0, 1, 0, 6, 0 )
@@ -1038 +1027 @@
              subdomain_size = ( nxr - nxl + 2 * nbgp + 1 ) * &
                               ( nyn - nys + 2 * nbgp + 1 ) * ( nzt - nzb + 2 )
-             gathered_size  = ( nxr_l - nxl_l + 3 ) * ( nyn_l - nys_l + 3 ) * &
+             gathered_size  = ( nxr_l - nxl_l + 3 ) * ( nyn_l - nys_l + 3 ) *  &
                               ( nzt_l - nzb + 2 )
 
 #else
-             message_string = 'multigrid gather/scatter impossible ' // &
+             message_string = 'multigrid gather/scatter impossible ' //        &
                           'in non parallel mode'
              CALL message( 'init_pegrid', 'PA0237', 1, 2, 0, 6, 0 )
@@ -1063 +1052 @@
 
 !
-!--    Temporary problem: Currently calculation of maxerror iin routine poismg crashes
+!--    Temporary problem: Currently calculation of maxerror in routine poismg crashes
 !--    if grid data are collected on PE0 already on the finest grid level.
 !--    To be solved later.