Changeset 1740 for palm/trunk/SOURCE/diffusion_u.f90

Timestamp:

Jan 13, 2016 8:19:40 AM (9 years ago)

Author:

raasch

Message:

unnecessary calculations of kmzm and kmzp removed from wall bounded parts

File:

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

palm/trunk/SOURCE/diffusion_u.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! unnecessary calculations of kmzm and kmzp in wall bounded parts removed
 ! 
 ! Former revisions:
@@ -225 +225 @@
                 kmzp = 0.25_wp *                                               &
                       ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
-                kmzm = 0.25_wp *                                               &
-                      ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
 
                 tend(k,j,i) = tend(k,j,i)                                      &
@@ -243 +241 @@
 !
 !--             Interpolate eddy diffusivities on staggered gridpoints
-                kmzp = 0.25_wp *                                               &
-                       ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
                 kmzm = 0.25_wp *                                               &
                        ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
@@ -410 +406 @@
                 kmzp = 0.25_wp *                                               &
                       ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
-                kmzm = 0.25_wp *                                               &
-                      ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
 
                 tend(k,j,i) = tend(k,j,i)                                      &
@@ -436 +430 @@
 !
 !--             Interpolate eddy diffusivities on staggered gridpoints
-                kmzp = 0.25_wp *                                               &
-                       ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
                 kmzm = 0.25_wp *                                               &
                        ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
@@ -574 +566 @@
 !--       Interpolate eddy diffusivities on staggered gridpoints
           kmzp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
-          kmzm = 0.25_wp * ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
 
           tend(k,j,i) = tend(k,j,i)                                            &
@@ -591 +582 @@
 !
 !--       Interpolate eddy diffusivities on staggered gridpoints
-          kmzp = 0.25_wp *                                                     &
-                 ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
           kmzm = 0.25_wp *                                                     &
                  ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )