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

Timestamp:

Oct 10, 2007 12:03:15 AM (17 years ago)

Author:

raasch

Message:

preliminary updates for implementing buildings in poismg

File:

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

palm/trunk/SOURCE/init_pegrid.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Allocation of wall flag arrays for multigrid solver
 ! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
 !
@@ -915 +916 @@
     ENDIF
 
+!
+!-- Allocate wall flag arrays used in the multigrid solver
+    IF ( psolver == 'multigrid' )  THEN
+
+       DO  i = maximum_grid_level, 1, -1
+
+           SELECT CASE ( i )
+
+              CASE ( 1 )
+                 ALLOCATE( wall_flags_1(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 2 )
+                 ALLOCATE( wall_flags_2(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 3 )
+                 ALLOCATE( wall_flags_3(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 4 )
+                 ALLOCATE( wall_flags_4(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 5 )
+                 ALLOCATE( wall_flags_5(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 6 )
+                 ALLOCATE( wall_flags_6(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 7 )
+                 ALLOCATE( wall_flags_7(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 8 )
+                 ALLOCATE( wall_flags_8(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 9 )
+                 ALLOCATE( wall_flags_9(nzb:nzt_mg(i)+1,         &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE ( 10 )
+                 ALLOCATE( wall_flags_10(nzb:nzt_mg(i)+1,        &
+                                        nys_mg(i)-1:nyn_mg(i)+1, &
+                                        nxl_mg(i)-1:nxr_mg(i)+1) )
+
+              CASE DEFAULT
+                 IF ( myid == 0 )  PRINT*, '+++ init_pegrid: more than 10 ', &
+                                           ' multigrid levels'
+                 CALL local_stop
+
+          END SELECT
+
+       ENDDO
+
+    ENDIF
+
  END SUBROUTINE init_pegrid