Ignore:
Timestamp:
Nov 7, 2018 2:45:23 PM (5 years ago)
Author:
suehring
Message:

Disable initialization of building roofs with ground-floor-level properties since this causes strong oscillations of surface temperature during the spin-up.

File:
1 edited

Legend:

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

    r3469 r3502  
    2828! -----------------
    2929! $Id$
     30! Disable initialization of building roofs with ground-floor-level properties,
     31! since this causes strong oscillations of surface temperature during the
     32! spinup.
     33!
     34! 3469 2018-10-30 20:05:07Z kanani
    3035! Add missing PUBLIC variables for new indoor model
    3136!
     
    10161021    REAL(wp)                                       :: window_inner_temperature = 295.0_wp  !< temperature of the inner window surface (~22 degrees C) (K)
    10171022
    1018         REAL(wp)                                       ::   m_total = 0.0_wp !< weighted total water content of the soil (m3/m3)
     1023    REAL(wp)                                       ::   m_total = 0.0_wp !< weighted total water content of the soil (m3/m3)
    10191024    INTEGER(iwp)                                   ::   soil_type
    10201025   
     
    49424947!--     use terrain height array from file, if available. This flag is later used
    49434948!--     to control initialization of surface attributes.
     4949!--     Todo: for the moment disable initialization of building roofs with
     4950!--     ground-floor-level properties.
    49444951        surf_usm_h%ground_level = .FALSE.
    4945         DO  m = 1, surf_usm_h%ns
    4946            i = surf_usm_h%i(m)
    4947            j = surf_usm_h%j(m)
    4948            k = surf_usm_h%k(m)
    4949 !
    4950 !--        Get local ground level. If no ground level is given in input file,
    4951 !--        use default value.
    4952            ground_floor_level_l = ground_floor_level
    4953            IF ( building_pars_f%from_file )  THEN
    4954               IF ( building_pars_f%pars_xy(ind_gflh,j,i) /=                    &
    4955                    building_pars_f%fill )  &
    4956                  ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)         
    4957            ENDIF
    4958 !
    4959 !--        Determine height of surface element above ground level
    4960            IF (  terrain_height_f%from_file )  THEN
    4961               z_agl = zw(k) - terrain_height_f%var(j,i)
    4962            ELSE
    4963               z_agl = zw(k)
    4964            ENDIF
    4965 !
    4966 !--        Set flag for ground level
    4967            IF ( z_agl <= ground_floor_level_l )                                &
    4968               surf_usm_h%ground_level(m) = .TRUE.
    4969         ENDDO
     4952!         DO  m = 1, surf_usm_h%ns
     4953!            i = surf_usm_h%i(m)
     4954!            j = surf_usm_h%j(m)
     4955!            k = surf_usm_h%k(m)
     4956! !
     4957! !--        Get local ground level. If no ground level is given in input file,
     4958! !--        use default value.
     4959!            ground_floor_level_l = ground_floor_level
     4960!            IF ( building_pars_f%from_file )  THEN
     4961!               IF ( building_pars_f%pars_xy(ind_gflh,j,i) /=                    &
     4962!                    building_pars_f%fill )  &
     4963!                  ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)         
     4964!            ENDIF
     4965! !
     4966! !--        Determine height of surface element above ground level
     4967!            IF (  terrain_height_f%from_file )  THEN
     4968!               z_agl = zw(k) - terrain_height_f%var(j,i)
     4969!            ELSE
     4970!               z_agl = zw(k)
     4971!            ENDIF
     4972! !
     4973! !--        Set flag for ground level
     4974!            IF ( z_agl <= ground_floor_level_l )                                &
     4975!               surf_usm_h%ground_level(m) = .TRUE.
     4976!         ENDDO
    49704977
    49714978        DO  l = 0, 3
Note: See TracChangeset for help on using the changeset viewer.