Ignore:
Timestamp:
May 30, 2017 5:47:52 PM (7 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_v_pw.f90

    r2101 r2232  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! topography representation via flags
    2323!
    2424! Former revisions:
     
    9696
    9797       USE indices,                                                            &
    98            ONLY:  nxl, nxr, nyn, nysv, nzb_v_inner, nzt
     98           ONLY:  nxl, nxr, nyn, nysv, nzb, nzt, wall_flags_0
    9999
    100100       USE kinds
     
    115115       DO  i = nxl, nxr
    116116          DO  j = nysv, nyn
    117              DO  k = nzb_v_inner(j,i)+1, nzt
     117             DO  k = nzb+1, nzt
    118118                tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                        &
    119119                         ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu )     &
     
    124124                         - v(k-1,j,i) * ( w(k-1,j-1,i) + w(k-1,j,i) ) )        &
    125125                                                                  * ddzw(k)    &
    126                                                       )
     126                                                      )                        &
     127                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     128                                               BTEST( wall_flags_0(k,j,i), 2 ) )
    127129             ENDDO
    128130          ENDDO
     
    149151
    150152       USE indices,                                                            &
    151            ONLY:  nzb_v_inner, nzt
     153           ONLY:  nzb, nzt, wall_flags_0
    152154
    153155       USE kinds
     
    166168       gu = 2.0_wp * u_gtrans
    167169       gv = 2.0_wp * v_gtrans
    168        DO  k = nzb_v_inner(j,i)+1, nzt
     170       DO  k = nzb+1, nzt
    169171          tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                              &
    170172                         ( v(k,j,i+1) * ( u(k,j-1,i+1) + u(k,j,i+1) - gu )     &
     
    175177                         - v(k-1,j,i) * ( w(k-1,j-1,i) + w(k-1,j,i) ) )        &
    176178                                                                  * ddzw(k)    &
    177                                                 )
     179                                                )                              &
     180                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     181                                               BTEST( wall_flags_0(k,j,i), 2 ) )
    178182       ENDDO
    179183 
Note: See TracChangeset for help on using the changeset viewer.