Changeset 1320 for palm/trunk/SOURCE/user_statistics.f90

Timestamp:

Mar 20, 2014 8:40:49 AM (10 years ago)

Author:

raasch

Message:

ONLY-attribute added to USE-statements,
kind-parameters added to all INTEGER and REAL declaration statements,
kinds are defined in new module kinds,
old module precision_kind is removed,
revision history before 2012 removed,
comment fields (!:) to be used for variable explanations added to all variable declaration statements

File:

• palm/trunk/SOURCE/user_statistics.f90 (modified) (5 diffs)

palm/trunk/SOURCE/user_statistics.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
+! kind-parameters added to all INTEGER and REAL declaration statements, 
+! kinds are defined in new module kinds, 
+! old module precision_kind is removed,
+! revision history before 2012 removed,
+! comment fields (!:) to be used for variable explanations added to
+! all variable declaration statements 
+!
 !
 ! Former revisions:
@@ -30 +37 @@
 ! 1036 2012-10-22 13:43:42Z raasch
 ! code put under GPL (PALM 3.9)
-!
-! 678 2011-02-02 14:31:56Z raasch
-! Bugfix: module pegrid was missing
-!
-! 622 2010-12-10 08:08:13Z raasch
-! optional barriers included in order to speed up collective operations
-!
-! 555 2010-09-07 07:32:53Z raasch
-! Bugfix: wrong dimension used for ts_value_l
-!
-! 493 2010-03-01 08:30:24Z raasch
-! Bugfix: timeseries data have to be collected by PE0
 !
 ! 211 2008-11-11 04:46:24Z raasch
@@ -57 +52 @@
 
     USE arrays_3d
+    
     USE indices
+    
+    USE kinds
+    
     USE netcdf_control
+    
     USE pegrid
+    
     USE statistics
+    
     USE user
 
     IMPLICIT NONE
 
-    CHARACTER (LEN=*) ::  mode
+    CHARACTER (LEN=*) ::  mode   !: 
 
-    INTEGER ::  i, j, k, sr, tn
+    INTEGER(iwp) ::  i    !: 
+    INTEGER(iwp) ::  j    !: 
+    INTEGER(iwp) ::  k    !: 
+    INTEGER(iwp) ::  sr   !: 
+    INTEGER(iwp) ::  tn   !: 
 
-    REAL, DIMENSION(dots_num_palm+1:dots_max) ::  ts_value_l
+    REAL(wp),                                                                  &
+       DIMENSION(dots_num_palm+1:dots_max) ::                                  &
+          ts_value_l   !: 
 
 
@@ -87 +95 @@
 !!--             Sample on how to calculate the profile of the resolved-scale
 !!--             horizontal momentum flux u*v*
-!                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +           &
-!                      ( 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) * &
-!                      ( 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) )   &
+!                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +               &
+!                      ( 0.5 * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) *     &
+!                      ( 0.5 * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) )       &
 !                                                 * rmask(j,i,sr)
 !!
 !!--             Further profiles can be defined and calculated by increasing
 !!--             the second index of array sums_l (replace ... appropriately)
-!                sums_l(k,pr_palm+2,tn) = sums_l(k,pr_palm+2,tn) + ... &
+!                sums_l(k,pr_palm+2,tn) = sums_l(k,pr_palm+2,tn) + ...           &
 !                                         * rmask(j,i,sr)
 !             ENDDO
@@ -121 +129 @@
 !#if defined( __parallel ) && ! defined ( __check )
 !       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
-!       CALL MPI_ALLREDUCE( ts_value_l(dots_num_palm+1),                       &
-!                           ts_value(dots_num_palm+1,sr),                      &
-!                           dots_max-dots_num_palm, MPI_REAL, MPI_SUM, comm2d, &
+!       CALL MPI_ALLREDUCE( ts_value_l(dots_num_palm+1),                         &
+!                           ts_value(dots_num_palm+1,sr),                        &
+!                           dots_max-dots_num_palm, MPI_REAL, MPI_SUM, comm2d,   &
 !                           ierr )
 !#else