Changeset 4671 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

Timestamp:

Sep 9, 2020 8:27:58 PM (4 years ago)

Author:

pavelkrc

Message:

Radiative transfer model RTM version 4.1

File:

• palm/trunk/SOURCE/surface_layer_fluxes_mod.f90 (modified) (8 diffs)

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implementation of downward facing USM and LSM surfaces
+!
+! 4594 2020-07-09 15:01:00Z suehring
 ! Include k index in OMP PRIVATE statements
 !
@@ -285 +288 @@
           CALL calc_ln
        ENDIF
+       IF ( surf_lsm_h(l)%ns >= 1 )  THEN
+          surf => surf_lsm_h(l)
+          CALL calc_ln
+       ENDIF
+       IF ( surf_usm_h(l)%ns >= 1 )  THEN
+          surf => surf_usm_h(l)
+          CALL calc_ln
+       ENDIF
     ENDDO
-    IF ( surf_lsm_h%ns >= 1 )  THEN
-       surf => surf_lsm_h
-       CALL calc_ln
-    ENDIF
-    IF ( surf_usm_h%ns >= 1 )  THEN
-       surf => surf_usm_h
-       CALL calc_ln
-    ENDIF
 
     DO  l = 0, 3
@@ -311 +314 @@
 !
 !-- Derive potential temperature and specific humidity at first grid level from the fields pt and q
-!-- First call for horizontal default-type surfaces (l=0 - upward facing, l=1 - downward facing)
     DO  l = 0, 1
+!--    First call for horizontal default-type surfaces (l=0 - upward facing, l=1 - downward facing)
        IF ( surf_def_h(l)%ns >= 1  )  THEN
           surf => surf_def_h(l)
@@ -324 +327 @@
           ENDIF
        ENDIF
+!
+!--    Call for natural-type horizontal surfaces
+       IF ( surf_lsm_h(l)%ns >= 1  )  THEN
+          surf => surf_lsm_h(l)
+          CALL calc_pt_q
+       ENDIF
+!
+!--    Call for urban-type horizontal surfaces
+       IF ( surf_usm_h(l)%ns >= 1  )  THEN
+          surf => surf_usm_h(l)
+          CALL calc_pt_q
+       ENDIF
     ENDDO
-!
-!-- Call for natural-type horizontal surfaces
-    IF ( surf_lsm_h%ns >= 1  )  THEN
-       surf => surf_lsm_h
-       CALL calc_pt_q
-    ENDIF
-
-!
-!-- Call for urban-type horizontal surfaces
-    IF ( surf_usm_h%ns >= 1  )  THEN
-       surf => surf_usm_h
-       CALL calc_pt_q
-    ENDIF
-
 !
 !-- Call for natural-type vertical surfaces
@@ -373 +374 @@
 !
 !-- Natural-type horizontal surfaces
-    IF ( surf_lsm_h%ns >= 1 )  THEN
-       surf => surf_lsm_h
+    IF ( surf_lsm_h(0)%ns >= 1 )  THEN
+       surf => surf_lsm_h(0)
        CALL calc_uvw_abs
        IF ( .NOT. neutral )  CALL calc_ol
@@ -383 +384 @@
 !
 !-- Urban-type horizontal surfaces
-    IF ( surf_usm_h%ns >= 1 )  THEN
-       surf => surf_usm_h
+    IF ( surf_usm_h(0)%ns >= 1 )  THEN
+       surf => surf_usm_h(0)
        CALL calc_uvw_abs
        IF ( .NOT. neutral )  CALL calc_ol
@@ -396 +397 @@
 
 !
-!-- Treat downward-facing horizontal surfaces. Note, so far, these are always default type.
+!-- Treat downward-facing horizontal surfaces.
 !-- Stratification is not considered in this case, hence, no further distinction between different
 !-- most_method is required.
+    downward = .TRUE.
+!-- Default type.
     IF ( surf_def_h(1)%ns >= 1 )  THEN
-       downward = .TRUE.
        surf => surf_def_h(1)
        CALL calc_uvw_abs
        CALL calc_us
        CALL calc_surface_fluxes
-       downward = .FALSE.
-    ENDIF
+    ENDIF
+!-- Natural surface type.
+    IF ( surf_lsm_h(1)%ns >= 1 )  THEN
+       surf => surf_lsm_h(1)
+       CALL calc_uvw_abs
+       CALL calc_us
+       CALL calc_surface_fluxes
+    ENDIF
+!-- Urban surface type.
+    IF ( surf_usm_h(1)%ns >= 1 )  THEN
+       surf => surf_usm_h(1)
+       CALL calc_uvw_abs
+       CALL calc_us
+       CALL calc_surface_fluxes
+    ENDIF
+    downward = .FALSE.
 !
 !-- Calculate surfaces fluxes at vertical surfaces for momentum and subgrid-scale TKE. No stability
@@ -664 +680 @@
 !-- In case of runs with neutral statification, set Obukhov length to a large value
     IF ( neutral )  THEN
-       IF ( surf_def_h(0)%ns >= 1 )  surf_def_h(0)%ol = 1.0E10_wp
-       IF ( surf_lsm_h%ns    >= 1 )  surf_lsm_h%ol    = 1.0E10_wp
-       IF ( surf_usm_h%ns    >= 1 )  surf_usm_h%ol    = 1.0E10_wp
-
+       DO  l = 0, 1
+          IF ( surf_def_h(l)%ns >= 1  .AND.                                                        &
+               ALLOCATED( surf_def_h(l)%ol ) )  surf_def_h(l)%ol = 1.0E10_wp
+          IF ( surf_lsm_h(l)%ns    >= 1  .AND.                                                     &
+               ALLOCATED( surf_lsm_h(l)%ol ) )  surf_lsm_h(l)%ol = 1.0E10_wp
+          IF ( surf_usm_h(l)%ns    >= 1  .AND.                                                     &
+               ALLOCATED( surf_usm_h(l)%ol ) )  surf_usm_h(l)%ol = 1.0E10_wp
+       ENDDO
        DO  l = 0, 3
           IF ( surf_def_v(l)%ns >= 1  .AND.                                                        &