Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

Ignore:
Timestamp:
Mar 20, 2014 4:38:49 PM (10 years ago)
Author:
raasch
Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:
1 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/SOURCE/diffusion_w.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
    … …  
    132132!
    133133!--             Interpolate eddy diffusivities on staggered gridpoints
    134                 kmxp = 0.25 * &
     134                kmxp = 0.25_wp *                                              &
    135135                       ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    136                 kmxm = 0.25 * &
     136                kmxm = 0.25_wp *                                              &
    137137                       ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    138                 kmyp = 0.25 * &
     138                kmyp = 0.25_wp *                                              &
    139139                       ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    140                 kmym = 0.25 * &
     140                kmym = 0.25_wp *                                              &
    141141                       ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    142142
    … …  
    152152                      &   - kmym * ( v(k+1,j,i) - v(k,j,i)   ) * ddzu(k+1)     &
    153153                      &   ) * ddy                                              &
    154                       & + 2.0 * (                                              &
     154                      & + 2.0_wp * (                                           &
    155155                      &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1)  &
    156156                      & - km(k,j,i)   * ( w(k,j,i)   - w(k-1,j,i) ) * ddzw(k)  &
    … …  
    165165!
    166166!--                Interpolate eddy diffusivities on staggered gridpoints
    167                    kmxp = 0.25 * &
     167                   kmxp = 0.25_wp *                                            &
    168168                          ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    169                    kmxm = 0.25 * &
     169                   kmxm = 0.25_wp *                                            &
    170170                          ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    171                    kmyp = 0.25 * &
     171                   kmyp = 0.25_wp *                                            &
    172172                          ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    173                    kmym = 0.25 * &
     173                   kmym = 0.25_wp *                                            &
    174174                          ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    175175
    … …  
    262262!
    263263!--                Interpolate eddy diffusivities on staggered gridpoints
    264                    kmxp = 0.25 * &
     264                   kmxp = 0.25_wp *                                            &
    265265                          ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    266                    kmxm = 0.25 * &
     266                   kmxm = 0.25_wp *                                            &
    267267                          ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    268                    kmyp = 0.25 * &
     268                   kmyp = 0.25_wp *                                            &
    269269                          ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    270                    kmym = 0.25 * &
     270                   kmym = 0.25_wp *                                            &
    271271                          ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    272272
    … …  
    282282                         &   - kmym * ( v(k+1,j,i) - v(k,j,i)   ) * ddzu(k+1)    &
    283283                         &   ) * ddy                                             &
    284                          & + 2.0 * (                                             &
     284                         & + 2.0_wp * (                                          &
    285285                         &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
    286286                         & - km(k,j,i)   * ( w(k,j,i)   - w(k-1,j,i) ) * ddzw(k) &
    … …  
    297297!
    298298!--                Interpolate eddy diffusivities on staggered gridpoints
    299                    kmxp = 0.25 * &
     299                   kmxp = 0.25_wp *                                            &
    300300                          ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    301                    kmxm = 0.25 * &
     301                   kmxm = 0.25_wp *                                            &
    302302                          ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    303                    kmyp = 0.25 * &
     303                   kmyp = 0.25_wp *                                            &
    304304                          ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    305                    kmym = 0.25 * &
     305                   kmym = 0.25_wp *                                            &
    306306                          ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    307307
    … …  
    378378!
    379379!--       Interpolate eddy diffusivities on staggered gridpoints
    380           kmxp = 0.25 * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    381           kmxm = 0.25 * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    382           kmyp = 0.25 * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    383           kmym = 0.25 * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
     380          kmxp = 0.25_wp * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
     381          kmxm = 0.25_wp * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
     382          kmyp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
     383          kmym = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    384384
    385385          tend(k,j,i) = tend(k,j,i)                                            &
    … …  
    394394                      &   - kmym * ( v(k+1,j,i) - v(k,j,i)   ) * ddzu(k+1)     &
    395395                      &   ) * ddy                                              &
    396                       & + 2.0 * (                                              &
     396                      & + 2.0_wp * (                                           &
    397397                      &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1)  &
    398398                      & - km(k,j,i)   * ( w(k,j,i)   - w(k-1,j,i) ) * ddzw(k)  &
    … …  
    423423!
    424424!--          Interpolate eddy diffusivities on staggered gridpoints
    425              kmxp = 0.25 * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    426              kmxm = 0.25 * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    427              kmyp = 0.25 * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    428              kmym = 0.25 * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
     425             kmxp = 0.25_wp * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
     426             kmxm = 0.25_wp * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
     427             kmyp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
     428             kmym = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    429429
    430430             tend(k,j,i) = tend(k,j,i)                                         &
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |