Changeset 1353 for palm/trunk/SOURCE/data_output_2d.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_2d.f90 (modified) (18 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! REAL constants provided with KIND-attribute
 ! 
 ! Former revisions:
@@ -235 +235 @@
              ns = ns - 1
              ALLOCATE( local_2d_sections(nxlg:nxrg,nysg:nyng,1:ns) )
-             local_2d_sections = 0.0
+             local_2d_sections = 0.0_wp
           ENDIF
 
@@ -266 +266 @@
              ALLOCATE( local_2d_sections(nxlg:nxrg,1:ns,nzb:nzt+1) )
              ALLOCATE( local_2d_sections_l(nxlg:nxrg,1:ns,nzb:nzt+1) )
-             local_2d_sections = 0.0; local_2d_sections_l = 0.0
+             local_2d_sections = 0.0_wp; local_2d_sections_l = 0.0_wp
           ENDIF
 
@@ -297 +297 @@
              ALLOCATE( local_2d_sections(1:ns,nysg:nyng,nzb:nzt+1) )
              ALLOCATE( local_2d_sections_l(1:ns,nysg:nyng,nzb:nzt+1) )
-             local_2d_sections = 0.0; local_2d_sections_l = 0.0
+             local_2d_sections = 0.0_wp; local_2d_sections_l = 0.0_wp
           ENDIF
 
@@ -400 +400 @@
                       CALL exchange_horiz( tend, nbgp )
                    ELSE
-                      tend = 0.0
+                      tend = 0.0_wp
                    ENDIF
                    DO  i = nxlg, nxrg
@@ -422 +422 @@
                             DO  k = nzb, nzt+1
                                psi = prt_start_index(k,j,i)
-                               s_r3 = 0.0
-                               s_r4 = 0.0
+                               s_r3 = 0.0_wp
+                               s_r4 = 0.0_wp
                                DO  n = psi, psi+prt_count(k,j,i)-1
                                   s_r3 = s_r3 + particles(n)%radius**3 *       &
@@ -430 +430 @@
                                                 particles(n)%weight_factor
                                ENDDO
-                               IF ( s_r3 /= 0.0 )  THEN
+                               IF ( s_r3 /= 0.0_wp )  THEN
                                   mean_r = s_r4 / s_r3
                                ELSE
-                                  mean_r = 0.0
+                                  mean_r = 0.0_wp
                                ENDIF
                                tend(k,j,i) = mean_r
@@ -441 +441 @@
                       CALL exchange_horiz( tend, nbgp )
                    ELSE
-                      tend = 0.0
+                      tend = 0.0_wp
                    END IF
                    DO  i = nxlg, nxrg
@@ -463 +463 @@
                    ENDDO
                 ENDDO
-                precipitation_amount = 0.0   ! reset for next integ. interval
+                precipitation_amount = 0.0_wp   ! reset for next integ. interval
                 resorted = .TRUE.
                 two_d = .TRUE.
@@ -608 +608 @@
                       CALL exchange_horiz( tend, nbgp )
                    ELSE
-                      tend = 0.0
+                      tend = 0.0_wp
                    END IF
                    DO  i = nxlg, nxrg
@@ -733 +733 @@
 !--             at the bottom boundary by the real surface values.
                 IF ( do2d(av,if) == 'u_xz'  .OR.  do2d(av,if) == 'u_yz' )  THEN
-                   IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0
+                   IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0_wp
                 ENDIF
 
@@ -765 +765 @@
 !--             at the bottom boundary by the real surface values.
                 IF ( do2d(av,if) == 'v_xz'  .OR.  do2d(av,if) == 'v_yz' )  THEN
-                   IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0
+                   IF ( ibc_uv_b == 0 )  local_pf(:,:,nzb) = 0.0_wp
                 ENDIF
 
@@ -898 +898 @@
                    IF ( section(is,s) == -1  .AND.  .NOT. two_d )  THEN
 
-                      local_2d = 0.0
+                      local_2d = 0.0_wp
 !
 !--                   Carry out the averaging (all data are on the PE)
@@ -909 +909 @@
                       ENDDO
 
-                      local_2d = local_2d / ( nzt -nzb + 2.0)
+                      local_2d = local_2d / ( nzt -nzb + 2.0_wp)
 
                    ELSE
@@ -1089 +1089 @@
 
                       ALLOCATE( local_2d_l(nxlg:nxrg,nzb:nzt+1) )
-                      local_2d_l = 0.0
+                      local_2d_l = 0.0_wp
                       ngp = ( nxrg-nxlg+1 ) * ( nzt-nzb+2 )
 !
@@ -1111 +1111 @@
                       local_2d = local_2d_l
 #endif
-                      local_2d = local_2d / ( ny + 1.0 )
+                      local_2d = local_2d / ( ny + 1.0_wp )
 
                       DEALLOCATE( local_2d_l )
@@ -1308 +1308 @@
 
                       ALLOCATE( local_2d_l(nysg:nyng,nzb:nzt+1) )
-                      local_2d_l = 0.0
+                      local_2d_l = 0.0_wp
                       ngp = ( nyng-nysg+1 ) * ( nzt-nzb+2 )
 !
@@ -1330 +1330 @@
                       local_2d = local_2d_l
 #endif
-                      local_2d = local_2d / ( nx + 1.0 )
+                      local_2d = local_2d / ( nx + 1.0_wp )
 
                       DEALLOCATE( local_2d_l )