Changeset 3901 for palm/trunk

Timestamp:

Apr 16, 2019 4:17:02 PM (6 years ago)

Author:

suehring

Message:

Set green-fraction on buildings to zero for now, in order to avoid crashes in the green-heat model

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Workaround - set green fraction to zero ( green-heat model crashes ).
+! 
+! 3896 2019-04-15 10:10:17Z suehring
+! 
+! 
+! 3896 2019-04-15 10:10:17Z suehring
 ! Bugfix, wrong index used for accessing building_pars from PIDS 
 ! 
@@ -5031 +5037 @@
           ENDDO
        ENDDO
+!
+!--    Just a work-around: Set green fraction on roof to zero (in favor of 
+!--    wall fraction). The green-heat model crashes due to some unknown reason. 
+!--    To be removed later.  
+       DO  m = 1, surf_usm_h%ns
+          IF ( surf_usm_h%frac(ind_pav_green,m) > 0.0_wp )  THEN
+             surf_usm_h%frac(ind_pav_green,m) = 0.0_wp
+             surf_usm_h%frac(ind_veg_wall,m)  = 1.0_wp
+          ENDIF
+       ENDDO
 
 !
@@ -5616 +5632 @@
 !> Green and substrate model as part of the urban surface model. The model predicts ground
 !> temperatures.
+!>
+!> Important: gree-heat model crashes due to unknown reason. Green fraction
+!> is thus set to zero (in favor of wall fraction). 
 !------------------------------------------------------------------------------!
     SUBROUTINE usm_green_heat_model
@@ -5648 +5667 @@
         !$OMP DO SCHEDULE (STATIC)
         DO  m = 1, surf_usm_h%ns
-
            IF (surf_usm_h%frac(ind_pav_green,m) > 0.0_wp) THEN
 !