Changeset 1353 for palm/trunk/SOURCE/data_output_ptseries.f90

Timestamp:

Apr 8, 2014 3:21:23 PM (10 years ago)

Author:

heinze

Message:

REAL constants provided with KIND-attribute

File:

• palm/trunk/SOURCE/data_output_ptseries.f90 (modified) (7 diffs)

palm/trunk/SOURCE/data_output_ptseries.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! REAL constants provided with KIND-attribute
 ! 
 ! Former revisions:
@@ -111 +111 @@
               pts_value_l(0:number_of_particle_groups,dopts_num) )
 
-    pts_value_l = 0.0
-    pts_value_l(:,16) = 9999999.9    ! for calculation of minimum radius
+    pts_value_l = 0.0_wp
+    pts_value_l(:,16) = 9999999.9_wp    ! for calculation of minimum radius
 
 !
@@ -135 +135 @@
           pts_value_l(0,11) = pts_value_l(0,11) + particles(n)%rvar3 ! mean sgsw
        ENDIF
-       IF ( particles(n)%speed_z > 0.0 )  THEN
-          pts_value_l(0,12) = pts_value_l(0,12) + 1.0  ! # of upward moving prts
+       IF ( particles(n)%speed_z > 0.0_wp )  THEN
+          pts_value_l(0,12) = pts_value_l(0,12) + 1.0_wp  ! # of upward moving prts
           pts_value_l(0,13) = pts_value_l(0,13) + &
                                             particles(n)%speed_z ! mean w upw.
@@ -170 +170 @@
              pts_value_l(j,11) = pts_value_l(j,11) + particles(n)%rvar3
           ENDIF
-          IF ( particles(n)%speed_z > 0.0 )  THEN
-             pts_value_l(j,12) = pts_value_l(j,12) + 1.0
+          IF ( particles(n)%speed_z > 0.0_wp )  THEN
+             pts_value_l(j,12) = pts_value_l(j,12) + 1.0_wp
              pts_value_l(j,13) = pts_value_l(j,13) + particles(n)%speed_z
           ELSE
@@ -179 +179 @@
           pts_value_l(j,16) = MIN( pts_value(j,16), particles(n)%radius )
           pts_value_l(j,17) = MAX( pts_value(j,17), particles(n)%radius )
-          pts_value_l(j,18) = pts_value_l(j,18) + 1.0
-          pts_value_l(j,19) = pts_value_l(j,19) + 1.0
+          pts_value_l(j,18) = pts_value_l(j,18) + 1.0_wp
+          pts_value_l(j,19) = pts_value_l(j,19) + 1.0_wp
 
        ENDIF
@@ -221 +221 @@
     DO  j = 0, inum
 
-       IF ( pts_value(j,1) > 0.0 )  THEN
+       IF ( pts_value(j,1) > 0.0_wp )  THEN
 
           pts_value(j,2:15) = pts_value(j,2:15) / pts_value(j,1)
-          IF ( pts_value(j,12) > 0.0  .AND.  pts_value(j,12) < 1.0 )  THEN
+          IF ( pts_value(j,12) > 0.0_wp  .AND.  pts_value(j,12) < 1.0_wp )  THEN
              pts_value(j,13) = pts_value(j,13) / pts_value(j,12)
-             pts_value(j,14) = pts_value(j,14) / ( 1.0 - pts_value(j,12) )
-          ELSEIF ( pts_value(j,12) == 0.0 )  THEN
-             pts_value(j,13) = -1.0
+             pts_value(j,14) = pts_value(j,14) / ( 1.0_wp - pts_value(j,12) )
+          ELSEIF ( pts_value(j,12) == 0.0_wp )  THEN
+             pts_value(j,13) = -1.0_wp
           ELSE
-             pts_value(j,14) = -1.0
+             pts_value(j,14) = -1.0_wp
           ENDIF
 
@@ -323 +323 @@
     DO  j = 0, inum
 
-       IF ( pts_value(j,1) > 0.0 )  THEN
+       IF ( pts_value(j,1) > 0.0_wp )  THEN
           pts_value(j,20:28) = pts_value(j,20:28) / pts_value(j,1)
        ENDIF