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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! REAL constants provided with KIND-attribute
 ! 
 ! Former revisions:
@@ -121 +121 @@
                    ALLOCATE( e_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                e_av = 0.0
+                e_av = 0.0_wp
 
              CASE ( 'lpt' )
@@ -127 +127 @@
                    ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                lpt_av = 0.0
+                lpt_av = 0.0_wp
 
              CASE ( 'lwp*' )
@@ -133 +133 @@
                    ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                lwp_av = 0.0
+                lwp_av = 0.0_wp
 
              CASE ( 'nr' )
@@ -139 +139 @@
                    ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                nr_av = 0.0
+                nr_av = 0.0_wp
 
              CASE ( 'p' )
@@ -145 +145 @@
                    ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                p_av = 0.0
+                p_av = 0.0_wp
 
              CASE ( 'pc' )
@@ -151 +151 @@
                    ALLOCATE( pc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                pc_av = 0.0
+                pc_av = 0.0_wp
 
              CASE ( 'pr' )
@@ -157 +157 @@
                    ALLOCATE( pr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                pr_av = 0.0
+                pr_av = 0.0_wp
 
              CASE ( 'prr' )
@@ -163 +163 @@
                    ALLOCATE( prr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                prr_av = 0.0
+                prr_av = 0.0_wp
 
              CASE ( 'prr*' )
@@ -169 +169 @@
                    ALLOCATE( precipitation_rate_av(nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                precipitation_rate_av = 0.0
+                precipitation_rate_av = 0.0_wp
 
              CASE ( 'pt' )
@@ -175 +175 @@
                    ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                pt_av = 0.0
+                pt_av = 0.0_wp
 
              CASE ( 'q' )
@@ -181 +181 @@
                    ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                q_av = 0.0
+                q_av = 0.0_wp
 
              CASE ( 'qc' )
@@ -187 +187 @@
                    ALLOCATE( qc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                qc_av = 0.0
+                qc_av = 0.0_wp
 
              CASE ( 'ql' )
@@ -193 +193 @@
                    ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                ql_av = 0.0
+                ql_av = 0.0_wp
 
              CASE ( 'ql_c' )
@@ -199 +199 @@
                    ALLOCATE( ql_c_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                ql_c_av = 0.0
+                ql_c_av = 0.0_wp
 
              CASE ( 'ql_v' )
@@ -205 +205 @@
                    ALLOCATE( ql_v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                ql_v_av = 0.0
+                ql_v_av = 0.0_wp
 
              CASE ( 'ql_vp' )
@@ -211 +211 @@
                    ALLOCATE( ql_vp_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                ql_vp_av = 0.0
+                ql_vp_av = 0.0_wp
 
              CASE ( 'qr' )
@@ -217 +217 @@
                    ALLOCATE( qr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                qr_av = 0.0
+                qr_av = 0.0_wp
 
              CASE ( 'qsws*' )
@@ -223 +223 @@
                    ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                qsws_av = 0.0
+                qsws_av = 0.0_wp
 
              CASE ( 'qv' )
@@ -229 +229 @@
                    ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                qv_av = 0.0
+                qv_av = 0.0_wp
 
              CASE ( 'rho' )
@@ -235 +235 @@
                    ALLOCATE( rho_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                rho_av = 0.0
+                rho_av = 0.0_wp
 
              CASE ( 's' )
@@ -241 +241 @@
                    ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                s_av = 0.0
+                s_av = 0.0_wp
 
              CASE ( 'sa' )
@@ -247 +247 @@
                    ALLOCATE( sa_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                sa_av = 0.0
+                sa_av = 0.0_wp
 
              CASE ( 'shf*' )
@@ -253 +253 @@
                    ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                shf_av = 0.0
+                shf_av = 0.0_wp
 
              CASE ( 't*' )
@@ -259 +259 @@
                    ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                ts_av = 0.0
+                ts_av = 0.0_wp
 
              CASE ( 'u' )
@@ -265 +265 @@
                    ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                u_av = 0.0
+                u_av = 0.0_wp
 
              CASE ( 'u*' )
@@ -271 +271 @@
                    ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                us_av = 0.0
+                us_av = 0.0_wp
 
              CASE ( 'v' )
@@ -277 +277 @@
                    ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                v_av = 0.0
+                v_av = 0.0_wp
 
              CASE ( 'vpt' )
@@ -283 +283 @@
                    ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                vpt_av = 0.0
+                vpt_av = 0.0_wp
 
              CASE ( 'w' )
@@ -289 +289 @@
                    ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                w_av = 0.0
+                w_av = 0.0_wp
 
              CASE ( 'z0*' )
@@ -295 +295 @@
                    ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                z0_av = 0.0
+                z0_av = 0.0_wp
 
              CASE ( 'z0h*' )
@@ -301 +301 @@
                    ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
                 ENDIF
-                z0h_av = 0.0
+                z0h_av = 0.0_wp
 
              CASE DEFAULT
@@ -380 +380 @@
                    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 * &
@@ -388 +388 @@
                                        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
                       pr_av(k,j,i) = pr_av(k,j,i) + mean_r