Ignore:
Timestamp:
Feb 2, 2011 2:31:56 PM (10 years ago)
Author:
raasch
Message:

New:
---

further adjustments on Tsubame and concerning openMP usage
(mrun, mbuild, subjob)

Changed:


Errors:


Bugfix in calculation of divergence of vertical flux of resolved scale
energy, pressure fluctuations, and flux of pressure fluctuation itself
(flow statistics)

Bugfix: module pegrid was missing. (user_statistics)

File:
1 edited

Legend:

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

    r674 r678  
    44! Current revisions:
    55! -----------------
     6! Bugfix in calculation of divergence of vertical flux of resolved scale
     7! energy, pressure fluctuations, and flux of pressure fluctuation itself
    68!
    79! Former revisions:
     
    802804                   sums_ll(k,1) = sums_ll(k,1) + 0.5 * w(k,j,i) * (        &
    803805                  ( 0.25 * ( u(k,j,i)+u(k+1,j,i)+u(k,j,i+1)+u(k+1,j,i+1)   &
    804                            - 2.0 * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) )     &
     806                           - 0.5 * ( hom(k,1,1,sr) + hom(k+1,1,1,sr) )     &
    805807                           ) )**2                                          &
    806808                + ( 0.25 * ( v(k,j,i)+v(k+1,j,i)+v(k,j+1,i)+v(k+1,j+1,i)   &
    807                            - 2.0 * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) )     &
     809                           - 0.5 * ( hom(k,1,2,sr) + hom(k+1,1,2,sr) )     &
    808810                           ) )**2                                          &
    809811                   + w(k,j,i)**2                                  )
     
    820822          sums_ll(nzt+1,2) = 0.0
    821823
    822           DO  k = nzb_s_inner(j,i)+1, nzt
     824          DO  k = nzb+1, nzt
    823825             sums_l(k,55,tn) = ( sums_ll(k,1) - sums_ll(k-1,1) ) * ddzw(k)
    824826             sums_l(k,56,tn) = ( sums_ll(k,2) - sums_ll(k-1,2) ) * ddzw(k)
Note: See TracChangeset for help on using the changeset viewer.