Ignore:
Timestamp:
Nov 29, 2011 3:11:20 AM (10 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/data_output_2d.f90

    r772 r790  
    44! Current revisions:
    55! -----------------
     6! bugfix: calculation of 'pr' must depend on the particle weighting factor
    67!
    78! Former revisions:
     
    303304                               s_r4 = 0.0
    304305                               DO  n = psi, psi+prt_count(k,j,i)-1
    305                                   s_r3 = s_r3 + particles(n)%radius**3
    306                                   s_r4 = s_r4 + particles(n)%radius**4
     306                                  s_r3 = s_r3 + particles(n)%radius**3 * &
     307                                                particles(n)%weight_factor
     308                                  s_r4 = s_r4 + particles(n)%radius**4 * &
     309                                                particles(n)%weight_factor
    307310                               ENDDO
    308311                               IF ( s_r3 /= 0.0 )  THEN
Note: See TracChangeset for help on using the changeset viewer.