Ignore:
Timestamp:
Dec 21, 2011 5:48:03 PM (9 years ago)
Author:
franke
Message:

Implementation of Wang collision kernel and bugfix for calculation of vpt, pt_p, and ec in case of cloud droplets

File:
1 edited

Legend:

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

    r484 r799  
    44! Current revisions:
    55! -----------------
    6 !
     6! Bugfix: ql is now included in calculation of vpt in case of
     7!         cloud droplets
    78!
    89! Former revisions:
     
    3233    INTEGER :: k
    3334
    34     IF ( .NOT. cloud_physics ) THEN
     35    IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN
    3536       vpt = pt * ( 1.0 + 0.61 * q )
    36     ELSE
     37    ELSE IF (cloud_physics) THEN
    3738       DO  k = nzb, nzt+1
    3839          vpt(k,:,:) = ( pt(k,:,:) + pt_d_t(k) * l_d_cp * ql(k,:,:) ) * &
    3940                       ( 1.0 + 0.61 * q(k,:,:) - 1.61 * ql(k,:,:) )
    4041       ENDDO
     42    ELSE
     43       vpt = pt * ( 1.0 + 0.61 * q - 1.61 * ql )
    4144    ENDIF
    4245
Note: See TracChangeset for help on using the changeset viewer.