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_s_pw.f90

    r2101 r2232  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! topography representation via flags
    2323!
    2424! Former revisions:
     
    104104
    105105       USE indices,                                                            &
    106            ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzb_s_inner,&
    107                   nzt
     106           ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb,             &
     107                  nzt, wall_flags_0
    108108
    109109       USE kinds
     
    125125       DO  i = nxl, nxr
    126126          DO  j = nys, nyn
    127              DO  k = nzb_s_inner(j,i)+1, nzt
    128                 tend(k,j,i) = tend(k,j,i)                                      &
    129               -0.5_wp * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &
     127             DO  k = nzb+1, nzt
     128                tend(k,j,i) = tend(k,j,i)                                      &
     129            ( -0.5_wp * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &
    130130                        - ( u(k,j,i)   - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) &
    131131                        ) * ddx                                                   &
     
    135135              -         (   w(k,j,i)                * ( sk(k+1,j,i) - sk(k,j,i) ) &
    136136                        -   w(k-1,j,i)              * ( sk(k-1,j,i) - sk(k,j,i) ) &
    137                         ) * dd2zu(k)
     137                        ) * dd2zu(k)                                              &
     138            ) * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    138139             ENDDO
    139140          ENDDO
     
    160161
    161162       USE indices,                                                            &
    162            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_s_inner, nzt
     163           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    163164
    164165       USE kinds
     
    178179
    179180
    180        DO  k = nzb_s_inner(j,i)+1, nzt
    181           tend(k,j,i) = tend(k,j,i)                                            &
    182               -0.5_wp * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &
     181       DO  k = nzb+1, nzt
     182          tend(k,j,i) = tend(k,j,i)                                            &
     183            ( -0.5_wp * ( ( u(k,j,i+1) - u_gtrans ) * ( sk(k,j,i+1) - sk(k,j,i) ) &
    183184                        - ( u(k,j,i)   - u_gtrans ) * ( sk(k,j,i-1) - sk(k,j,i) ) &
    184185                        ) * ddx                                                   &
     
    188189              -         (   w(k,j,i)                * ( sk(k+1,j,i) - sk(k,j,i) ) &
    189190                        -   w(k-1,j,i)              * ( sk(k-1,j,i) - sk(k,j,i) ) &
    190                         ) * dd2zu(k)
     191                        ) * dd2zu(k)                                              &
     192            ) * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    191193       ENDDO
    192194
Note: See TracChangeset for help on using the changeset viewer.