Changeset 4346 for palm/trunk/SOURCE/nesting_offl_mod.f90
- Timestamp:
- Dec 18, 2019 11:55:56 AM (4 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/nesting_offl_mod.f90 (modified) (22 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/nesting_offl_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 ! … … 191 195 nysv, nysg, nyn, nyng, nzb, nz, nzt, & 192 196 topo_top_ind, & 193 wall_flags_ static_0197 wall_flags_total_0 194 198 195 199 USE kinds … … 866 870 volume_flow_l(1) = volume_flow_l(1) + u(k,j,i) * dzw(k) * dy & 867 871 * MERGE( 1.0_wp, 0.0_wp, & 868 BTEST( wall_flags_static_0(k,j,i), 1 ) )872 BTEST( wall_flags_total_0(k,j,i), 1 ) ) 869 873 ENDDO 870 874 ENDDO … … 876 880 volume_flow_l(1) = volume_flow_l(1) - u(k,j,i) * dzw(k) * dy & 877 881 * MERGE( 1.0_wp, 0.0_wp, & 878 BTEST( wall_flags_static_0(k,j,i), 1 ) )882 BTEST( wall_flags_total_0(k,j,i), 1 ) ) 879 883 ENDDO 880 884 ENDDO … … 886 890 volume_flow_l(2) = volume_flow_l(2) + v(k,j,i) * dzw(k) * dx & 887 891 * MERGE( 1.0_wp, 0.0_wp, & 888 BTEST( wall_flags_static_0(k,j,i), 2 ) )892 BTEST( wall_flags_total_0(k,j,i), 2 ) ) 889 893 ENDDO 890 894 ENDDO … … 896 900 volume_flow_l(2) = volume_flow_l(2) - v(k,j,i) * dzw(k) * dx & 897 901 * MERGE( 1.0_wp, 0.0_wp, & 898 BTEST( wall_flags_static_0(k,j,i), 2 ) )902 BTEST( wall_flags_total_0(k,j,i), 2 ) ) 899 903 ENDDO 900 904 ENDDO … … 1001 1005 fac_dt ) * & 1002 1006 MERGE( 1.0_wp, 0.0_wp, & 1003 BTEST( wall_flags_ static_0(k,j,0), 1 ) )1007 BTEST( wall_flags_total_0(k,j,0), 1 ) ) 1004 1008 u(k,j,-1) = u(k,j,0) 1005 1009 ENDDO … … 1013 1017 fac_dt ) * & 1014 1018 MERGE( 1.0_wp, 0.0_wp, & 1015 BTEST( wall_flags_ static_0(k,j,-1), 3 ) )1019 BTEST( wall_flags_total_0(k,j,-1), 3 ) ) 1016 1020 ENDDO 1017 1021 w(nzt,j,-1) = w(nzt-1,j,-1) … … 1024 1028 fac_dt ) * & 1025 1029 MERGE( 1.0_wp, 0.0_wp, & 1026 BTEST( wall_flags_ static_0(k,j,-1), 2 ) )1030 BTEST( wall_flags_total_0(k,j,-1), 2 ) ) 1027 1031 ENDDO 1028 1032 v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j,-1) … … 1080 1084 fac_dt ) * & 1081 1085 MERGE( 1.0_wp, 0.0_wp, & 1082 BTEST( wall_flags_ static_0(k,j,nxr+1), 1 ) )1086 BTEST( wall_flags_total_0(k,j,nxr+1), 1 ) ) 1083 1087 ENDDO 1084 1088 u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,j,nxr+1) … … 1090 1094 fac_dt ) * & 1091 1095 MERGE( 1.0_wp, 0.0_wp, & 1092 BTEST( wall_flags_ static_0(k,j,nxr+1), 3 ) )1096 BTEST( wall_flags_total_0(k,j,nxr+1), 3 ) ) 1093 1097 ENDDO 1094 1098 w(nzt,j,nxr+1) = w(nzt-1,j,nxr+1) … … 1101 1105 fac_dt ) * & 1102 1106 MERGE( 1.0_wp, 0.0_wp, & 1103 BTEST( wall_flags_ static_0(k,j,nxr+1), 2 ) )1107 BTEST( wall_flags_total_0(k,j,nxr+1), 2 ) ) 1104 1108 ENDDO 1105 1109 v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,j,nxr+1) … … 1158 1162 fac_dt ) * & 1159 1163 MERGE( 1.0_wp, 0.0_wp, & 1160 BTEST( wall_flags_ static_0(k,0,i), 2 ) )1164 BTEST( wall_flags_total_0(k,0,i), 2 ) ) 1161 1165 v(k,-1,i) = v(k,0,i) 1162 1166 ENDDO … … 1170 1174 fac_dt ) * & 1171 1175 MERGE( 1.0_wp, 0.0_wp, & 1172 BTEST( wall_flags_ static_0(k,-1,i), 3 ) )1176 BTEST( wall_flags_total_0(k,-1,i), 3 ) ) 1173 1177 ENDDO 1174 1178 w(nzt,-1,i) = w(nzt-1,-1,i) … … 1181 1185 fac_dt ) * & 1182 1186 MERGE( 1.0_wp, 0.0_wp, & 1183 BTEST( wall_flags_ static_0(k,-1,i), 1 ) )1187 BTEST( wall_flags_total_0(k,-1,i), 1 ) ) 1184 1188 ENDDO 1185 1189 u_ref_l(nzb+1:nzt) = u_ref_l(nzb+1:nzt) + u(nzb+1:nzt,-1,i) … … 1239 1243 fac_dt ) * & 1240 1244 MERGE( 1.0_wp, 0.0_wp, & 1241 BTEST( wall_flags_static_0(k,nyn+1,i), 2 ) )1245 BTEST( wall_flags_total_0(k,nyn+1,i), 2 ) ) 1242 1246 ENDDO 1243 1247 v_ref_l(nzb+1:nzt) = v_ref_l(nzb+1:nzt) + v(nzb+1:nzt,nyn+1,i) … … 1249 1253 fac_dt ) * & 1250 1254 MERGE( 1.0_wp, 0.0_wp, & 1251 BTEST( wall_flags_static_0(k,nyn+1,i), 3 ) )1255 BTEST( wall_flags_total_0(k,nyn+1,i), 3 ) ) 1252 1256 ENDDO 1253 1257 w(nzt,nyn+1,i) = w(nzt-1,nyn+1,i) … … 1260 1264 fac_dt ) * & 1261 1265 MERGE( 1.0_wp, 0.0_wp, & 1262 BTEST( wall_flags_static_0(k,nyn+1,i), 1 ) )1266 BTEST( wall_flags_total_0(k,nyn+1,i), 1 ) ) 1263 1267 1264 1268 ENDDO … … 1318 1322 fac_dt ) * & 1319 1323 MERGE( 1.0_wp, 0.0_wp, & 1320 BTEST( wall_flags_ static_0(nzt+1,j,i), 1 ) )1324 BTEST( wall_flags_total_0(nzt+1,j,i), 1 ) ) 1321 1325 u_ref_l(nzt+1) = u_ref_l(nzt+1) + u(nzt+1,j,i) 1322 1326 ENDDO … … 1334 1338 fac_dt ) * & 1335 1339 MERGE( 1.0_wp, 0.0_wp, & 1336 BTEST( wall_flags_ static_0(nzt+1,j,i), 2 ) )1340 BTEST( wall_flags_total_0(nzt+1,j,i), 2 ) ) 1337 1341 v_ref_l(nzt+1) = v_ref_l(nzt+1) + v(nzt+1,j,i) 1338 1342 ENDDO … … 1350 1354 fac_dt ) * & 1351 1355 MERGE( 1.0_wp, 0.0_wp, & 1352 BTEST( wall_flags_ static_0(nzt,j,i), 3 ) )1356 BTEST( wall_flags_total_0(nzt,j,i), 3 ) ) 1353 1357 w(nzt+1,j,i) = w(nzt,j,i) 1354 1358 ENDDO … … 2133 2137 !-- After boundary data is initialized, mask topography at the 2134 2138 !-- boundaries for the velocity components. 2135 u = MERGE( u, 0.0_wp, BTEST( wall_flags_ static_0, 1 ) )2136 v = MERGE( v, 0.0_wp, BTEST( wall_flags_ static_0, 2 ) )2137 w = MERGE( w, 0.0_wp, BTEST( wall_flags_ static_0, 3 ) )2139 u = MERGE( u, 0.0_wp, BTEST( wall_flags_total_0, 1 ) ) 2140 v = MERGE( v, 0.0_wp, BTEST( wall_flags_total_0, 2 ) ) 2141 w = MERGE( w, 0.0_wp, BTEST( wall_flags_total_0, 3 ) ) 2138 2142 ! 2139 2143 !-- Initial calculation of the boundary layer depth from the prescribed
Note: See TracChangeset
for help on using the changeset viewer.
