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

Timestamp:

Nov 16, 2012 3:28:04 PM (11 years ago)

Author:

raasch

Message:

bugfix for mg-solver

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Indices for arrays n.._mg start from zero due to definition of arrays f2 and
+! p2 as automatic arrays in recursive subroutine next_mg_level
 !
 ! Former revisions:
@@ -860 +861 @@
        ENDIF
 
-       ALLOCATE( grid_level_count(maximum_grid_level),                   &
-                 nxl_mg(maximum_grid_level), nxr_mg(maximum_grid_level), &
-                 nyn_mg(maximum_grid_level), nys_mg(maximum_grid_level), &
-                 nzt_mg(maximum_grid_level) )
+       ALLOCATE( grid_level_count(maximum_grid_level),                       &
+                 nxl_mg(0:maximum_grid_level), nxr_mg(0:maximum_grid_level), &
+                 nyn_mg(0:maximum_grid_level), nys_mg(0:maximum_grid_level), &
+                 nzt_mg(0:maximum_grid_level) )
 
        grid_level_count = 0
+!
+!--    Index zero required as dummy due to definition of arrays f2 and p2 in
+!--    recursive subroutine next_mg_level
+       nxl_mg(0) = 0; nxr_mg(0) = 0; nyn_mg(0) = 0; nys_mg(0) = 0; nzt_mg(0) = 0
 
        nxl_l = nxl; nxr_l = nxr; nys_l = nys; nyn_l = nyn; nzt_l = nzt