Ignore:
Timestamp:
Dec 21, 2011 5:48:03 PM (10 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/interaction_droplets_ptq.f90

    r484 r799  
    44! Current revisions:
    55! -----------------
    6 !
     6! Bugfix: pt_d_t(k) was missing in calculation of pt_p
    77!
    88! Former revisions:
     
    5353             DO  k = nzb_2d(j,i)+1, nzt
    5454                q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i)
    55                 pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i)
     55                pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k)
    5656             ENDDO
    5757          ENDDO
     
    8080       DO  k = nzb_2d(j,i)+1, nzt
    8181          q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i)
    82           pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i)
     82          pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k)
    8383       ENDDO
    8484
Note: See TracChangeset for help on using the changeset viewer.