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

Last change on this file since 1 was 1, checked in by raasch, 17 years ago

Initial repository layout and content

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