Ignore:
Timestamp:
Jun 12, 2019 4:50:15 PM (2 years ago)
Author:
suehring
Message:

Improve albedo initialization in land-surface model by setting always some default albedo types, in order to be independent on any user settings; Mask topography at boundary grid cells in nesting mass conservation

File:
1 edited

Legend:

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

    r4011 r4026  
    2525! -----------------
    2626! $Id$
     27! Masked topography at boundary grid points in mass conservation, in order to
     28! avoid that mean velocities within topography are imposed
     29!
     30! 4011 2019-05-31 14:34:03Z hellstea
    2731! Mass (volume) flux correction included to ensure global mass conservation for child domains.
    2832!
     
    52885292          DO  j = nys, nyn
    52895293             DO  k = nzb + 1, nzt
    5290                 u(k,j,i) = u(k,j,i) + u_corr_left
     5294                u(k,j,i) = u(k,j,i) + u_corr_left                              &
     5295                     * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 1 ) )
    52915296             ENDDO
    52925297          ENDDO
     
    52995304          DO  j = nys, nyn
    53005305             DO  k = nzb + 1, nzt
    5301                 u(k,j,i) = u(k,j,i) + u_corr_right
     5306                u(k,j,i) = u(k,j,i) + u_corr_right                              &
     5307                      * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 1 ) )
    53025308             ENDDO
    53035309          ENDDO
     
    53105316          DO  j = nysg, nys
    53115317             DO  k = nzb + 1, nzt
    5312                 v(k,j,i) = v(k,j,i) + v_corr_south
     5318                v(k,j,i) = v(k,j,i) + v_corr_south                              &
     5319                      * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 2 ) )
    53135320             ENDDO
    53145321          ENDDO
     
    53215328          DO  j = nyn+1, nyng
    53225329             DO  k = nzb + 1, nzt
    5323                 v(k,j,i) = v(k,j,i) + v_corr_north
     5330                v(k,j,i) = v(k,j,i) + v_corr_north                              &
     5331                      * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 2 ) )
    53245332             ENDDO
    53255333          ENDDO
Note: See TracChangeset for help on using the changeset viewer.