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

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

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

    r1321 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    102102             DO  k = nzb_s_inner(j,i)+1, nzt
    103103                tend(k,j,i) = tend(k,j,i)                                      &
    104               -0.5 * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &
    105                      - ( u(k,j,i)   - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) &
    106                      ) * ddx                                                   &
    107               -0.5 * ( ( v(k,j+1,i) - v_gtrans ) * ( sk(k,j+1,i) - sk(k,j,i) ) &
    108                      - ( v(k,j,i)   - v_gtrans ) * ( sk(k,j-1,i) - sk(k,j,i) ) &
    109                      ) * ddy                                                   &
    110               -      (   w(k,j,i)                * ( sk(k+1,j,i) - sk(k,j,i) ) &
    111                      -   w(k-1,j,i)              * ( sk(k-1,j,i) - sk(k,j,i) ) &
    112                      ) * dd2zu(k)
     104              -0.5_wp * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &
     105                        - ( u(k,j,i)   - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) &
     106                        ) * ddx                                                   &
     107              -0.5_wp * ( ( v(k,j+1,i) - v_gtrans ) * ( sk(k,j+1,i) - sk(k,j,i) ) &
     108                        - ( v(k,j,i)   - v_gtrans ) * ( sk(k,j-1,i) - sk(k,j,i) ) &
     109                        ) * ddy                                                   &
     110              -         (   w(k,j,i)                * ( sk(k+1,j,i) - sk(k,j,i) ) &
     111                        -   w(k-1,j,i)              * ( sk(k-1,j,i) - sk(k,j,i) ) &
     112                        ) * dd2zu(k)
    113113             ENDDO
    114114          ENDDO
     
    153153       DO  k = nzb_s_inner(j,i)+1, nzt
    154154          tend(k,j,i) = tend(k,j,i)                                            &
    155               -0.5 * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &
    156                      - ( u(k,j,i)   - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) &
    157                      ) * ddx                                                   &
    158               -0.5 * ( ( v(k,j+1,i) - v_gtrans ) * ( sk(k,j+1,i) - sk(k,j,i) ) &
    159                      - ( v(k,j,i)   - v_gtrans ) * ( sk(k,j-1,i) - sk(k,j,i) ) &
    160                      ) * ddy                                                   &
    161               -      (   w(k,j,i)                * ( sk(k+1,j,i) - sk(k,j,i) ) &
    162                      -   w(k-1,j,i)              * ( sk(k-1,j,i) - sk(k,j,i) ) &
    163                      ) * dd2zu(k)
     155              -0.5_wp * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &
     156                        - ( u(k,j,i)   - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) &
     157                        ) * ddx                                                   &
     158              -0.5_wp * ( ( v(k,j+1,i) - v_gtrans ) * ( sk(k,j+1,i) - sk(k,j,i) ) &
     159                        - ( v(k,j,i)   - v_gtrans ) * ( sk(k,j-1,i) - sk(k,j,i) ) &
     160                        ) * ddy                                                   &
     161              -         (   w(k,j,i)                * ( sk(k+1,j,i) - sk(k,j,i) ) &
     162                        -   w(k-1,j,i)              * ( sk(k-1,j,i) - sk(k,j,i) ) &
     163                        ) * dd2zu(k)
    164164       ENDDO
    165165
Note: See TracChangeset for help on using the changeset viewer.