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

Timestamp:

Sep 24, 2018 3:42:55 PM (6 years ago)

Author:

knoop

Message:

Modularization of all bulk cloud physics code components

File:

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

palm/trunk/SOURCE/interaction_droplets_ptq.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Modularization of all bulk cloud physics code components
+! 
+! 3040 2018-05-25 10:22:08Z schwenkel
 ! Changed the name specific humidity to mixing ratio
 ! 
@@ -74 +77 @@
 !
 ! 799 2011-12-21 17:48:03Z franke
-! Bugfix: pt_d_t(k) was missing in calculation of pt_p
+! Bugfix: d_exner(k) was missing in calculation of pt_p
 !
 ! RCS Log replace by Id keyword, revision history cleaned up
@@ -89 +92 @@
  MODULE interaction_droplets_ptq_mod
  
+
+    USE arrays_3d,                                                             &
+        ONLY:  pt_p, ql_c, q_p, d_exner
+
+    USE basic_constants_and_equations_mod,                                     &
+        ONLY:  lv_d_cp
+
+    USE indices,                                                               &
+        ONLY:  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
+
+    USE kinds
+
+    USE pegrid
 
     PRIVATE
@@ -111 +127 @@
     SUBROUTINE i_droplets_ptq
 
-       USE arrays_3d,                                                          &
-           ONLY:  pt_p, ql_c, q_p
-           
-       USE cloud_parameters,                                                   &
-           ONLY:  l_d_cp, pt_d_t
-           
-       USE indices,                                                            &
-           ONLY:  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
-           
-       USE kinds
-
-       USE pegrid
-
        IMPLICIT NONE
 
@@ -140 +143 @@
 
                 q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i) * flag
-                pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k)   &
+                pt_p(k,j,i) = pt_p(k,j,i) + lv_d_cp * ql_c(k,j,i) * d_exner(k) &
                                                         * flag
              ENDDO
@@ -155 +158 @@
 !------------------------------------------------------------------------------!
     SUBROUTINE i_droplets_ptq_ij( i, j )
-
-       USE arrays_3d,                                                          &
-           ONLY:  pt_p, ql_c, q_p
-
-       USE cloud_parameters,                                                   &
-           ONLY:  l_d_cp, pt_d_t
-
-       USE indices,                                                            &
-           ONLY:  nzb, nzt, wall_flags_0
-
-       USE kinds,                                                              &
-           ONLY:  iwp, wp
-
-       USE pegrid
 
        IMPLICIT NONE
@@ -185 +174 @@
 
           q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i) * flag
-          pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k) * flag
+          pt_p(k,j,i) = pt_p(k,j,i) + lv_d_cp * ql_c(k,j,i) * d_exner(k) * flag
        ENDDO