Ignore:
Timestamp:
May 30, 2017 5:47:52 PM (4 years ago)
Author:
suehring
Message:

Adjustments according new topography and surface-modelling concept implemented

File:
1 edited

Legend:

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

    r2101 r2232  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! topography representation via flags
    2323!
    2424! Former revisions:
     
    9696
    9797       USE indices,                                                            &
    98            ONLY:  nxlu, nxr, nyn, nys, nzb_u_inner, nzt
     98           ONLY:  nxlu, nxr, nyn, nys, nzb, nzt, wall_flags_0
    9999
    100100       USE kinds
     
    114114       DO  i = nxlu, nxr
    115115          DO  j = nys, nyn
    116              DO  k = nzb_u_inner(j,i)+1, nzt
     116             DO  k = nzb+1, nzt
    117117                tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                        &
    118118                         ( u(k,j,i+1) * ( u(k,j,i+1) + u(k,j,i) - gu )         &
     
    123123                         - u(k-1,j,i) * ( w(k-1,j,i) + w(k-1,j,i-1) ) )        &
    124124                                                                  * ddzw(k)    &
    125                                                       )
     125                                                      )                        &
     126                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     127                                               BTEST( wall_flags_0(k,j,i), 1 ) )
    126128             ENDDO
    127129          ENDDO
     
    148150
    149151       USE indices,                                                            &
    150            ONLY:  nzb_u_inner, nzt
     152           ONLY:  nzb, nzt, wall_flags_0
    151153
    152154       USE kinds
     
    164166       gu = 2.0_wp * u_gtrans
    165167       gv = 2.0_wp * v_gtrans
    166        DO  k = nzb_u_inner(j,i)+1, nzt
     168       DO  k = nzb+1, nzt
    167169          tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                              &
    168170                         ( u(k,j,i+1) * ( u(k,j,i+1) + u(k,j,i) - gu )         &
     
    173175                         - u(k-1,j,i) * ( w(k-1,j,i) + w(k-1,j,i-1) ) )        &
    174176                                                                  * ddzw(k)    &
    175                                                 )
     177                                                )                              &
     178                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     179                                               BTEST( wall_flags_0(k,j,i), 1 ) )
    176180       ENDDO
    177181
Note: See TracChangeset for help on using the changeset viewer.