Ignore:
Timestamp:
Aug 30, 2018 1:35:35 PM (3 years ago)
Author:
suehring
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r3203 r3222  
    2828! -----------------
    2929! $Id$
     30! Introduction of surface array for type and its name
     31!
     32! 3203 2018-08-23 10:48:36Z suehring
    3033! Revise bulk parameter for emissivity at ground-floor level
    3134!
     
    886889       
    887890!--     allocate array of wall types and wall parameters
    888         ALLOCATE ( surf_usm_h%surface_types(1:surf_usm_h%ns) )
     891        ALLOCATE ( surf_usm_h%surface_types(1:surf_usm_h%ns)      )
     892        ALLOCATE ( surf_usm_h%building_type(1:surf_usm_h%ns)      )
     893        ALLOCATE ( surf_usm_h%building_type_name(1:surf_usm_h%ns) )
     894        surf_usm_h%building_type = 0
     895        surf_usm_h%building_type = 'none'
    889896        DO  l = 0, 3
    890897           ALLOCATE( surf_usm_v(l)%surface_types(1:surf_usm_v(l)%ns) )
     898           ALLOCATE ( surf_usm_v(l)%building_type(1:surf_usm_v(l)%ns)      )
     899           ALLOCATE ( surf_usm_v(l)%building_type_name(1:surf_usm_v(l)%ns) )
     900           surf_usm_v(l)%building_type = 0
     901           surf_usm_v(l)%building_type = 'none'
    891902        ENDDO
    892903!
     
    37523763                                     surf_usm_h%ground_level(m) )
    37533764!
     3765!--        Store building type and its name on each surface element
     3766           surf_usm_h%building_type(m)      = building_type
     3767           surf_usm_h%building_type_name(m) = building_type_name(building_type)
     3768!
    37543769!--        Initialize relatvie wall- (0), green- (1) and window (2) fractions
    37553770           surf_usm_h%frac(ind_veg_wall,m)  = building_pars(ind_wall_frac,building_type)   
     
    38693884              ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
    38703885                                        surf_usm_v(l)%ground_level(m) )
    3871 
     3886!
     3887!--           Store building type and its name on each surface element
     3888              surf_usm_v(l)%building_type(m)      = building_type
     3889              surf_usm_v(l)%building_type_name(m) = building_type_name(building_type)
    38723890!
    38733891!--           Initialize relatvie wall- (0), green- (1) and window (2) fractions
     
    39994017                 ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
    40004018                                           surf_usm_h%ground_level(m) )
    4001 
     4019!
     4020!--              Store building type and its name on each surface element
     4021                 surf_usm_h%building_type(m)      = st
     4022                 surf_usm_h%building_type_name(m) = building_type_name(st)
    40024023!
    40034024!--              Initialize relatvie wall- (0), green- (1) and window (2) fractions
     
    41274148                    ind_z0qh         = MERGE( ind_z0qh_gfl,         ind_z0qh_agfl,         &
    41284149                                              surf_usm_v(l)%ground_level(m) )
    4129 
     4150!
     4151!--                 Store building type and its name on each surface element
     4152                    surf_usm_v(l)%building_type(m)      = st
     4153                    surf_usm_v(l)%building_type_name(m) = building_type_name(st)
    41304154!
    41314155!--                 Initialize relatvie wall- (0), green- (1) and window (2) fractions
     
    42104234           ENDDO
    42114235        ENDIF
    4212 
     4236       
    42134237!
    42144238!--     Level 3 - initialization via building_pars read from file
Note: See TracChangeset for help on using the changeset viewer.