Changeset 1353 for palm/trunk/SOURCE/compute_vpt.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/compute_vpt.f90
r1321 r1353 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 68 68 69 69 IF ( .NOT. cloud_physics .AND. .NOT. cloud_droplets ) THEN 70 vpt = pt * ( 1.0 + 0.61* q )70 vpt = pt * ( 1.0_wp + 0.61_wp * q ) 71 71 ELSE IF (cloud_physics) THEN 72 72 DO k = nzb, nzt+1 73 73 vpt(k,:,:) = ( pt(k,:,:) + pt_d_t(k) * l_d_cp * ql(k,:,:) ) * & 74 ( 1.0 + 0.61 * q(k,:,:) - 1.61* ql(k,:,:) )74 ( 1.0_wp + 0.61_wp * q(k,:,:) - 1.61_wp * ql(k,:,:) ) 75 75 ENDDO 76 76 ELSE 77 vpt = pt * ( 1.0 + 0.61* q - ql )77 vpt = pt * ( 1.0_wp + 0.61_wp * q - ql ) 78 78 ENDIF 79 79
Note: See TracChangeset
for help on using the changeset viewer.