Changeset 3222 for palm/trunk/SOURCE/urban_surface_mod.f90

Timestamp:

Aug 30, 2018 1:35:35 PM (7 years ago)

Author:

suehring

Message:

Introduction of addtional surface variables indicating type and name of the surface elements; Bugfix in LSM

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Introduction of surface array for type and its name
+! 
+! 3203 2018-08-23 10:48:36Z suehring
 ! Revise bulk parameter for emissivity at ground-floor level
 ! 
@@ -886 +889 @@
        
 !--     allocate array of wall types and wall parameters
-        ALLOCATE ( surf_usm_h%surface_types(1:surf_usm_h%ns) )
+        ALLOCATE ( surf_usm_h%surface_types(1:surf_usm_h%ns)      )
+        ALLOCATE ( surf_usm_h%building_type(1:surf_usm_h%ns)      )
+        ALLOCATE ( surf_usm_h%building_type_name(1:surf_usm_h%ns) )
+        surf_usm_h%building_type = 0
+        surf_usm_h%building_type = 'none'
         DO  l = 0, 3
            ALLOCATE( surf_usm_v(l)%surface_types(1:surf_usm_v(l)%ns) )
+           ALLOCATE ( surf_usm_v(l)%building_type(1:surf_usm_v(l)%ns)      )
+           ALLOCATE ( surf_usm_v(l)%building_type_name(1:surf_usm_v(l)%ns) )
+           surf_usm_v(l)%building_type = 0
+           surf_usm_v(l)%building_type = 'none'
         ENDDO
 !
@@ -3752 +3763 @@
                                      surf_usm_h%ground_level(m) )
 !
+!--        Store building type and its name on each surface element
+           surf_usm_h%building_type(m)      = building_type
+           surf_usm_h%building_type_name(m) = building_type_name(building_type)
+!
 !--        Initialize relatvie wall- (0), green- (1) and window (2) fractions
            surf_usm_h%frac(ind_veg_wall,m)  = building_pars(ind_wall_frac,building_type)   
@@ -3869 +3884 @@
               ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
                                         surf_usm_v(l)%ground_level(m) )
-
+!
+!--           Store building type and its name on each surface element
+              surf_usm_v(l)%building_type(m)      = building_type
+              surf_usm_v(l)%building_type_name(m) = building_type_name(building_type)
 !
 !--           Initialize relatvie wall- (0), green- (1) and window (2) fractions
@@ -3999 +4017 @@
                  ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
                                            surf_usm_h%ground_level(m) )
-
+!
+!--              Store building type and its name on each surface element
+                 surf_usm_h%building_type(m)      = st
+                 surf_usm_h%building_type_name(m) = building_type_name(st)
 !
 !--              Initialize relatvie wall- (0), green- (1) and window (2) fractions
@@ -4127 +4148 @@
                     ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
                                               surf_usm_v(l)%ground_level(m) )
-
+!
+!--                 Store building type and its name on each surface element
+                    surf_usm_v(l)%building_type(m)      = st
+                    surf_usm_v(l)%building_type_name(m) = building_type_name(st)
 !
 !--                 Initialize relatvie wall- (0), green- (1) and window (2) fractions
@@ -4210 +4234 @@
            ENDDO
         ENDIF 
-
+        
 !
 !--     Level 3 - initialization via building_pars read from file