Changeset 1776


Ignore:
Timestamp:
Mar 2, 2016 5:54:58 PM (6 years ago)
Author:
hoffmann
Message:

Bugfix in computation of collection efficiencies

File:
1 edited

Legend:

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

    r1683 r1776  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Bugfix: Collection efficiencies must be calculated for the larger droplet.
    2222!
    2323! Former revisions:
     
    875875             IF ( ir < 16 )  THEN
    876876                IF ( ir >= 2 )  THEN
    877                    pp = ( ( radclass(j) * 1.0E06_wp ) - r0(ir-1) ) / &
     877                   pp = ( ( MAX(radclass(j),radclass(i)) * 1.0E06_wp ) - r0(ir-1) ) / &
    878878                        ( r0(ir) - r0(ir-1) )
    879879                   qq = ( rq - rat(iq-1) ) / ( rat(iq) - rat(iq-1) )
     
    10411041             IF ( ir < 8 )  THEN
    10421042                IF ( ir >= 2 )  THEN
    1043                    pp = ( radclass(j)*1.0E6_wp - r0(ir-1) ) / ( r0(ir) - r0(ir-1) )
     1043                   pp = ( MAX(radclass(j),radclass(i))*1.0E6_wp - r0(ir-1) ) / &
     1044                        ( r0(ir) - r0(ir-1) )
    10441045                   qq = ( rq - rat(iq-1) ) / ( rat(iq) - rat(iq-1) )
    10451046                   y2 = ( 1.0_wp - pp ) * ( 1.0_wp - qq ) * ecoll_100(ir-1,iq-1) + &
Note: See TracChangeset for help on using the changeset viewer.