Changeset 1322 for palm/trunk/SOURCE/timestep.f90

Timestamp:

Mar 20, 2014 4:38:49 PM (10 years ago)

Author:

raasch

Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:

• palm/trunk/SOURCE/timestep.f90 (modified) (3 diffs)

palm/trunk/SOURCE/timestep.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! REAL functions provided with KIND-attribute
 ! 
 ! Former revisions:
@@ -156 +156 @@
              v_gtrans = v_gtrans + hom(k,1,2,0)
           ENDDO
-          u_gtrans = u_gtrans / REAL( nzt - nzb )
-          v_gtrans = v_gtrans / REAL( nzt - nzb )
+          u_gtrans = u_gtrans / REAL( nzt - nzb, KIND=wp )
+          v_gtrans = v_gtrans / REAL( nzt - nzb, KIND=wp )
        ELSE
 !
@@ -173 +173 @@
           ENDDO
           !$acc end parallel
-          uv_gtrans_l(1) = u_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
-          uv_gtrans_l(2) = v_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
+          uv_gtrans_l(1) = u_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
+          uv_gtrans_l(2) = v_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
 #if defined( __parallel )
           IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
           CALL MPI_ALLREDUCE( uv_gtrans_l, uv_gtrans, 2, MPI_REAL, MPI_SUM, &
                               comm2d, ierr )
-          u_gtrans = uv_gtrans(1) / REAL( numprocs )
-          v_gtrans = uv_gtrans(2) / REAL( numprocs )
+          u_gtrans = uv_gtrans(1) / REAL( numprocs, KIND=wp )
+          v_gtrans = uv_gtrans(2) / REAL( numprocs, KIND=wp )
 #else
           u_gtrans = uv_gtrans_l(1)