Changeset 2197

Timestamp:

Mar 24, 2017 2:25:00 AM (7 years ago)

Author:

raasch

Message:

bugfix: do not allow odd values for nz at the coarsest grid level in case of optimized multigrid solver

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! bugfix: do not allow odd values for nz at the coarsest grid level in case of
+! optimized multigrid solver
 ! 
 ! Former revisions:
@@ -862 +863 @@
           mg_levels_z = mg_levels_z + 1
        ENDDO
+!
+!--    The optimized MG-solver does not allow odd values for nz at the coarsest
+!--    grid level
+       IF ( TRIM( psolver ) /= 'multigrid_noopt' )  THEN
+          IF ( MOD( k, 2 ) /= 0 )  mg_levels_z = mg_levels_z - 1
+       ENDIF
 
        maximum_grid_level = MIN( mg_levels_x, mg_levels_y, mg_levels_z )