Changeset 4346 for palm/trunk/SOURCE/diffusion_u.f90

Timestamp:

Dec 18, 2019 11:55:56 AM (5 years ago)

Author:

motisi

Message:

Introduction of wall_flags_total_0, which currently sets bits based on static topography information used in wall_flags_static_0

File:

• palm/trunk/SOURCE/diffusion_u.f90 (modified) (8 diffs)

palm/trunk/SOURCE/diffusion_u.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Introduction of wall_flags_total_0, which currently sets bits based on static
+! topography information used in wall_flags_static_0
+! 
+! 4329 2019-12-10 15:46:36Z motisi
 ! Renamed wall_flags_0 to wall_flags_static_0
 ! 
@@ -75 +79 @@
        
        USE indices,                                                            &
-           ONLY:  nxlu, nxr, nyn, nys, nzb, nzt, wall_flags_static_0
+           ONLY:  nxlu, nxr, nyn, nys, nzb, nzt, wall_flags_total_0
      
        USE kinds
@@ -108 +112 @@
        !$ACC PRIVATE(surf_e, surf_s, flag, kmym, kmyp, kmzm, kmzp) &
        !$ACC PRIVATE(mask_bottom, mask_north, mask_south, mask_top) &
-       !$ACC PRESENT(wall_flags_static_0, km) &
+       !$ACC PRESENT(wall_flags_total_0, km) &
        !$ACC PRESENT(u, v, w) &
        !$ACC PRESENT(ddzu, ddzw, drho_air, rho_air_zw) &
@@ -124 +128 @@
 !--             It is sufficient to masked only north- and south-facing surfaces, which
 !--             need special treatment for the u-component. 
-                flag       = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i),   1 ) ) 
-                mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j-1,i), 1 ) )
-                mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j+1,i), 1 ) )
+                flag       = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i),   1 ) ) 
+                mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j-1,i), 1 ) )
+                mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j+1,i), 1 ) )
 !
 !--             Interpolate eddy diffusivities on staggered gridpoints
@@ -199 +203 @@
 !--             momentum flux at model top.  
                 mask_bottom = MERGE( 1.0_wp, 0.0_wp,                           &
-                                     BTEST( wall_flags_static_0(k-1,j,i), 8 ) ) 
+                                 BTEST( wall_flags_total_0(k-1,j,i), 8 ) ) 
                 mask_top    = MERGE( 1.0_wp, 0.0_wp,                           &
-                                     BTEST( wall_flags_static_0(k+1,j,i), 8 ) ) *     &
+                                 BTEST( wall_flags_total_0(k+1,j,i), 8 ) ) *   &
                               MERGE( 1.0_wp, 0.0_wp,                           &
-                                     BTEST( wall_flags_static_0(k+1,j,i), 9 ) ) 
+                                 BTEST( wall_flags_total_0(k+1,j,i), 9 ) ) 
                 flag        = MERGE( 1.0_wp, 0.0_wp,                           &
-                                     BTEST( wall_flags_static_0(k,j,i), 1 ) ) 
+                                 BTEST( wall_flags_total_0(k,j,i), 1 ) ) 
 !
 !--             Interpolate eddy diffusivities on staggered gridpoints
@@ -322 +326 @@
        
        USE indices,                                                            &
-           ONLY:  nzb, nzt, wall_flags_static_0
+           ONLY:  nzb, nzt, wall_flags_total_0
       
        USE kinds
@@ -356 +360 @@
 !--       It is sufficient to masked only north- and south-facing surfaces, which
 !--       need special treatment for the u-component. 
-          flag       = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i),   1 ) ) 
-          mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j-1,i), 1 ) )
-          mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j+1,i), 1 ) )
+          flag       = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i),   1 ) ) 
+          mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j-1,i), 1 ) )
+          mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j+1,i), 1 ) )
 !
 !--       Interpolate eddy diffusivities on staggered gridpoints
@@ -425 +429 @@
 !--       momentum flux at model top. 
           mask_bottom = MERGE( 1.0_wp, 0.0_wp,                                 &
-                               BTEST( wall_flags_static_0(k-1,j,i), 8 ) ) 
+                               BTEST( wall_flags_total_0(k-1,j,i), 8 ) ) 
           mask_top    = MERGE( 1.0_wp, 0.0_wp,                                 &
-                               BTEST( wall_flags_static_0(k+1,j,i), 8 ) ) *           &
+                               BTEST( wall_flags_total_0(k+1,j,i), 8 ) ) *     &
                         MERGE( 1.0_wp, 0.0_wp,                                 &
-                               BTEST( wall_flags_static_0(k+1,j,i), 9 ) )
+                               BTEST( wall_flags_total_0(k+1,j,i), 9 ) )
           flag        = MERGE( 1.0_wp, 0.0_wp,                                 &
-                               BTEST( wall_flags_static_0(k,j,i), 1 ) ) 
+                               BTEST( wall_flags_total_0(k,j,i), 1 ) ) 
 !
 !--       Interpolate eddy diffusivities on staggered gridpoints