Changeset 3502 for palm

Timestamp:

Nov 7, 2018 2:45:23 PM (6 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:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Disable initialization of building roofs with ground-floor-level properties, 
+! since this causes strong oscillations of surface temperature during the 
+! spinup.
+! 
+! 3469 2018-10-30 20:05:07Z kanani
 ! Add missing PUBLIC variables for new indoor model
 ! 
@@ -1016 +1021 @@
     REAL(wp)                                       :: window_inner_temperature = 295.0_wp  !< temperature of the inner window surface (~22 degrees C) (K)
 
-        REAL(wp)                                       ::   m_total = 0.0_wp !< weighted total water content of the soil (m3/m3)
+    REAL(wp)                                       ::   m_total = 0.0_wp !< weighted total water content of the soil (m3/m3)
     INTEGER(iwp)                                   ::   soil_type
     
@@ -4942 +4947 @@
 !--     use terrain height array from file, if available. This flag is later used
 !--     to control initialization of surface attributes.
+!--     Todo: for the moment disable initialization of building roofs with
+!--     ground-floor-level properties. 
         surf_usm_h%ground_level = .FALSE. 
-        DO  m = 1, surf_usm_h%ns
-           i = surf_usm_h%i(m)
-           j = surf_usm_h%j(m)
-           k = surf_usm_h%k(m)
-!
-!--        Get local ground level. If no ground level is given in input file,
-!--        use default value.
-           ground_floor_level_l = ground_floor_level
-           IF ( building_pars_f%from_file )  THEN
-              IF ( building_pars_f%pars_xy(ind_gflh,j,i) /=                    &
-                   building_pars_f%fill )  &
-                 ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)          
-           ENDIF
-!
-!--        Determine height of surface element above ground level
-           IF (  terrain_height_f%from_file )  THEN
-              z_agl = zw(k) - terrain_height_f%var(j,i)
-           ELSE
-              z_agl = zw(k)
-           ENDIF
-!
-!--        Set flag for ground level
-           IF ( z_agl <= ground_floor_level_l )                                &
-              surf_usm_h%ground_level(m) = .TRUE.
-        ENDDO
+!         DO  m = 1, surf_usm_h%ns
+!            i = surf_usm_h%i(m)
+!            j = surf_usm_h%j(m)
+!            k = surf_usm_h%k(m)
+! !
+! !--        Get local ground level. If no ground level is given in input file,
+! !--        use default value.
+!            ground_floor_level_l = ground_floor_level
+!            IF ( building_pars_f%from_file )  THEN
+!               IF ( building_pars_f%pars_xy(ind_gflh,j,i) /=                    &
+!                    building_pars_f%fill )  &
+!                  ground_floor_level_l = building_pars_f%pars_xy(ind_gflh,j,i)          
+!            ENDIF
+! !
+! !--        Determine height of surface element above ground level
+!            IF (  terrain_height_f%from_file )  THEN
+!               z_agl = zw(k) - terrain_height_f%var(j,i)
+!            ELSE
+!               z_agl = zw(k)
+!            ENDIF
+! !
+! !--        Set flag for ground level
+!            IF ( z_agl <= ground_floor_level_l )                                &
+!               surf_usm_h%ground_level(m) = .TRUE.
+!         ENDDO
 
         DO  l = 0, 3