Changeset 4750 for palm/trunk/SOURCE/surface_layer_fluxes_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/surface_layer_fluxes_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Revision of 10-cm temperature at vertical walls - assume grid-cell temperature rather than employ
+! MOST
+!
+! 4717 2020-09-30 22:27:40Z pavelkrc
 ! Fixes and optimizations of OpenMP parallelization, formatting of OpenMP
 ! directives (J. Resler)
@@ -2044 +2048 @@
 
        CASE ( '10cm' )
-
-          DO  m = 1, surf%ns
-             surf%pt_10cm(m) = surf%pt_surface(m) + surf%ts(m) / kappa                             &
+!
+!--       For horizontal upward-facing surfaces 10-cm temperature can be calculated using MOST.
+          IF ( .NOT.  downward  .AND.  .NOT. surf_vertical )  THEN
+             DO  m = 1, surf%ns
+                surf%pt_10cm(m) = surf%pt_surface(m) + surf%ts(m) / kappa                          &
                                * ( LOG( 0.1_wp /  surf%z0h(m) ) - psi_h( 0.1_wp / surf%ol(m) )     &
                                    + psi_h( surf%z0h(m) / surf%ol(m) ) )
-          ENDDO
+             ENDDO
+!
+!--       At vertical surfaces 10-cm temperature cannot be calculated via MOST as the Obukhov length
+!--       and temperature scaling parameter are not calculated. Hence, set 10-cm temperature to
+!--       the grid-cell temperature.
+          ELSE
+             DO  m = 1, surf%ns
+                surf%pt_10cm(m) = pt(surf%k(m)+surf%koff,surf%j(m)+surf%joff,surf%i(m)+surf%ioff)
+             ENDDO
+          ENDIF
 
     END SELECT