Ignore:
Timestamp:
Mar 8, 2007 4:12:41 AM (15 years ago)
Author:
raasch
Message:

small bugs removed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/wall_fluxes.f90

    r54 r55  
    190190
    191191    INTEGER ::  i, j, k, kk, nzb_w, nzt_w, wall_index
    192     REAL    ::  a, b, c1, c2, h1, h2, zp
     192    REAL    ::  a, b, c1, c2, h1, h2, vel_zp, zp
    193193
    194194    REAL ::  rifs
     
    224224!
    225225!--    (5) Compute wall_flux (u'v', v'u', w'v', or w'u')
     226       vel_zp = 0.5 * (       a * ( u(k,j,i) + u(k,j,i+1) ) +  &
     227                              b * ( v(k,j,i) + v(k,j+1,i) ) +  &
     228                        (c1+c2) * ( w(k,j,i) + w(k-1,j,i) )    &
     229                      )
     230
    226231       IF ( rifs >= 0.0 )  THEN
    227232
    228233!
    229234!--       Stable stratification (and neutral)
    230           wall_flux(k) = kappa *                                          &
    231                          ( a*u(k,j,i) + b*v(k,j,i) + (c1+c2)*w(k,j,i) ) / &
    232                          (  LOG( zp / z0(j,i) ) +                         &
    233                             5.0 * rifs * ( zp - z0(j,i) ) / zp            &
    234                          )
     235          wall_flux(k) = kappa *  vel_zp / &
     236                         ( LOG( zp/z0(j,i) ) + 5.0*rifs * ( zp-z0(j,i) ) / zp )
    235237       ELSE
    236238
     
    245247!--       argument of the logarithm.
    246248          IF ( h1 == 1.0  .OR.  h2 == 1.0 )  THEN
    247              wall_flux(k) = kappa *                                          &
    248                             ( a*u(k,j,i) + b*v(k,j,i) + (c1+c2)*w(k,j,i) ) / &
     249             wall_flux(k) = kappa * vel_zp /                                 &
    249250                            ( LOG( zp / z0(j,i) ) +                          &
    250251                              5.0 * rifs * ( zp - z0(j,i) ) / zp             &
    251252                            )
    252253          ELSE
    253              wall_flux(k) = kappa *                                            &
    254                             ( a*u(k,j,i) + b*v(k,j,i) + (c1+c2)*w(k,j,i) ) /   &
     254             wall_flux(k) = kappa * vel_zp /                                   &
    255255                            ( LOG( (1.0+h2) / (1.0-h2) * (1.0-h1) / (1.0+h1) ) &
    256256                              + 2.0 * ( ATAN( h2 ) - ATAN( h1 ) )              &
Note: See TracChangeset for help on using the changeset viewer.