Home


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

Ignore:
Timestamp:
Apr 8, 2014 3:21:23 PM (10 years ago)
Author:
heinze
Message:

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/SOURCE/coriolis.f90

    r1321 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
    … …  
    106106                DO  j = nys, nyn
    107107                   DO  k = nzb_u_inner(j,i)+1, nzt
    108                       tend(k,j,i) = tend(k,j,i) + f  *    ( 0.25 *             &
     108                      tend(k,j,i) = tend(k,j,i) + f  *    ( 0.25_wp *          &
    109109                                   ( v(k,j,i-1) + v(k,j,i) + v(k,j+1,i-1) +    &
    110110                                     v(k,j+1,i) ) - vg(k) )                    &
    111                                              - fs *    ( 0.25 *                &
     111                                                - fs *    ( 0.25_wp *          &
    112112                                   ( w(k-1,j,i-1) + w(k-1,j,i) + w(k,j,i-1) +  &
    113113                                     w(k,j,i)   ) &
    … …  
    123123                DO  j = nysv, nyn
    124124                   DO  k = nzb_v_inner(j,i)+1, nzt
    125                       tend(k,j,i) = tend(k,j,i) - f *     ( 0.25 *             &
     125                      tend(k,j,i) = tend(k,j,i) - f *     ( 0.25_wp *          &
    126126                                   ( u(k,j-1,i) + u(k,j,i) + u(k,j-1,i+1) +    &
    127127                                     u(k,j,i+1) ) - ug(k) )
    … …  
    136136                DO  j = nys, nyn
    137137                   DO  k = nzb_w_inner(j,i)+1, nzt
    138                       tend(k,j,i) = tend(k,j,i) + fs * 0.25 *                  &
     138                      tend(k,j,i) = tend(k,j,i) + fs * 0.25_wp *               &
    139139                                   ( u(k,j,i) + u(k+1,j,i) + u(k,j,i+1) +      &
    140140                                     u(k+1,j,i+1) )
    … …  
    190190                   DO  k = 1, nzt
    191191                      IF  ( k > nzb_u_inner(j,i) )  THEN
    192                          tend(k,j,i) = tend(k,j,i) + f  *    ( 0.25 *          &
     192                         tend(k,j,i) = tend(k,j,i) + f  *    ( 0.25_wp *       &
    193193                                      ( v(k,j,i-1) + v(k,j,i) + v(k,j+1,i-1) + &
    194194                                        v(k,j+1,i) ) - vg(k) )                 &
    195                                                 - fs *    ( 0.25 *             &
     195                                                   - fs *    ( 0.25_wp *       &
    196196                                      ( w(k-1,j,i-1) + w(k-1,j,i) + w(k,j,i-1) &
    197197                                        + w(k,j,i)   )                         &
    … …  
    211211                   DO  k = 1, nzt
    212212                      IF  ( k > nzb_v_inner(j,i) )  THEN
    213                          tend(k,j,i) = tend(k,j,i) - f *     ( 0.25 *          &
     213                         tend(k,j,i) = tend(k,j,i) - f *     ( 0.25_wp *       &
    214214                                      ( u(k,j-1,i) + u(k,j,i) + u(k,j-1,i+1) + &
    215215                                        u(k,j,i+1) ) - ug(k) )
    … …  
    228228                   DO  k = 1, nzt
    229229                      IF  ( k > nzb_w_inner(j,i) )  THEN
    230                          tend(k,j,i) = tend(k,j,i) + fs * 0.25 *               &
     230                         tend(k,j,i) = tend(k,j,i) + fs * 0.25_wp *            &
    231231                                      ( u(k,j,i) + u(k+1,j,i) + u(k,j,i+1) +   &
    232232                                        u(k+1,j,i+1) )
    … …  
    278278          CASE ( 1 )
    279279             DO  k = nzb_u_inner(j,i)+1, nzt
    280                 tend(k,j,i) = tend(k,j,i) + f  *    ( 0.25 *                   &
     280                tend(k,j,i) = tend(k,j,i) + f  *    ( 0.25_wp *                &
    281281                                ( v(k,j,i-1) + v(k,j,i) + v(k,j+1,i-1) +       &
    282282                                  v(k,j+1,i) ) - vg(k) )                       &
    283                                           - fs *    ( 0.25 *                   &
     283                                          - fs *    ( 0.25_wp *                &
    284284                                ( w(k-1,j,i-1) + w(k-1,j,i) + w(k,j,i-1) +     &
    285285                                  w(k,j,i)   ) )
    … …  
    290290          CASE ( 2 )
    291291             DO  k = nzb_v_inner(j,i)+1, nzt
    292                 tend(k,j,i) = tend(k,j,i) - f *     ( 0.25 *                   &
     292                tend(k,j,i) = tend(k,j,i) - f *     ( 0.25_wp *                &
    293293                                ( u(k,j-1,i) + u(k,j,i) + u(k,j-1,i+1) +       &
    294294                                  u(k,j,i+1) ) - ug(k) )
    … …  
    299299          CASE ( 3 )
    300300             DO  k = nzb_w_inner(j,i)+1, nzt
    301                 tend(k,j,i) = tend(k,j,i) + fs * 0.25 *                        &
     301                tend(k,j,i) = tend(k,j,i) + fs * 0.25_wp *                     &
    302302                                ( u(k,j,i) + u(k+1,j,i) + u(k,j,i+1) +         &
    303303                                  u(k+1,j,i+1) )
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |