Changeset 4026 for palm/trunk/SOURCE/pmc_interface_mod.f90

Timestamp:

Jun 12, 2019 4:50:15 PM (5 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:

• palm/trunk/SOURCE/pmc_interface_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Masked topography at boundary grid points in mass conservation, in order to 
+! avoid that mean velocities within topography are imposed
+! 
+! 4011 2019-05-31 14:34:03Z hellstea
 ! Mass (volume) flux correction included to ensure global mass conservation for child domains.
 ! 
@@ -5288 +5292 @@
           DO  j = nys, nyn
              DO  k = nzb + 1, nzt
-                u(k,j,i) = u(k,j,i) + u_corr_left
+                u(k,j,i) = u(k,j,i) + u_corr_left                              &
+                     * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 1 ) )
              ENDDO
           ENDDO
@@ -5299 +5304 @@
           DO  j = nys, nyn
              DO  k = nzb + 1, nzt
-                u(k,j,i) = u(k,j,i) + u_corr_right
+                u(k,j,i) = u(k,j,i) + u_corr_right                              &
+                      * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 1 ) )
              ENDDO
           ENDDO
@@ -5310 +5316 @@
           DO  j = nysg, nys
              DO  k = nzb + 1, nzt
-                v(k,j,i) = v(k,j,i) + v_corr_south
+                v(k,j,i) = v(k,j,i) + v_corr_south                              &
+                      * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 2 ) )
              ENDDO
           ENDDO
@@ -5321 +5328 @@
           DO  j = nyn+1, nyng
              DO  k = nzb + 1, nzt
-                v(k,j,i) = v(k,j,i) + v_corr_north
+                v(k,j,i) = v(k,j,i) + v_corr_north                              &
+                      * MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 2 ) )
              ENDDO
           ENDDO