Changeset 1320 for palm/trunk/SOURCE/compute_vpt.f90
- Timestamp:
- Mar 20, 2014 8:40:49 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/compute_vpt.f90
r1310 r1320 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! ONLY-attribute added to USE-statements, 23 ! kind-parameters added to all INTEGER and REAL declaration statements, 24 ! kinds are defined in new module kinds, 25 ! old module precision_kind is removed, 26 ! revision history before 2012 removed, 27 ! comment fields (!:) to be used for variable explanations added to 28 ! all variable declaration statements 23 29 ! 24 30 ! Former revisions: … … 32 38 ! Bugfix: wrong factor in calculation of vpt in case of cloud droplets 33 39 ! 34 ! 799 2011-12-21 17:48:03Z franke35 ! Bugfix: ql is now included in calculation of vpt in case of36 ! cloud droplets37 !38 ! RCS Log replace by Id keyword, revision history cleaned up39 !40 ! Revision 1.5 2001/03/30 06:58:52 raasch41 ! Translation of remaining German identifiers (variables, subroutines, etc.)42 !43 40 ! Revision 1.1 2000/04/13 14:40:53 schroeter 44 41 ! Initial revision … … 48 45 ! ------------- 49 46 ! Computation of the virtual potential temperature 50 !------------------------------------------------------------------------------ -!47 !------------------------------------------------------------------------------! 51 48 52 USE arrays_3d 53 USE indices 54 USE cloud_parameters 55 USE control_parameters 49 USE arrays_3d, & 50 ONLY: pt, q, ql, vpt 51 52 USE indices, & 53 ONLY: nzb, nzt 54 55 USE cloud_parameters, & 56 ONLY: l_d_cp, pt_d_t 57 58 USE control_parameters, & 59 ONLY: cloud_droplets, cloud_physics 60 61 USE kinds 56 62 57 63 IMPLICIT NONE 58 64 59 INTEGER :: k65 INTEGER(iwp) :: k !: 60 66 61 IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets )THEN67 IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN 62 68 vpt = pt * ( 1.0 + 0.61 * q ) 63 ELSE IF (cloud_physics) THEN69 ELSE IF (cloud_physics) THEN 64 70 DO k = nzb, nzt+1 65 vpt(k,:,:) = ( pt(k,:,:) + pt_d_t(k) * l_d_cp * ql(k,:,:) ) * &71 vpt(k,:,:) = ( pt(k,:,:) + pt_d_t(k) * l_d_cp * ql(k,:,:) ) * & 66 72 ( 1.0 + 0.61 * q(k,:,:) - 1.61 * ql(k,:,:) ) 67 73 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.