Changeset 1353 for palm/trunk/SOURCE/data_output_mask.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_mask.f90
r1329 r1353 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 218 218 DO k = nzb, nzt+1 219 219 psi = prt_start_index(k,j,i) 220 s_r3 = 0.0 221 s_r4 = 0.0 220 s_r3 = 0.0_wp 221 s_r4 = 0.0_wp 222 222 DO n = psi, psi+prt_count(k,j,i)-1 223 223 s_r3 = s_r3 + particles(n)%radius**3 * & … … 226 226 particles(n)%weight_factor 227 227 ENDDO 228 IF ( s_r3 /= 0.0 ) THEN228 IF ( s_r3 /= 0.0_wp ) THEN 229 229 mean_r = s_r4 / s_r3 230 230 ELSE 231 mean_r = 0.0 231 mean_r = 0.0_wp 232 232 ENDIF 233 233 tend(k,j,i) = mean_r
Note: See TracChangeset
for help on using the changeset viewer.