Ignore:
Timestamp:
Sep 17, 2015 8:12:24 AM (6 years ago)
Author:
raasch
Message:

bugfix in calculation of vertical profile of energy production by turbulent transport of TKE
bugfix for output of configuration settings (mbuild)
output of login init command (mbuild)

File:
1 edited

Legend:

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

    r1594 r1652  
    2121! Current revisions:
    2222! -----------------
    23 !
     23! bugfix in calculation of energy production by turbulent transport of TKE
    2424!
    2525! Former revisions:
     
    993993
    994994                   sums_ll(k,1) = sums_ll(k,1) + 0.5_wp * w(k,j,i) * (         &
    995                   ( 0.25_wp * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1)    &
    996                             - 0.5_wp * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) )     &
    997                               ) )**2                                           &
    998                 + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i)    &
    999                             - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) )     &
    1000                               ) )**2                                           &
     995                  ( 0.25_wp * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1) )  &
     996                            - 0.5_wp * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) ) )**2&
     997                + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i) )  &
     998                            - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) ) )**2&                                           &
    1001999                + w(k,j,i)**2                                        )
    10021000
     
    28672865
    28682866                   sums_ll(k,1) = sums_ll(k,1) + 0.5_wp * w(k,j,i) * (         &
    2869                   ( 0.25_wp * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1)    &
    2870                               - 0.5_wp * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) )   &
    2871                               ) )**2                                           &
    2872                 + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i)    &
    2873                               - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) )   &
    2874                               ) )**2                                           &
     2867                  ( 0.25_wp * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1) )  &
     2868                            - 0.5_wp * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) ) )**2&                                           &
     2869                + ( 0.25_wp * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i) )  &
     2870                            - 0.5_wp * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) ) )**2&                                          &
    28752871                + w(k,j,i)**2                                        )
    28762872
Note: See TracChangeset for help on using the changeset viewer.