Changeset 55 for palm/trunk/SOURCE/wall_fluxes.f90

Timestamp:

Mar 8, 2007 4:12:41 AM (17 years ago)

Author:

raasch

Message:

small bugs removed

File:

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

palm/trunk/SOURCE/wall_fluxes.f90

@@ -190 +190 @@
 
     INTEGER ::  i, j, k, kk, nzb_w, nzt_w, wall_index
-    REAL    ::  a, b, c1, c2, h1, h2, zp
+    REAL    ::  a, b, c1, c2, h1, h2, vel_zp, zp
 
     REAL ::  rifs
@@ -224 +224 @@
 !
 !--    (5) Compute wall_flux (u'v', v'u', w'v', or w'u')
+       vel_zp = 0.5 * (       a * ( u(k,j,i) + u(k,j,i+1) ) +  &
+                              b * ( v(k,j,i) + v(k,j+1,i) ) +  &
+                        (c1+c2) * ( w(k,j,i) + w(k-1,j,i) )    &
+                      )
+
        IF ( rifs >= 0.0 )  THEN
 
 !
 !--       Stable stratification (and neutral)
-          wall_flux(k) = kappa *                                          &
-                         ( a*u(k,j,i) + b*v(k,j,i) + (c1+c2)*w(k,j,i) ) / &
-                         (  LOG( zp / z0(j,i) ) +                         &
-                            5.0 * rifs * ( zp - z0(j,i) ) / zp            &
-                         )
+          wall_flux(k) = kappa *  vel_zp / &
+                         ( LOG( zp/z0(j,i) ) + 5.0*rifs * ( zp-z0(j,i) ) / zp )
        ELSE
 
@@ -245 +247 @@
 !--       argument of the logarithm. 
           IF ( h1 == 1.0  .OR.  h2 == 1.0 )  THEN
-             wall_flux(k) = kappa *                                          &
-                            ( a*u(k,j,i) + b*v(k,j,i) + (c1+c2)*w(k,j,i) ) / &
+             wall_flux(k) = kappa * vel_zp /                                 &
                             ( LOG( zp / z0(j,i) ) +                          &
                               5.0 * rifs * ( zp - z0(j,i) ) / zp             &
                             )
           ELSE
-             wall_flux(k) = kappa *                                            &
-                            ( a*u(k,j,i) + b*v(k,j,i) + (c1+c2)*w(k,j,i) ) /   &
+             wall_flux(k) = kappa * vel_zp /                                   &
                             ( LOG( (1.0+h2) / (1.0-h2) * (1.0-h1) / (1.0+h1) ) &
                               + 2.0 * ( ATAN( h2 ) - ATAN( h1 ) )              &