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_w_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, nys, nzb_w_inner, nzt
     98           ONLY:  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
    9999
    100100       USE kinds
     
    115115       DO  i = nxl, nxr
    116116          DO  j = nys, nyn
    117              DO  k = nzb_w_inner(j,i)+1, nzt
     117             DO  k = nzb+1, nzt
    118118                tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                        &
    119119                         ( w(k,j,i+1) * ( u(k+1,j,i+1) + u(k,j,i+1) - gu )     &
     
    124124                         - w(k-1,j,i) * ( w(k,j,i) + w(k-1,j,i) ) )            &
    125125                                                                  * ddzu(k+1)  &
    126                                                       )
     126                                                      )                        &
     127                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     128                                               BTEST( wall_flags_0(k,j,i), 3 ) )
    127129             ENDDO
    128130          ENDDO
     
    149151
    150152       USE indices,                                                            &
    151            ONLY:  nzb_w_inner, nzt
     153           ONLY:  nzb, nzt, wall_flags_0
    152154
    153155       USE kinds
     
    165167       gu = 2.0_wp * u_gtrans
    166168       gv = 2.0_wp * v_gtrans
    167        DO  k = nzb_w_inner(j,i)+1, nzt
     169       DO  k = nzb+1, nzt
    168170          tend(k,j,i) = tend(k,j,i) - 0.25_wp * (                              &
    169171                         ( w(k,j,i+1) * ( u(k+1,j,i+1) + u(k,j,i+1) - gu )     &
     
    174176                         - w(k-1,j,i) * ( w(k,j,i) + w(k-1,j,i) ) )            &
    175177                                                                  * ddzu(k+1)  &
    176                                                 )
     178                                                )                              &
     179                                      * MERGE( 1.0_wp, 0.0_wp,                 &
     180                                               BTEST( wall_flags_0(k,j,i), 3 ) )
    177181       ENDDO
    178182    END SUBROUTINE advec_w_pw_ij
Note: See TracChangeset for help on using the changeset viewer.