Changeset 2968 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Apr 13, 2018 11:52:24 AM (6 years ago)

Author:

suehring

Message:

Bugfixes in initalization of land-surface model as well as timeseries data output in case of elevated model surfaces.

File:

• palm/trunk/SOURCE/flow_statistics.f90 (modified) (3 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in output of timeseries quantities in case of elevated model surfaces.
+! 
+! 2817 2018-02-19 16:32:21Z knoop
 ! Preliminary gust module interface implemented
 ! 
@@ -292 +295 @@
     USE indices,                                                               &
         ONLY:   ngp_2dh, ngp_2dh_s_inner, ngp_3d, ngp_3d_inner, ngp_sums,      &
-                ngp_sums_ls, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
+                ngp_sums_ls, nxl, nxr, nyn, nys, nzb, nzt, topo_min_level,     &
+                wall_flags_0
         
     USE kinds
@@ -2102 +2106 @@
 
 !
-!--    Collect the time series quantities
-       ts_value(1,sr) = hom(nzb+4,1,pr_palm,sr)     ! E
-       ts_value(2,sr) = hom(nzb+5,1,pr_palm,sr)     ! E*
+!--    Collect the time series quantities. Please note, timeseries quantities 
+!--    which are collected from horizontally averaged profiles, e.g. wpt 
+!--    or pt(zp), are treated specially. In case of elevated model surfaces, 
+!--    index nzb+1 might be within topography and data will be zero. Therefore,
+!--    take value for the first atmosphere index, which is topo_min_level+1. 
+       ts_value(1,sr) = hom(nzb+4,1,pr_palm,sr)        ! E
+       ts_value(2,sr) = hom(nzb+5,1,pr_palm,sr)        ! E*
        ts_value(3,sr) = dt_3d
-       ts_value(4,sr) = hom(nzb,1,pr_palm,sr)       ! u*
-       ts_value(5,sr) = hom(nzb+3,1,pr_palm,sr)     ! th*
+       ts_value(4,sr) = hom(nzb,1,pr_palm,sr)          ! u*
+       ts_value(5,sr) = hom(nzb+3,1,pr_palm,sr)        ! th*
        ts_value(6,sr) = u_max
        ts_value(7,sr) = v_max
        ts_value(8,sr) = w_max
-       ts_value(9,sr) = hom(nzb+10,1,pr_palm,sr)    ! new divergence
-       ts_value(10,sr) = hom(nzb+9,1,pr_palm,sr)    ! old Divergence
-       ts_value(11,sr) = hom(nzb+6,1,pr_palm,sr)    ! z_i(1)
-       ts_value(12,sr) = hom(nzb+7,1,pr_palm,sr)    ! z_i(2)
-       ts_value(13,sr) = hom(nzb+8,1,pr_palm,sr)    ! w*
-       ts_value(14,sr) = hom(nzb,1,16,sr)           ! w'pt'   at k=0
-       ts_value(15,sr) = hom(nzb+1,1,16,sr)         ! w'pt'   at k=1
-       ts_value(16,sr) = hom(nzb+1,1,18,sr)         ! wpt     at k=1
-       ts_value(17,sr) = hom(nzb+14,1,pr_palm,sr)   ! pt(0)
-       ts_value(18,sr) = hom(nzb+1,1,4,sr)          ! pt(zp)
-       ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr)    ! u'w'    at k=0
-       ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)    ! v'w'    at k=0
-       ts_value(21,sr) = hom(nzb,1,48,sr)           ! w"q"    at k=0
+       ts_value(9,sr) = hom(nzb+10,1,pr_palm,sr)       ! new divergence
+       ts_value(10,sr) = hom(nzb+9,1,pr_palm,sr)       ! old Divergence
+       ts_value(11,sr) = hom(nzb+6,1,pr_palm,sr)       ! z_i(1)
+       ts_value(12,sr) = hom(nzb+7,1,pr_palm,sr)       ! z_i(2)
+       ts_value(13,sr) = hom(nzb+8,1,pr_palm,sr)       ! w*
+       ts_value(14,sr) = hom(nzb,1,16,sr)              ! w'pt'   at k=0
+       ts_value(15,sr) = hom(topo_min_level+1,1,16,sr) ! w'pt'   at k=1
+       ts_value(16,sr) = hom(topo_min_level+1,1,18,sr) ! wpt     at k=1
+       ts_value(17,sr) = hom(nzb+14,1,pr_palm,sr)      ! pt(0)
+       ts_value(18,sr) = hom(topo_min_level+1,1,4,sr)  ! pt(zp)
+       ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr)       ! u'w'    at k=0
+       ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)       ! v'w'    at k=0
+       ts_value(21,sr) = hom(nzb,1,48,sr)              ! w"q"    at k=0
 
        IF ( .NOT. neutral )  THEN