Changeset 1353 for palm/trunk/SOURCE/data_output_3d.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/data_output_3d.f90
r1329 r1353 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 257 257 DO k = nzb, nz_do3d 258 258 psi = prt_start_index(k,j,i) 259 s_r3 = 0.0 260 s_r4 = 0.0 259 s_r3 = 0.0_wp 260 s_r4 = 0.0_wp 261 261 DO n = psi, psi+prt_count(k,j,i)-1 262 262 s_r3 = s_r3 + particles(n)%radius**3 * & … … 265 265 particles(n)%weight_factor 266 266 ENDDO 267 IF ( s_r3 /= 0.0 ) THEN267 IF ( s_r3 /= 0.0_wp ) THEN 268 268 mean_r = s_r4 / s_r3 269 269 ELSE 270 mean_r = 0.0 270 mean_r = 0.0_wp 271 271 ENDIF 272 272 tend(k,j,i) = mean_r
Note: See TracChangeset
for help on using the changeset viewer.