Ignore:
Timestamp:
Apr 13, 2018 11:52:24 AM (3 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/flow_statistics.f90

    r2817 r2968  
    2525! -----------------
    2626! $Id$
     27! Bugfix in output of timeseries quantities in case of elevated model surfaces.
     28!
     29! 2817 2018-02-19 16:32:21Z knoop
    2730! Preliminary gust module interface implemented
    2831!
     
    292295    USE indices,                                                               &
    293296        ONLY:   ngp_2dh, ngp_2dh_s_inner, ngp_3d, ngp_3d_inner, ngp_sums,      &
    294                 ngp_sums_ls, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
     297                ngp_sums_ls, nxl, nxr, nyn, nys, nzb, nzt, topo_min_level,     &
     298                wall_flags_0
    295299       
    296300    USE kinds
     
    21022106
    21032107!
    2104 !--    Collect the time series quantities
    2105        ts_value(1,sr) = hom(nzb+4,1,pr_palm,sr)     ! E
    2106        ts_value(2,sr) = hom(nzb+5,1,pr_palm,sr)     ! E*
     2108!--    Collect the time series quantities. Please note, timeseries quantities
     2109!--    which are collected from horizontally averaged profiles, e.g. wpt
     2110!--    or pt(zp), are treated specially. In case of elevated model surfaces,
     2111!--    index nzb+1 might be within topography and data will be zero. Therefore,
     2112!--    take value for the first atmosphere index, which is topo_min_level+1.
     2113       ts_value(1,sr) = hom(nzb+4,1,pr_palm,sr)        ! E
     2114       ts_value(2,sr) = hom(nzb+5,1,pr_palm,sr)        ! E*
    21072115       ts_value(3,sr) = dt_3d
    2108        ts_value(4,sr) = hom(nzb,1,pr_palm,sr)       ! u*
    2109        ts_value(5,sr) = hom(nzb+3,1,pr_palm,sr)     ! th*
     2116       ts_value(4,sr) = hom(nzb,1,pr_palm,sr)          ! u*
     2117       ts_value(5,sr) = hom(nzb+3,1,pr_palm,sr)        ! th*
    21102118       ts_value(6,sr) = u_max
    21112119       ts_value(7,sr) = v_max
    21122120       ts_value(8,sr) = w_max
    2113        ts_value(9,sr) = hom(nzb+10,1,pr_palm,sr)    ! new divergence
    2114        ts_value(10,sr) = hom(nzb+9,1,pr_palm,sr)    ! old Divergence
    2115        ts_value(11,sr) = hom(nzb+6,1,pr_palm,sr)    ! z_i(1)
    2116        ts_value(12,sr) = hom(nzb+7,1,pr_palm,sr)    ! z_i(2)
    2117        ts_value(13,sr) = hom(nzb+8,1,pr_palm,sr)    ! w*
    2118        ts_value(14,sr) = hom(nzb,1,16,sr)           ! w'pt'   at k=0
    2119        ts_value(15,sr) = hom(nzb+1,1,16,sr)        ! w'pt'   at k=1
    2120        ts_value(16,sr) = hom(nzb+1,1,18,sr)        ! wpt     at k=1
    2121        ts_value(17,sr) = hom(nzb+14,1,pr_palm,sr)   ! pt(0)
    2122        ts_value(18,sr) = hom(nzb+1,1,4,sr)          ! pt(zp)
    2123        ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr)    ! u'w'    at k=0
    2124        ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)    ! v'w'    at k=0
    2125        ts_value(21,sr) = hom(nzb,1,48,sr)           ! w"q"    at k=0
     2121       ts_value(9,sr) = hom(nzb+10,1,pr_palm,sr)       ! new divergence
     2122       ts_value(10,sr) = hom(nzb+9,1,pr_palm,sr)       ! old Divergence
     2123       ts_value(11,sr) = hom(nzb+6,1,pr_palm,sr)       ! z_i(1)
     2124       ts_value(12,sr) = hom(nzb+7,1,pr_palm,sr)       ! z_i(2)
     2125       ts_value(13,sr) = hom(nzb+8,1,pr_palm,sr)       ! w*
     2126       ts_value(14,sr) = hom(nzb,1,16,sr)              ! w'pt'   at k=0
     2127       ts_value(15,sr) = hom(topo_min_level+1,1,16,sr) ! w'pt'   at k=1
     2128       ts_value(16,sr) = hom(topo_min_level+1,1,18,sr) ! wpt     at k=1
     2129       ts_value(17,sr) = hom(nzb+14,1,pr_palm,sr)      ! pt(0)
     2130       ts_value(18,sr) = hom(topo_min_level+1,1,4,sr)  ! pt(zp)
     2131       ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr)       ! u'w'    at k=0
     2132       ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)       ! v'w'    at k=0
     2133       ts_value(21,sr) = hom(nzb,1,48,sr)              ! w"q"    at k=0
    21262134
    21272135       IF ( .NOT. neutral )  THEN
Note: See TracChangeset for help on using the changeset viewer.