source: palm/trunk/SOURCE/interaction_droplets_ptq.f90 @ 829

Last change on this file since 829 was 800, checked in by franke, 12 years ago

last commit documented

  • Property svn:keywords set to Id
File size: 2.2 KB
Line 
1 MODULE interaction_droplets_ptq_mod
2
3!------------------------------------------------------------------------------!
4! Current revisions:
5! -----------------
6!
7!
8! Former revisions:
9! -----------------
10! $Id: interaction_droplets_ptq.f90 800 2011-12-21 18:13:51Z raasch $
11!
12! 799 2011-12-21 17:48:03Z franke
13! Bugfix: pt_d_t(k) was missing in calculation of pt_p
14!
15! RCS Log replace by Id keyword, revision history cleaned up
16!
17! Revision 1.1  2005/06/26 19:57:47  raasch
18! Initial revision
19!
20!
21! Description:
22! ------------
23! Release of latent heat and change of specific humidity due to condensation /
24! evaporation of droplets.
25!------------------------------------------------------------------------------!
26
27    PRIVATE
28    PUBLIC interaction_droplets_ptq
29
30    INTERFACE interaction_droplets_ptq
31       MODULE PROCEDURE interaction_droplets_ptq
32       MODULE PROCEDURE interaction_droplets_ptq_ij
33    END INTERFACE interaction_droplets_ptq
34 
35 CONTAINS
36
37
38!------------------------------------------------------------------------------!
39! Call for all grid points
40!------------------------------------------------------------------------------!
41    SUBROUTINE interaction_droplets_ptq
42
43       USE arrays_3d
44       USE cloud_parameters
45       USE control_parameters
46       USE indices
47
48       USE pegrid
49
50       IMPLICIT NONE
51
52       INTEGER ::  i, j, k
53
54 
55       DO  i = nxl, nxr
56          DO  j = nys, nyn
57             DO  k = nzb_2d(j,i)+1, nzt
58                q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i)
59                pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k)
60             ENDDO
61          ENDDO
62       ENDDO
63
64    END SUBROUTINE interaction_droplets_ptq
65
66
67!------------------------------------------------------------------------------!
68! Call for grid point i,j
69!------------------------------------------------------------------------------!
70    SUBROUTINE interaction_droplets_ptq_ij( i, j )
71
72       USE arrays_3d
73       USE cloud_parameters
74       USE control_parameters
75       USE indices
76
77       USE pegrid
78
79       IMPLICIT NONE
80
81       INTEGER ::  i, j, k
82
83
84       DO  k = nzb_2d(j,i)+1, nzt
85          q_p(k,j,i)  = q_p(k,j,i)  - ql_c(k,j,i)
86          pt_p(k,j,i) = pt_p(k,j,i) + l_d_cp * ql_c(k,j,i) * pt_d_t(k)
87       ENDDO
88
89    END SUBROUTINE interaction_droplets_ptq_ij
90
91 END MODULE interaction_droplets_ptq_mod
Note: See TracBrowser for help on using the repository browser.