source: palm/tags/release-3.8/SOURCE/compute_vpt.f90 @ 4343

Last change on this file since 4343 was 484, checked in by raasch, 11 years ago

typo in file headers removed

  • Property svn:keywords set to Id
File size: 1.0 KB
Line 
1 SUBROUTINE compute_vpt
2
3!-------------------------------------------------------------------------------!
4! Current revisions:
5! -----------------
6!
7!
8! Former revisions:
9! -----------------
10! $Id: compute_vpt.f90 484 2010-02-05 07:36:54Z oliver.maas $
11! RCS Log replace by Id keyword, revision history cleaned up
12!
13! Revision 1.5  2001/03/30 06:58:52  raasch
14! Translation of remaining German identifiers (variables, subroutines, etc.)
15!
16! Revision 1.1  2000/04/13 14:40:53  schroeter
17! Initial revision
18!
19!
20! Description:
21! -------------
22! Computation of the virtual potential temperature
23!-------------------------------------------------------------------------------!
24
25    USE arrays_3d
26    USE indices
27    USE cloud_parameters
28    USE control_parameters
29
30    IMPLICIT NONE
31
32    INTEGER :: k
33
34    IF ( .NOT. cloud_physics ) THEN
35       vpt = pt * ( 1.0 + 0.61 * q )
36    ELSE
37       DO  k = nzb, nzt+1
38          vpt(k,:,:) = ( pt(k,:,:) + pt_d_t(k) * l_d_cp * ql(k,:,:) ) * &
39                       ( 1.0 + 0.61 * q(k,:,:) - 1.61 * ql(k,:,:) ) 
40       ENDDO
41    ENDIF
42
43 END SUBROUTINE compute_vpt
Note: See TracBrowser for help on using the repository browser.