Ignore:
Timestamp:
May 30, 2017 5:47:52 PM (4 years ago)
Author:
suehring
Message:

Adjustments according new topography and surface-modelling concept implemented

File:
1 edited

Legend:

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

    r2101 r2232  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Ajdustments for new topography masking
    2323!
    2424! Former revisions:
     
    111111!       DO  i = nxl, nxr
    112112!          DO  j = nys, nyn
    113 !             DO  k = nzb_s_inner(j,i)+1, nzt
     113!             DO  k = nzb+1, nzt
    114114!!
    115115!!--             Sample on how to calculate the profile of the resolved-scale
    116116!!--             horizontal momentum flux u*v*
    117 !                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +               &
    118 !                      ( 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) *     &
    119 !                      ( 0.5_wp * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) )       &
    120 !                                                 * rmask(j,i,sr)
     117!                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +             &
     118!                      ( 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) *&
     119!                      ( 0.5_wp * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) )  &
     120!                                     * rmask(j,i,sr)                          &
     121!                                     * MERGE( 1.0_wp, 0.0_wp,                 &
     122!                                              BTEST( wall_flags_0(k,j,i), 0 ) )
    121123!!
    122124!!--             Further profiles can be defined and calculated by increasing
Note: See TracChangeset for help on using the changeset viewer.