Changeset 1738 for palm/trunk/SOURCE/calc_mean_profile.f90

Timestamp:

Dec 18, 2015 1:56:05 PM (9 years ago)

Author:

raasch

Message:

bugfixes for calculations in statistical regions which do not contain grid points in the lowest vertical levels, mean surface level height considered in the calculation of the characteristic vertical velocity

File:

• palm/trunk/SOURCE/calc_mean_profile.f90 (modified) (2 diffs)

palm/trunk/SOURCE/calc_mean_profile.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! bugfix: if a layer is completely filled with topography, no mean is calculated
 ! 
 ! Former revisions:
@@ -126 +126 @@
 #endif
 
-          hom(:,1,pr,0) = sums(:,pr) / ngp_2dh_s_inner(:,0)
+          DO  k = nzb, nzt+1
+             IF ( ngp_2dh_s_inner(k,0) /= 0 )  THEN
+                hom(k,1,pr,0) = sums(k,pr) / ngp_2dh_s_inner(k,0)
+             ENDIF
+          ENDDO
 
        ENDIF