Ignore:
Timestamp:
Jan 13, 2016 8:19:40 AM (6 years ago)
Author:
raasch
Message:

unnecessary calculations of kmzm and kmzp removed from wall bounded parts

File:
1 edited

Legend:

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

    r1683 r1740  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! unnecessary calculations of kmzm and kmzp in wall bounded parts removed
    2222!
    2323! Former revisions:
     
    220220                kmzp = 0.25_wp *                                               &
    221221                       ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    222                 kmzm = 0.25_wp *                                               &
    223                        ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )
    224222
    225223                tend(k,j,i) = tend(k,j,i)                                      &
     
    238236!
    239237!--             Interpolate eddy diffusivities on staggered gridpoints
    240                 kmzp = 0.25_wp *                                               &
    241                        ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    242238                kmzm = 0.25_wp *                                               &
    243239                       ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )
     
    405401                kmzp = 0.25_wp *                                               &
    406402                       ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    407                 kmzm = 0.25_wp *                                               &
    408                        ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )
    409403
    410404                tend(k,j,i) = tend(k,j,i)                                      &
     
    431425!
    432426!--             Interpolate eddy diffusivities on staggered gridpoints
    433                 kmzp = 0.25_wp *                                               &
    434                        ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    435427                kmzm = 0.25_wp *                                               &
    436428                       ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )
     
    571563!--       Interpolate eddy diffusivities on staggered gridpoints
    572564          kmzp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    573           kmzm = 0.25_wp * ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )
    574565
    575566          tend(k,j,i) = tend(k,j,i)                                            &
     
    588579!
    589580!--       Interpolate eddy diffusivities on staggered gridpoints
    590           kmzp = 0.25_wp * &
    591                  ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    592581          kmzm = 0.25_wp * &
    593582                 ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )
Note: See TracChangeset for help on using the changeset viewer.