Changeset 3057

Timestamp:

Jun 5, 2018 9:03:41 AM (6 years ago)

Author:

raasch

Message:

bugfix: check that nz is even in case that optimized multigrid is used

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix: check that nz is even in case that optimized multigrid is used
+! 
+! 3049 2018-05-29 13:52:36Z Giersch
 ! Error messages revised
 ! 
@@ -911 +914 @@
        IF ( TRIM( psolver ) /= 'multigrid_noopt' )  THEN
           IF ( MOD( k, 2 ) /= 0 )  mg_levels_z = mg_levels_z - 1
+!
+!--       An odd value of nz does not work. The finest level must have an even
+!--       value.
+          IF (  mg_levels_z == 0 )  THEN
+             message_string = 'optimized multigrid method requires nz to be even'
+             CALL message( 'init_pegrid', 'PA0028', 1, 2, 0, 6, 0 )
+          ENDIF
        ENDIF
 
@@ -949 +959 @@
 
           ELSE
+
              mg_switch_to_pe0_level_l = 0
              maximum_grid_level_l = maximum_grid_level