Ignore:
Timestamp:
Nov 29, 2011 3:11:20 AM (12 years ago)
Author:
raasch
Message:

Bugfix for output of mean particle radius + preliminary works for implementing the Wang collision kernel

File:
1 edited

Legend:

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

    r772 r790  
    44! Current revisions:
    55! -----------------
     6! bugfix: calculation of 'pr' must depend on the particle weighting factor
    67!
    78! Former revisions:
     
    295296                      s_r4 = 0.0
    296297                      DO  n = psi, psi+prt_count(k,j,i)-1
    297                          s_r3 = s_r3 + particles(n)%radius**3
    298                          s_r4 = s_r4 + particles(n)%radius**4
     298                         s_r3 = s_r3 + particles(n)%radius**3 * &
     299                                       particles(n)%weight_factor
     300                         s_r4 = s_r4 + particles(n)%radius**4 * &
     301                                       particles(n)%weight_factor
    299302                      ENDDO
    300303                      IF ( s_r3 /= 0.0 )  THEN
Note: See TracChangeset for help on using the changeset viewer.