Changeset 1740


Ignore:
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

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r1692 r1740  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! unnecessary calculations of kmzm and kmzp in wall bounded parts removed
    2222!
    2323! Former revisions:
     
    225225                kmzp = 0.25_wp *                                               &
    226226                      ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
    227                 kmzm = 0.25_wp *                                               &
    228                       ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
    229227
    230228                tend(k,j,i) = tend(k,j,i)                                      &
     
    243241!
    244242!--             Interpolate eddy diffusivities on staggered gridpoints
    245                 kmzp = 0.25_wp *                                               &
    246                        ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
    247243                kmzm = 0.25_wp *                                               &
    248244                       ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
     
    410406                kmzp = 0.25_wp *                                               &
    411407                      ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
    412                 kmzm = 0.25_wp *                                               &
    413                       ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
    414408
    415409                tend(k,j,i) = tend(k,j,i)                                      &
     
    436430!
    437431!--             Interpolate eddy diffusivities on staggered gridpoints
    438                 kmzp = 0.25_wp *                                               &
    439                        ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
    440432                kmzm = 0.25_wp *                                               &
    441433                       ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
     
    574566!--       Interpolate eddy diffusivities on staggered gridpoints
    575567          kmzp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
    576           kmzm = 0.25_wp * ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
    577568
    578569          tend(k,j,i) = tend(k,j,i)                                            &
     
    591582!
    592583!--       Interpolate eddy diffusivities on staggered gridpoints
    593           kmzp = 0.25_wp *                                                     &
    594                  ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )
    595584          kmzm = 0.25_wp *                                                     &
    596585                 ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )
  • 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.