Changeset 4346 for palm/trunk/SOURCE/large_scale_forcing_nudging_mod.f90
- Timestamp:
- Dec 18, 2019 11:55:56 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/large_scale_forcing_nudging_mod.f90
r4329 r4346 25 25 ! ----------------- 26 26 ! $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 27 31 ! Renamed wall_flags_0 to wall_flags_static_0 28 32 ! … … 74 78 USE indices, & 75 79 ONLY: nbgp, ngp_sums_ls, nx, nxl, nxlg, nxlu, nxr, nxrg, ny, nys, & 76 nysv, nysg, nyn, nyng, nzb, nz, nzt, wall_flags_ static_080 nysv, nysg, nyn, nyng, nzb, nz, nzt, wall_flags_total_0 77 81 78 82 USE kinds … … 802 806 ( td_lsa_lpt(k,nt+1) - td_lsa_lpt(k,nt) ) *& 803 807 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 ) ) 805 809 ENDDO 806 810 ENDDO … … 815 819 ( td_lsa_q(k,nt+1) - td_lsa_q(k,nt) ) * & 816 820 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 ) ) 818 822 ENDDO 819 823 ENDDO … … 836 840 ( td_sub_lpt(k,nt+1) - td_sub_lpt(k,nt) )*& 837 841 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 ) ) 839 843 ENDDO 840 844 ENDDO … … 849 853 ( td_sub_q(k,nt+1) - td_sub_q(k,nt) ) * & 850 854 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 ) ) 852 856 ENDDO 853 857 ENDDO … … 902 906 + fac * ( td_lsa_lpt(k,nt+1) - td_lsa_lpt(k,nt) )*& 903 907 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 ) ) 905 909 ENDDO 906 910 … … 911 915 + fac * ( td_lsa_q(k,nt+1) - td_lsa_q(k,nt) ) * & 912 916 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 ) ) 914 918 ENDDO 915 919 … … 928 932 ( td_sub_lpt(k,nt+1) - td_sub_lpt(k,nt) ) * & 929 933 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 ) ) 931 935 ENDDO 932 936 … … 937 941 ( td_sub_q(k,nt+1) - td_sub_q(k,nt) ) * & 938 942 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 ) ) 940 944 ENDDO 941 945 … … 1181 1185 tend(k,j,i) = tend(k,j,i) + tmp_tend * & 1182 1186 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 ) ) 1184 1188 1185 1189 sums_ls_l(k,6) = sums_ls_l(k,6) + tmp_tend * & … … 1204 1208 tend(k,j,i) = tend(k,j,i) + tmp_tend * & 1205 1209 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 ) ) 1207 1211 1208 1212 sums_ls_l(k,7) = sums_ls_l(k,7) + tmp_tend * & … … 1227 1231 tend(k,j,i) = tend(k,j,i) + tmp_tend * & 1228 1232 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 ) ) 1230 1234 1231 1235 sums_ls_l(k,4) = sums_ls_l(k,4) + tmp_tend * & … … 1250 1254 tend(k,j,i) = tend(k,j,i) + tmp_tend * & 1251 1255 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 ) ) 1253 1257 1254 1258 sums_ls_l(k,5) = sums_ls_l(k,5) + tmp_tend * & … … 1316 1320 tend(k,j,i) = tend(k,j,i) + tmp_tend * & 1317 1321 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 ) ) 1319 1323 1320 1324 sums_ls_l(k,6) = sums_ls_l(k,6) + tmp_tend & … … 1333 1337 tend(k,j,i) = tend(k,j,i) + tmp_tend * & 1334 1338 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 ) ) 1336 1340 1337 1341 sums_ls_l(k,7) = sums_ls_l(k,7) + tmp_tend & … … 1350 1354 tend(k,j,i) = tend(k,j,i) + tmp_tend * & 1351 1355 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 ) ) 1353 1357 1354 1358 sums_ls_l(k,4) = sums_ls_l(k,4) + tmp_tend & … … 1368 1372 tend(k,j,i) = tend(k,j,i) + tmp_tend * & 1369 1373 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 ) ) 1371 1375 1372 1376 sums_ls_l(k,5) = sums_ls_l(k,5) + tmp_tend &
Note: See TracChangeset
for help on using the changeset viewer.