Changeset 2269 for palm/trunk/SOURCE/surface_mod.f90

Timestamp:

Jun 9, 2017 11:57:32 AM (8 years ago)

Author:

suehring

Message:

Enable restarts with USM with different number of PEs; some bugfixes in new surface structure in USM; formatting adjustments and descriptions in surface_mod

File:

• palm/trunk/SOURCE/surface_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/surface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Formatting and description adjustments
+! 
+! 2256 2017-06-07 13:58:08Z suehring
 ! Enable heating at downward-facing surfaces
 ! 
@@ -33 +36 @@
 ! Description:
 ! ------------
-!> Surface module contains defines derived data structures to treat surface-
+!> Surface module defines derived data structures to treat surface-
 !> bounded grid cells. Three different types of surfaces are defined: 
-!> default surfaces, natural surfaces, and urban surfaces. Moreover, the module
-!> encompasses the initialization of near-surface grid cells, and handles reading 
-!> and writing restart data. 
+!> default surfaces, natural surfaces, and urban surfaces. The module
+!> encompasses the allocation and initialization of surface arrays, and handles 
+!> reading and writing restart data. 
 !> In addition, a further derived data structure is defined, in order to set
-!> boundary conditions at surfaces.   
+!> boundary conditions at surfaces.  
 !------------------------------------------------------------------------------!
  MODULE surface_mod
@@ -1602 +1605 @@
        IMPLICIT NONE
 
-       CHARACTER(LEN=1)             ::  dum
-
-       INTEGER(iwp)                 ::  i
-       INTEGER(iwp)                 ::  j
-       INTEGER(iwp)                 ::  l
-       INTEGER(iwp)                 ::  m
-       INTEGER(iwp), DIMENSION(0:3) ::  mm
-
-       TYPE(surf_type), DIMENSION(0:2) ::  surf_h
-       TYPE(surf_type), DIMENSION(0:3) ::  surf_v
+       CHARACTER(LEN=1)             ::  dum  !< dummy string to create output-variable name
+
+       INTEGER(iwp)                 ::  i    !< running index x-direction
+       INTEGER(iwp)                 ::  j    !< running index y-direction
+       INTEGER(iwp)                 ::  l    !< index surface type orientation
+       INTEGER(iwp)                 ::  m    !< running index for surface elements on individual surface array
+       INTEGER(iwp), DIMENSION(0:3) ::  mm   !< running index for surface elements on gathered surface array
+
+       TYPE(surf_type), DIMENSION(0:2) ::  surf_h !< gathered horizontal surfaces, contains all surface types
+       TYPE(surf_type), DIMENSION(0:3) ::  surf_v !< gathered vertical surfaces, contains all surface types
 
 !
@@ -2088 +2091 @@
 !> respective surface types within this routine. This allows e.g. changing the 
 !> surface type after reading the restart data, which might be required in case 
-!> of cyclic-filling a simulation. 
+!> of cyclic_fill mode. 
 !------------------------------------------------------------------------------!
     SUBROUTINE surface_read_restart_data( ii,                                  &
@@ -2135 +2138 @@
 
        LOGICAL                         ::  horizontal_surface !< flag indicating horizontal surfaces
-       LOGICAL                         ::  vertical_surface   !< flag indicating vertical surfaces
        LOGICAL                         ::  surf_match_def     !< flag indicating that surface element is of default type
        LOGICAL                         ::  surf_match_lsm     !< flag indicating that surface element is of natural type
        LOGICAL                         ::  surf_match_usm     !< flag indicating that surface element is of urban type
+       LOGICAL                         ::  vertical_surface   !< flag indicating vertical surfaces
 
        TYPE(surf_type), DIMENSION(0:2) ::  surf_h             !< horizontal surface type on file