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

Timestamp:

Oct 16, 2020 2:27:48 PM (4 years ago)

Author:

suehring

Message:

Indoor-model new: Namelist parameter added to switch-off/on the indoor model during wall/soil spinup; bugfixes concerning indoor model: bugfix in window-wall treatment during spinup - in the urban-surface model the window fraction is set to zero during spinup, so it is done here also; bugfix in wall treatment - inner wall temperature was too low due to wrong weighting of wall/green/window fractions; Revision of 10-cm temperature at vertical walls - assume grid-cell temperature rather than employ MOST; call hourly-based indoor model only once per hour during spinup, not every timestep; add missing dependency in Makefile; urban-surface model: bugfix in openmp directive

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! - bugfix in openmp directive
+! - make t_green_h and t_green_v public (required in indoor model)
+! 
+! 4747 2020-10-16 09:19:57Z pavelkrc
 ! Fix window absorptivity calculation (correctly account for 2-sided reflection)
 ! 
@@ -880 +884 @@
            nzb_wall,                                                                               &
            nzt_wall,                                                                               &
+           t_green_h,                                                                              &
+           t_green_v,                                                                              &
            t_wall_h,                                                                               &
            t_wall_v,                                                                               &
@@ -5149 +5155 @@
     IMPLICIT NONE
 
+    LOGICAL       ::  during_spinup  !< if true, no calculation of window temperatures
     LOGICAL       ::  horizontal     !< Flag indicating horizontal or vertical surfaces
+
+    INTEGER(iwp)  ::  kw             !< grid index - wall depth
     INTEGER(iwp)  ::  l              !< direction index
-    LOGICAL       ::  during_spinup  !< if true, no calculation of window temperatures
-
-    INTEGER(iwp)  ::  i,j,k,kw, m    !< running indices
+    INTEGER(iwp)  ::  m              !< running index for surface elements
 
     REAL(wp)  ::  win_absorp        !< absorption coefficient from transmissivity
@@ -5197 +5204 @@
 !
 !-- Cycle for all surfaces in given direction
-    !$OMP PARALLEL DO PRIVATE (m, i, j, k, kw, wtend, wintend, win_absorp, wall_mod) SCHEDULE (STATIC)
+    !$OMP PARALLEL DO PRIVATE (m, kw, wtend, wintend, win_absorp) SCHEDULE (STATIC)
     DO  m = 1, surf%ns
-!
-!--    Obtain indices
-       i = surf%i(m)
-       j = surf%j(m)
-       k = surf%k(m)
 !
 !--    Prognostic equation for ground/roof temperature t_wall