Ignore:
Timestamp:
Dec 18, 2019 11:55:56 AM (22 months 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:
1 edited

Legend:

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

    r4329 r4346  
    2525! -----------------
    2626! $Id$
     27! Introduction of wall_flags_total_0, which currently sets bits based on static
     28! topography information used in wall_flags_static_0
     29!
     30! 4329 2019-12-10 15:46:36Z motisi
    2731! Renamed wall_flags_0 to wall_flags_static_0
    2832!
     
    6973           
    7074       USE indices,                                                            &           
    71            ONLY :  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_static_0
     75           ONLY :  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_total_0
    7276           
    7377       USE kinds
     
    101105       !$ACC PRIVATE(surf_e, surf_s, flag, kmxm, kmxp, kmym, kmyp) &
    102106       !$ACC PRIVATE(mask_west, mask_east, mask_south, mask_north) &
    103        !$ACC PRESENT(wall_flags_static_0, km) &
     107       !$ACC PRESENT(wall_flags_total_0, km) &
    104108       !$ACC PRESENT(u, v, w) &
    105109       !$ACC PRESENT(ddzu, ddzw, rho_air, drho_air_zw) &
     
    114118!--             Predetermine flag to mask topography and wall-bounded grid points.
    115119                flag       = MERGE( 1.0_wp, 0.0_wp,                            &
    116                                     BTEST( wall_flags_static_0(k,j,i),   3 ) )
     120                                    BTEST( wall_flags_total_0(k,j,i),   3 ) )
    117121                mask_east  = MERGE( 1.0_wp, 0.0_wp,                            &
    118                                     BTEST( wall_flags_static_0(k,j,i+1), 3 ) )
     122                                    BTEST( wall_flags_total_0(k,j,i+1), 3 ) )
    119123                mask_west  = MERGE( 1.0_wp, 0.0_wp,                            &
    120                                     BTEST( wall_flags_static_0(k,j,i-1), 3 ) )
     124                                    BTEST( wall_flags_total_0(k,j,i-1), 3 ) )
    121125                mask_south = MERGE( 1.0_wp, 0.0_wp,                            &
    122                                     BTEST( wall_flags_static_0(k,j-1,i), 3 ) )
     126                                    BTEST( wall_flags_total_0(k,j-1,i), 3 ) )
    123127                mask_north = MERGE( 1.0_wp, 0.0_wp,                            &
    124                                     BTEST( wall_flags_static_0(k,j+1,i), 3 ) )
     128                                    BTEST( wall_flags_total_0(k,j+1,i), 3 ) )
    125129!
    126130!--             Interpolate eddy diffusivities on staggered gridpoints
     
    254258           
    255259       USE indices,                                                            &           
    256            ONLY :  nzb, nzt, wall_flags_static_0
     260           ONLY :  nzb, nzt, wall_flags_total_0
    257261           
    258262       USE kinds
     
    286290!
    287291!--       Predetermine flag to mask topography and wall-bounded grid points.
    288           flag       = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i),   3 ) )
    289           mask_east  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i+1), 3 ) )
    290           mask_west  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j,i-1), 3 ) )
    291           mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j-1,i), 3 ) )
    292           mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_static_0(k,j+1,i), 3 ) )
     292          flag       = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i),   3 ) )
     293          mask_east  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i+1), 3 ) )
     294          mask_west  = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j,i-1), 3 ) )
     295          mask_south = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j-1,i), 3 ) )
     296          mask_north = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(k,j+1,i), 3 ) )
    293297!
    294298!--       Interpolate eddy diffusivities on staggered gridpoints
Note: See TracChangeset for help on using the changeset viewer.