Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/large_scale_forcing_nudging_mod.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!
     
    7478    USE indices,                                                               &
    7579        ONLY:  nbgp, ngp_sums_ls, nx, nxl, nxlg, nxlu, nxr, nxrg, ny, nys,     &
    76                nysv, nysg, nyn, nyng, nzb, nz, nzt, wall_flags_static_0
     80               nysv, nysg, nyn, nyng, nzb, nz, nzt, wall_flags_total_0
    7781
    7882    USE kinds
     
    802806                                    ( td_lsa_lpt(k,nt+1) - td_lsa_lpt(k,nt) ) *&
    803807                                        MERGE( 1.0_wp, 0.0_wp,                 &
    804                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     808                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    805809                   ENDDO
    806810                ENDDO
     
    815819                                    ( td_lsa_q(k,nt+1) - td_lsa_q(k,nt) ) *    &
    816820                                        MERGE( 1.0_wp, 0.0_wp,                 &
    817                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     821                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    818822                   ENDDO
    819823                ENDDO
     
    836840                                     ( td_sub_lpt(k,nt+1) - td_sub_lpt(k,nt) )*&
    837841                                        MERGE( 1.0_wp, 0.0_wp,                 &
    838                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     842                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    839843                      ENDDO
    840844                   ENDDO
     
    849853                                       ( td_sub_q(k,nt+1) - td_sub_q(k,nt) ) * &
    850854                                        MERGE( 1.0_wp, 0.0_wp,                 &
    851                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     855                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    852856                      ENDDO
    853857                   ENDDO
     
    902906                             + fac * ( td_lsa_lpt(k,nt+1) - td_lsa_lpt(k,nt) )*&
    903907                                        MERGE( 1.0_wp, 0.0_wp,                 &
    904                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     908                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    905909             ENDDO
    906910
     
    911915                              + fac * ( td_lsa_q(k,nt+1) - td_lsa_q(k,nt) ) *  &
    912916                                        MERGE( 1.0_wp, 0.0_wp,                 &
    913                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     917                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    914918             ENDDO
    915919
     
    928932                                 ( td_sub_lpt(k,nt+1) - td_sub_lpt(k,nt) ) *   &
    929933                                        MERGE( 1.0_wp, 0.0_wp,                 &
    930                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     934                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    931935                ENDDO
    932936 
     
    937941                                 ( td_sub_q(k,nt+1) - td_sub_q(k,nt) ) *       &
    938942                                        MERGE( 1.0_wp, 0.0_wp,                 &
    939                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     943                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    940944                ENDDO
    941945
     
    11811185                      tend(k,j,i) = tend(k,j,i) + tmp_tend *                   &
    11821186                                        MERGE( 1.0_wp, 0.0_wp,                 &
    1183                                                BTEST( wall_flags_static_0(k,j,i), 1 ) )
     1187                                        BTEST( wall_flags_total_0(k,j,i), 1 ) )
    11841188
    11851189                      sums_ls_l(k,6) = sums_ls_l(k,6) + tmp_tend *             &
     
    12041208                      tend(k,j,i) = tend(k,j,i) + tmp_tend *                   &
    12051209                                        MERGE( 1.0_wp, 0.0_wp,                 &
    1206                                                BTEST( wall_flags_static_0(k,j,i), 2 ) )
     1210                                        BTEST( wall_flags_total_0(k,j,i), 2 ) )
    12071211
    12081212                      sums_ls_l(k,7) = sums_ls_l(k,7) + tmp_tend *             &
     
    12271231                      tend(k,j,i) = tend(k,j,i) + tmp_tend *                   &
    12281232                                        MERGE( 1.0_wp, 0.0_wp,                 &
    1229                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     1233                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    12301234
    12311235                      sums_ls_l(k,4) = sums_ls_l(k,4) + tmp_tend *             &
     
    12501254                      tend(k,j,i) = tend(k,j,i) + tmp_tend *                   &
    12511255                                        MERGE( 1.0_wp, 0.0_wp,                 &
    1252                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     1256                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    12531257
    12541258                      sums_ls_l(k,5) = sums_ls_l(k,5) + tmp_tend *             &
     
    13161320                tend(k,j,i) = tend(k,j,i) + tmp_tend *                         &
    13171321                                        MERGE( 1.0_wp, 0.0_wp,                 &
    1318                                                BTEST( wall_flags_static_0(k,j,i), 1 ) )
     1322                                        BTEST( wall_flags_total_0(k,j,i), 1 ) )
    13191323
    13201324                sums_ls_l(k,6) = sums_ls_l(k,6) + tmp_tend                     &
     
    13331337                tend(k,j,i) = tend(k,j,i) + tmp_tend *                         &
    13341338                                        MERGE( 1.0_wp, 0.0_wp,                 &
    1335                                                BTEST( wall_flags_static_0(k,j,i), 2 ) )
     1339                                        BTEST( wall_flags_total_0(k,j,i), 2 ) )
    13361340
    13371341                sums_ls_l(k,7) = sums_ls_l(k,7) + tmp_tend                     &
     
    13501354                tend(k,j,i) = tend(k,j,i) + tmp_tend *                         &
    13511355                                        MERGE( 1.0_wp, 0.0_wp,                 &
    1352                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     1356                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    13531357
    13541358                sums_ls_l(k,4) = sums_ls_l(k,4) + tmp_tend                     &
     
    13681372                tend(k,j,i) = tend(k,j,i) + tmp_tend *                         &
    13691373                                        MERGE( 1.0_wp, 0.0_wp,                 &
    1370                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     1374                                        BTEST( wall_flags_total_0(k,j,i), 0 ) )
    13711375
    13721376                sums_ls_l(k,5) = sums_ls_l(k,5) + tmp_tend                     &
Note: See TracChangeset for help on using the changeset viewer.