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

Timestamp:

Dec 18, 2019 11:55:56 AM (4 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_w.f90 (modified) (6 diffs)

palm/trunk/SOURCE/diffusion_w.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
 ! 
@@ -69 +73 @@
            
        USE indices,                                                            &            
-           ONLY :  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_static_0
+           ONLY :  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_total_0
            
        USE kinds
@@ -101 +105 @@
        !$ACC PRIVATE(surf_e, surf_s, flag, kmxm, kmxp, kmym, kmyp) &
        !$ACC PRIVATE(mask_west, mask_east, mask_south, mask_north) &
-       !$ACC PRESENT(wall_flags_static_0, km) &
+       !$ACC PRESENT(wall_flags_total_0, km) &
        !$ACC PRESENT(u, v, w) &
        !$ACC PRESENT(ddzu, ddzw, rho_air, drho_air_zw) &
@@ -114 +118 @@
 !--             Predetermine flag to mask topography and wall-bounded grid points. 
                 flag       = MERGE( 1.0_wp, 0.0_wp,                            &
-                                    BTEST( wall_flags_static_0(k,j,i),   3 ) ) 
+                                    BTEST( wall_flags_total_0(k,j,i),   3 ) ) 
                 mask_east  = MERGE( 1.0_wp, 0.0_wp,                            &
-                                    BTEST( wall_flags_static_0(k,j,i+1), 3 ) )
+                                    BTEST( wall_flags_total_0(k,j,i+1), 3 ) )
                 mask_west  = MERGE( 1.0_wp, 0.0_wp,                            &
-                                    BTEST( wall_flags_static_0(k,j,i-1), 3 ) )
+                                    BTEST( wall_flags_total_0(k,j,i-1), 3 ) )
                 mask_south = MERGE( 1.0_wp, 0.0_wp,                            &
-                                    BTEST( wall_flags_static_0(k,j-1,i), 3 ) )
+                                    BTEST( wall_flags_total_0(k,j-1,i), 3 ) )
                 mask_north = MERGE( 1.0_wp, 0.0_wp,                            &
-                                    BTEST( wall_flags_static_0(k,j+1,i), 3 ) )
+                                    BTEST( wall_flags_total_0(k,j+1,i), 3 ) )
 !
 !--             Interpolate eddy diffusivities on staggered gridpoints
@@ -254 +258 @@
            
        USE indices,                                                            &            
-           ONLY :  nzb, nzt, wall_flags_static_0
+           ONLY :  nzb, nzt, wall_flags_total_0
            
        USE kinds
@@ -286 +290 @@
 !
 !--       Predetermine flag to mask topography and wall-bounded grid points. 
-          flag       = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i),   3 ) ) 
-          mask_east  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i+1), 3 ) )
-          mask_west  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i-1), 3 ) )
-          mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j-1,i), 3 ) )
-          mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j+1,i), 3 ) )
+          flag       = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i),   3 ) ) 
+          mask_east  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i+1), 3 ) )
+          mask_west  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i-1), 3 ) )
+          mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j-1,i), 3 ) )
+          mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j+1,i), 3 ) )
 !
 !--       Interpolate eddy diffusivities on staggered gridpoints