Ignore:
Timestamp:
Aug 4, 2008 3:39:12 PM (16 years ago)
Author:
letzel
Message:
  • mrun: (optionally) link local input
  • diffusion_s: bugfix for calculation of horizontal fluxes at vertical walls
  • NCL scripts in trunk/SCRIPTS/NCL updated
File:
1 edited

Legend:

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

    r139 r183  
    44! Actual revisions:
    55! -----------------
    6 !
     6! bugfix: calculation of fluxes at vertical surfaces
    77!
    88! Former revisions:
     
    9090
    9191                   tend(k,j,i) = tend(k,j,i)                                  &
    92                                           + 0.5 * ( fwxp(j,i) *               &
     92                                                + ( fwxp(j,i) * 0.5 *         &
    9393                        ( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
    9494                        + ( 1.0 - fwxp(j,i) ) * wall_s_flux(1)                &
    95                                                    -fwxm(j,i) *               &
     95                                                   -fwxm(j,i) * 0.5 *         &
    9696                        ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
    9797                        + ( 1.0 - fwxm(j,i) ) * wall_s_flux(2)                &
    9898                                                  ) * ddx2                    &
    99                                           + 0.5 * ( fwyp(j,i) *               &
     99                                                + ( fwyp(j,i) * 0.5 *         &
    100100                        ( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
    101101                        + ( 1.0 - fwyp(j,i) ) * wall_s_flux(3)                &
    102                                                    -fwym(j,i) *               &
     102                                                   -fwym(j,i) * 0.5 *         &
    103103                        ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
    104104                        + ( 1.0 - fwym(j,i) ) * wall_s_flux(4)                &
     
    200200
    201201             tend(k,j,i) = tend(k,j,i)                                        &
    202                                           + 0.5 * ( fwxp(j,i) *               &
     202                                                + ( fwxp(j,i) * 0.5 *         &
    203203                        ( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
    204204                        + ( 1.0 - fwxp(j,i) ) * wall_s_flux(1)                &
    205                                                    -fwxm(j,i) *               &
     205                                                   -fwxm(j,i) * 0.5 *         &
    206206                        ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
    207207                        + ( 1.0 - fwxm(j,i) ) * wall_s_flux(2)                &
    208208                                                  ) * ddx2                    &
    209                                           + 0.5 * ( fwyp(j,i) *               &
     209                                                + ( fwyp(j,i) * 0.5 *         &
    210210                        ( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
    211211                        + ( 1.0 - fwyp(j,i) ) * wall_s_flux(3)                &
    212                                                    -fwym(j,i) *               &
     212                                                   -fwym(j,i) * 0.5 *         &
    213213                        ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
    214214                        + ( 1.0 - fwym(j,i) ) * wall_s_flux(4)                &
Note: See TracChangeset for help on using the changeset viewer.