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

Timestamp:

Apr 17, 2020 10:12:45 AM (4 years ago)

Author:

suehring

Message:

Surface output: correct output of ground/wall-heat flux at USM surfaces; add conversion factor to heat- and momentum-flux outputs; data_output_2d: Unify output conversion of sensible and latent heat flux; data-output module: avoid uninitialized variables; restart_data_mpi_io: fix overlong lines

File:

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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Allocate array for wall heat flux, which is further used to aggregate tile
+! fractions in the surface output
+! 
+! 4495 2020-04-13 20:11:20Z raasch
 ! restart data handling with MPI-IO added
 ! 
@@ -1122 +1126 @@
         ENDDO
 !
-!--     allocate wall heat flux output array and set initial values. For horizontal surfaces
+!--     Allocate wall heat flux output arrays and set initial values. For horizontal surfaces
 !        ALLOCATE ( surf_usm_h%wshf(1:surf_usm_h%ns)    )  !can be removed 
         ALLOCATE ( surf_usm_h%wshf_eb(1:surf_usm_h%ns) )
@@ -1130 +1134 @@
         ALLOCATE ( surf_usm_h%iwghf_eb(1:surf_usm_h%ns) )
         ALLOCATE ( surf_usm_h%iwghf_eb_window(1:surf_usm_h%ns) )
+        ALLOCATE ( surf_usm_h%ghf(1:surf_usm_h%ns) )
         IF ( ALLOCATED( surf_usm_h%wshf    ) )  surf_usm_h%wshf    = 0.0_wp
         IF ( ALLOCATED( surf_usm_h%wshf_eb ) )  surf_usm_h%wshf_eb = 0.0_wp
@@ -1137 +1142 @@
         IF ( ALLOCATED( surf_usm_h%iwghf_eb ) )  surf_usm_h%iwghf_eb = 0.0_wp
         IF ( ALLOCATED( surf_usm_h%iwghf_eb_window ) )  surf_usm_h%iwghf_eb_window = 0.0_wp
+        IF ( ALLOCATED( surf_usm_h%ghf ) )  surf_usm_h%ghf = 0.0_wp
 !
 !--     Now, for vertical surfaces
@@ -1147 +1153 @@
            ALLOCATE ( surf_usm_v(l)%iwghf_eb(1:surf_usm_v(l)%ns) )
            ALLOCATE ( surf_usm_v(l)%iwghf_eb_window(1:surf_usm_v(l)%ns) )
+           ALLOCATE ( surf_usm_v(l)%ghf(1:surf_usm_v(l)%ns) )
            IF ( ALLOCATED( surf_usm_v(l)%wshf    ) )  surf_usm_v(l)%wshf    = 0.0_wp
            IF ( ALLOCATED( surf_usm_v(l)%wshf_eb ) )  surf_usm_v(l)%wshf_eb = 0.0_wp
@@ -1154 +1161 @@
            IF ( ALLOCATED( surf_usm_v(l)%iwghf_eb ) )  surf_usm_v(l)%iwghf_eb = 0.0_wp
            IF ( ALLOCATED( surf_usm_v(l)%iwghf_eb_window ) )  surf_usm_v(l)%iwghf_eb_window = 0.0_wp
+           IF ( ALLOCATED( surf_usm_v(l)%ghf ) )  surf_usm_v(l)%ghf = 0.0_wp
         ENDDO
 !