Changeset 1931 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Jun 10, 2016 12:06:59 PM (9 years ago)

Author:

suehring

Message:

Rename multigrid into multigrid_noopt and multigrid_fast into multigrid, subroutines poismg is renamed into poismg_noopt and poismg_fast_mod into poismg_mod

File:

• palm/trunk/SOURCE/init_grid.f90 (modified) (7 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename multigrid into multigrid_noopt and multigrid_fast into multigrid
 !
 ! Former revisions:
@@ -225 +225 @@
     INTEGER(iwp) ::  cyn     !< index for north canyon wall
     INTEGER(iwp) ::  cys     !< index for south canyon wall
-    INTEGER(iwp) ::  gls     !< number of lateral ghost points at total model
-                             !< domain boundaries required for multigrid solver
+    INTEGER(iwp) ::  gls     !< number of lateral ghost points at total model domain boundaries required for multigrid solver
     INTEGER(iwp) ::  i       !< index variable along x
     INTEGER(iwp) ::  ii      !< loop variable for reading topography file
@@ -242 +241 @@
 
     INTEGER(iwp), DIMENSION(:), ALLOCATABLE   ::                               &
-                     vertical_influence  !< number of vertical grid points above
-                                         !< obstacle where adjustment of near-
-                                         !< wall mixing length is required
+                     vertical_influence  !< number of vertical grid points above obstacle where adjustment of near-wall mixing length is required
                                          
-    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  corner_nl  !< index of
-                                         !< north-left corner location to limit 
-                                         !< near-wall mixing length
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  corner_nl  !< index of north-left corner location to limit near-wall mixing length
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  corner_nr  !< north-right
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  corner_sl  !< south-left
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  corner_sr  !< south-right
-    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  wall_l     !< distance to 
-                                                             !< adjacent left-facing 
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  wall_l     !< distance to adjacent left-facing 
                                                              !< wall
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  wall_n     !< north-facing
@@ -260 +254 @@
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_local  !< index for topography
                                                              !< top at cell-center
-    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_tmp    !< dummy to calculate
-                                                             !< topography indices
-                                                             !< on u- and v-grid
+    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_tmp    !< dummy to calculate topography indices on u- and v-grid
 
     LOGICAL  :: flag_set = .FALSE.  !< steering variable for advection flags
@@ -270 +262 @@
     REAL(wp) ::  dz_stretched  !< stretched vertical grid spacing
 
-    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  topo_height  !< input variable for 
-                                                           !< topography height
+    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  topo_height  !< input variable for topography height
 
     
@@ -1138 +1129 @@
 
 !
-!-- Calculate wall flag arrays for the multigrid method
-    IF ( psolver(1:9) == 'multigrid' )  THEN
+!-- Calculate wall flag arrays for the multigrid method.
+!-- Please note, wall flags are only applied in the not cache-optimized 
+!-- version.
+    IF ( psolver == 'multigrid_noopt' )  THEN
 !
 !--    Gridpoint increment of the current level
@@ -1193 +1186 @@
 !--       In case of masking method, flags are not set and multigrid method
 !--       works like FFT-solver
-          IF ( psolver == 'multigrid' .AND. .NOT. masking_method )  THEN
+          IF ( .NOT. masking_method )  THEN
 
              DO  i = nxl_l-1, nxr_l+1