Changeset 4329 for palm/trunk/SOURCE/prognostic_equations.f90
- Timestamp:
- Dec 10, 2019 3:46:36 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/prognostic_equations.f90
r4182 r4329 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Renamed wall_flags_0 to wall_flags_static_0 28 ! 29 ! 4182 2019-08-22 15:20:23Z scharf 27 30 ! Corrected "Former revisions" section 28 31 ! … … 203 206 ONLY: advc_flags_s, & 204 207 nbgp, nxl, nxlg, nxlu, nxr, nxrg, nyn, nyng, nys, nysg, nysv, & 205 nzb, nzt, wall_flags_ 0208 nzb, nzt, wall_flags_static_0 206 209 207 210 USE kinds … … 373 376 * ( u(k,j,i) - u_init(k) ) & 374 377 ) * MERGE( 1.0_wp, 0.0_wp, & 375 BTEST( wall_flags_ 0(k,j,i), 1 )&378 BTEST( wall_flags_static_0(k,j,i), 1 )& 376 379 ) 377 380 ENDDO … … 451 454 * ( v(k,j,i) - v_init(k) )& 452 455 ) * MERGE( 1.0_wp, 0.0_wp, & 453 BTEST( wall_flags_ 0(k,j,i), 2 )&456 BTEST( wall_flags_static_0(k,j,i), 2 )& 454 457 ) 455 458 ENDDO … … 525 528 - tsc(5) * rdf(k) * w(k,j,i) & 526 529 ) * MERGE( 1.0_wp, 0.0_wp, & 527 BTEST( wall_flags_ 0(k,j,i), 3 )&530 BTEST( wall_flags_static_0(k,j,i), 3 )& 528 531 ) 529 532 ENDDO … … 623 626 ) & 624 627 * MERGE( 1.0_wp, 0.0_wp, & 625 BTEST( wall_flags_ 0(k,j,i), 0 )&628 BTEST( wall_flags_static_0(k,j,i), 0 )& 626 629 ) 627 630 ENDDO … … 713 716 ) & 714 717 * MERGE( 1.0_wp, 0.0_wp, & 715 BTEST( wall_flags_ 0(k,j,i), 0 )&718 BTEST( wall_flags_static_0(k,j,i), 0 )& 716 719 ) 717 720 IF ( q_p(k,j,i) < 0.0_wp ) q_p(k,j,i) = 0.1_wp * q(k,j,i) … … 811 814 ) & 812 815 * MERGE( 1.0_wp, 0.0_wp, & 813 BTEST( wall_flags_ 0(k,j,i), 0 )&816 BTEST( wall_flags_static_0(k,j,i), 0 )& 814 817 ) 815 818 IF ( s_p(k,j,i) < 0.0_wp ) s_p(k,j,i) = 0.1_wp * s(k,j,i) … … 924 927 !-- Prognostic equation for u-velocity component 925 928 !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) & 926 !$ACC PRESENT(u, tend, tu_m, u_init, rdf, wall_flags_ 0) &929 !$ACC PRESENT(u, tend, tu_m, u_init, rdf, wall_flags_static_0) & 927 930 !$ACC PRESENT(tsc(2:5)) & 928 931 !$ACC PRESENT(u_p) … … 935 938 ( u(k,j,i) - u_init(k) ) & 936 939 ) * MERGE( 1.0_wp, 0.0_wp, & 937 BTEST( wall_flags_ 0(k,j,i), 1 ) &940 BTEST( wall_flags_static_0(k,j,i), 1 ) & 938 941 ) 939 942 ENDDO … … 1027 1030 !-- Prognostic equation for v-velocity component 1028 1031 !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) & 1029 !$ACC PRESENT(v, tend, tv_m, v_init, rdf, wall_flags_ 0) &1032 !$ACC PRESENT(v, tend, tv_m, v_init, rdf, wall_flags_static_0) & 1030 1033 !$ACC PRESENT(tsc(2:5)) & 1031 1034 !$ACC PRESENT(v_p) … … 1038 1041 ( v(k,j,i) - v_init(k) ) & 1039 1042 ) * MERGE( 1.0_wp, 0.0_wp, & 1040 BTEST( wall_flags_ 0(k,j,i), 2 )&1043 BTEST( wall_flags_static_0(k,j,i), 2 )& 1041 1044 ) 1042 1045 ENDDO … … 1126 1129 !-- Prognostic equation for w-velocity component 1127 1130 !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) & 1128 !$ACC PRESENT(w, tend, tw_m, v_init, rdf, wall_flags_ 0) &1131 !$ACC PRESENT(w, tend, tw_m, v_init, rdf, wall_flags_static_0) & 1129 1132 !$ACC PRESENT(tsc(2:5)) & 1130 1133 !$ACC PRESENT(w_p) … … 1136 1139 - tsc(5) * rdf(k) * w(k,j,i) & 1137 1140 ) * MERGE( 1.0_wp, 0.0_wp, & 1138 BTEST( wall_flags_ 0(k,j,i), 3 )&1141 BTEST( wall_flags_static_0(k,j,i), 3 )& 1139 1142 ) 1140 1143 ENDDO … … 1264 1267 !-- Prognostic equation for potential temperature 1265 1268 !$ACC PARALLEL LOOP COLLAPSE(3) PRIVATE(i, j, k) & 1266 !$ACC PRESENT(pt, tend, tpt_m, wall_flags_ 0) &1269 !$ACC PRESENT(pt, tend, tpt_m, wall_flags_static_0) & 1267 1270 !$ACC PRESENT(pt_init, rdf_sc, ptdf_x, ptdf_y) & 1268 1271 !$ACC PRESENT(tsc(3:5)) & … … 1278 1281 ) & 1279 1282 * MERGE( 1.0_wp, 0.0_wp, & 1280 BTEST( wall_flags_ 0(k,j,i), 0 ) &1283 BTEST( wall_flags_static_0(k,j,i), 0 ) & 1281 1284 ) 1282 1285 ENDDO … … 1399 1402 ( q(k,j,i) - q_init(k) ) & 1400 1403 ) * MERGE( 1.0_wp, 0.0_wp, & 1401 BTEST( wall_flags_ 0(k,j,i), 0 ) &1404 BTEST( wall_flags_static_0(k,j,i), 0 ) & 1402 1405 ) 1403 1406 IF ( q_p(k,j,i) < 0.0_wp ) q_p(k,j,i) = 0.1_wp * q(k,j,i) … … 1520 1523 ) & 1521 1524 * MERGE( 1.0_wp, 0.0_wp, & 1522 BTEST( wall_flags_ 0(k,j,i), 0 ) &1525 BTEST( wall_flags_static_0(k,j,i), 0 ) & 1523 1526 ) 1524 1527 IF ( s_p(k,j,i) < 0.0_wp ) s_p(k,j,i) = 0.1_wp * s(k,j,i)
Note: See TracChangeset
for help on using the changeset viewer.