Changeset 4671 for palm/trunk/SOURCE/init_3d_model.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/init_3d_model.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Implementation of downward facing USM and LSM surfaces
+!
+! 4648 2020-08-25 07:52:08Z raasch
 ! file re-formatted to follow the PALM coding standard
 !
@@ -1631 +1634 @@
        IF ( use_surface_fluxes  .AND.  constant_heatflux  .AND.  random_heatflux )  THEN
           IF ( surf_def_h(0)%ns >= 1 )  CALL disturb_heatflux( surf_def_h(0) )
-          IF ( surf_lsm_h%ns    >= 1 )  CALL disturb_heatflux( surf_lsm_h    )
-          IF ( surf_usm_h%ns    >= 1 )  CALL disturb_heatflux( surf_usm_h    )
+          IF ( surf_lsm_h(0)%ns >= 1 )  CALL disturb_heatflux( surf_lsm_h(0) )
+          IF ( surf_usm_h(0)%ns >= 1 )  CALL disturb_heatflux( surf_usm_h(0) )
        ENDIF
     ENDIF
@@ -1673 +1676 @@
                    mean_surface_level_height_l(sr) = mean_surface_level_height_l(sr) + zw(k-1)
                 ENDIF
-                IF ( surf_lsm_h%start_index(j,i) <= surf_lsm_h%end_index(j,i) )  THEN
-                   m = surf_lsm_h%start_index(j,i)
-                   k = surf_lsm_h%k(m)
+                IF ( surf_lsm_h(0)%start_index(j,i) <= surf_lsm_h(0)%end_index(j,i) )  THEN
+                   m = surf_lsm_h(0)%start_index(j,i)
+                   k = surf_lsm_h(0)%k(m)
                    mean_surface_level_height_l(sr) = mean_surface_level_height_l(sr) + zw(k-1)
                 ENDIF
-                IF ( surf_usm_h%start_index(j,i) <= surf_usm_h%end_index(j,i) )  THEN
-                   m = surf_usm_h%start_index(j,i)
-                   k = surf_usm_h%k(m)
+                IF ( surf_usm_h(0)%start_index(j,i) <= surf_usm_h(0)%end_index(j,i) )  THEN
+                   m = surf_usm_h(0)%start_index(j,i)
+                   k = surf_usm_h(0)%k(m)
                    mean_surface_level_height_l(sr) = mean_surface_level_height_l(sr) + zw(k-1)
                 ENDIF