Index: palm/trunk/SOURCE/init_pegrid.f90
===================================================================
--- palm/trunk/SOURCE/init_pegrid.f90	(revision 3543)
+++ palm/trunk/SOURCE/init_pegrid.f90	(revision 3552)
@@ -20,5 +20,5 @@
 ! Current revisions:
 ! ------------------
-! 
+! variables documented
 ! 
 ! Former revisions:
@@ -289,42 +289,42 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  i                        !<
-    INTEGER(iwp) ::  id_inflow_l              !<
+    INTEGER(iwp) ::  i                        !< running index over number of processors or number of multigrid level
+    INTEGER(iwp) ::  id_inflow_l              !< ID indicating processors located at the left inflow boundary
     INTEGER(iwp) ::  id_outflow_l             !< local value of id_outflow
     INTEGER(iwp) ::  id_outflow_source_l      !< local value of id_outflow_source
-    INTEGER(iwp) ::  id_recycling_l           !<
+    INTEGER(iwp) ::  id_recycling_l           !< ID indicating processors located at the recycling plane
     INTEGER(iwp) ::  id_stg_left_l            !< left lateral boundary local core id in case of turbulence generator  
     INTEGER(iwp) ::  id_stg_north_l           !< north lateral boundary local core id in case of turbulence generator  
     INTEGER(iwp) ::  id_stg_right_l           !< right lateral boundary local core id in case of turbulence generator  
     INTEGER(iwp) ::  id_stg_south_l           !< south lateral boundary local core id in case of turbulence generator  
-    INTEGER(iwp) ::  ind(5)                   !<
-    INTEGER(iwp) ::  j                        !<
-    INTEGER(iwp) ::  k                        !<
-    INTEGER(iwp) ::  maximum_grid_level_l     !<
-    INTEGER(iwp) ::  mg_levels_x              !<
-    INTEGER(iwp) ::  mg_levels_y              !<
-    INTEGER(iwp) ::  mg_levels_z              !<
-    INTEGER(iwp) ::  mg_switch_to_pe0_level_l !<
-    INTEGER(iwp) ::  nnx_y                    !<
-    INTEGER(iwp) ::  nny_x                    !<
-    INTEGER(iwp) ::  nny_z                    !<
-    INTEGER(iwp) ::  nnz_x                    !<
-    INTEGER(iwp) ::  nnz_y                    !<
-    INTEGER(iwp) ::  numproc_sqr              !<
-    INTEGER(iwp) ::  nxl_l                    !<
-    INTEGER(iwp) ::  nxr_l                    !<
-    INTEGER(iwp) ::  nyn_l                    !<
-    INTEGER(iwp) ::  nys_l                    !<
-    INTEGER(iwp) ::  nzb_l                    !<
-    INTEGER(iwp) ::  nzt_l                    !<
-!$  INTEGER(iwp) ::  omp_get_num_threads      !<
-
-    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  ind_all !<
-    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  nxlf    !<
-    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  nxrf    !<
-    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  nynf    !<
-    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  nysf    !<
-
-    INTEGER(iwp), DIMENSION(2) ::  pdims_remote         !<
+    INTEGER(iwp) ::  ind(5)                   !< array containing the subdomain bounds
+    INTEGER(iwp) ::  j                        !< running index, used for various loops
+    INTEGER(iwp) ::  k                        !< number of vertical grid points in different multigrid level
+    INTEGER(iwp) ::  maximum_grid_level_l     !< maximum number of grid level without switching to PE 0
+    INTEGER(iwp) ::  mg_levels_x              !< maximum number of grid level allowed along x-direction
+    INTEGER(iwp) ::  mg_levels_y              !< maximum number of grid level allowed along y-direction
+    INTEGER(iwp) ::  mg_levels_z              !< maximum number of grid level allowed along z-direction
+    INTEGER(iwp) ::  mg_switch_to_pe0_level_l !< maximum number of grid level with switching to PE 0
+    INTEGER(iwp) ::  nnx_y                    !< quotient of number of grid points along x-direction and number of PEs used along y-direction
+    INTEGER(iwp) ::  nny_x                    !< quotient of number of grid points along y-direction and number of PEs used along x-direction
+    INTEGER(iwp) ::  nny_z                    !< quotient of number of grid points along y-direction and number of PEs used along x-direction
+    INTEGER(iwp) ::  nnz_x                    !< quotient of number of grid points along z-direction and number of PEs used along x-direction
+    INTEGER(iwp) ::  nnz_y                    !< quotient of number of grid points along z-direction and number of PEs used along x-direction
+    INTEGER(iwp) ::  numproc_sqr              !< square root of the number of processors
+    INTEGER(iwp) ::  nxl_l                    !< lower index bound along x-direction on subdomain and different multigrid level
+    INTEGER(iwp) ::  nxr_l                    !< upper index bound along x-direction on subdomain and different multigrid level
+    INTEGER(iwp) ::  nyn_l                    !< lower index bound along y-direction on subdomain and different multigrid level
+    INTEGER(iwp) ::  nys_l                    !< upper index bound along y-direction on subdomain and different multigrid level
+    INTEGER(iwp) ::  nzb_l                    !< lower index bound along z-direction on subdomain and different multigrid level
+    INTEGER(iwp) ::  nzt_l                    !< upper index bound along z-direction on subdomain and different multigrid level
+!$  INTEGER(iwp) ::  omp_get_num_threads      !< number of OpenMP threads
+
+    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  ind_all !< dummy array containing index bounds on subdomain, used for gathering
+    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  nxlf    !< lower index bound allong x-direction for every PE
+    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  nxrf    !< upper index bound allong x-direction for every PE
+    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  nynf    !< lower index bound allong y-direction for every PE
+    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  nysf    !< lower index bound allong y-direction for every PE
+
+    INTEGER(iwp), DIMENSION(2) ::  pdims_remote         !< number of PEs used for coupled model (only in atmospher-ocean coupling)
     INTEGER(iwp)               ::  lcoord(2)            !< PE coordinates of left neighbor along x and y
     INTEGER(iwp)               ::  rcoord(2)            !< PE coordinates of right neighbor along x and y
