Changeset 799 for palm/trunk/SOURCE/interaction_droplets_ptq.f90

Timestamp:

Dec 21, 2011 5:48:03 PM (12 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:

• palm/trunk/SOURCE/interaction_droplets_ptq.f90 (modified) (3 diffs)

palm/trunk/SOURCE/interaction_droplets_ptq.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Bugfix: pt_d_t(k) was missing in calculation of pt_p
 !
 ! Former revisions:
@@ -53 +53 @@
              DO  k = nzb_2d(j,i)+1, nzt
                 q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i)
-                pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i)
+                pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k)
              ENDDO
           ENDDO
@@ -80 +80 @@
        DO  k = nzb_2d(j,i)+1, nzt
           q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i)
-          pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i)
+          pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k)
        ENDDO