Changeset 4329 for palm/trunk/SOURCE/init_3d_model.f90
- Timestamp:
- Dec 10, 2019 3:46:36 PM (4 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/init_3d_model.f90 (modified) (14 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_3d_model.f90
r4286 r4329 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Renamed wall_flags_0 to wall_flags_static_0 28 ! 29 ! 4286 2019-10-30 16:01:14Z resler 27 30 ! implement new palm_date_time_mod 28 31 ! … … 774 777 !-- Set velocity components at non-atmospheric / oceanic grid points to 775 778 !-- zero. 776 u = MERGE( u, 0.0_wp, BTEST( wall_flags_ 0, 1 ) )777 v = MERGE( v, 0.0_wp, BTEST( wall_flags_ 0, 2 ) )778 w = MERGE( w, 0.0_wp, BTEST( wall_flags_ 0, 3 ) )779 u = MERGE( u, 0.0_wp, BTEST( wall_flags_static_0, 1 ) ) 780 v = MERGE( v, 0.0_wp, BTEST( wall_flags_static_0, 2 ) ) 781 w = MERGE( w, 0.0_wp, BTEST( wall_flags_static_0, 3 ) ) 779 782 ! 780 783 !-- Initialize surface variables, e.g. friction velocity, momentum … … 824 827 ! 825 828 !-- Set velocities back to zero 826 u = MERGE( u, 0.0_wp, BTEST( wall_flags_ 0, 1 ) )827 v = MERGE( v, 0.0_wp, BTEST( wall_flags_ 0, 2 ) )829 u = MERGE( u, 0.0_wp, BTEST( wall_flags_static_0, 1 ) ) 830 v = MERGE( v, 0.0_wp, BTEST( wall_flags_static_0, 2 ) ) 828 831 ! 829 832 !-- WARNING: The extra boundary conditions set after running the … … 868 871 ! 869 872 !-- Mask topography 870 u = MERGE( u, 0.0_wp, BTEST( wall_flags_ 0, 1 ) )871 v = MERGE( v, 0.0_wp, BTEST( wall_flags_ 0, 2 ) )873 u = MERGE( u, 0.0_wp, BTEST( wall_flags_static_0, 1 ) ) 874 v = MERGE( v, 0.0_wp, BTEST( wall_flags_static_0, 2 ) ) 872 875 ! 873 876 !-- Set initial horizontal velocities at the lowest computational grid … … 884 887 DO k = nzb, nzt 885 888 u(k,j,i) = MERGE( u(k,j,i), 0.0_wp, & 886 BTEST( wall_flags_ 0(k,j,i), 20 ) )889 BTEST( wall_flags_static_0(k,j,i), 20 ) ) 887 890 v(k,j,i) = MERGE( v(k,j,i), 0.0_wp, & 888 BTEST( wall_flags_ 0(k,j,i), 21 ) )891 BTEST( wall_flags_static_0(k,j,i), 21 ) ) 889 892 ENDDO 890 893 ENDDO … … 1310 1313 DO k = nzb, nzt 1311 1314 u(k,j,i) = MERGE( u(k,j,i), 0.0_wp, & 1312 BTEST( wall_flags_ 0(k,j,i), 1 ) )1315 BTEST( wall_flags_static_0(k,j,i), 1 ) ) 1313 1316 v(k,j,i) = MERGE( v(k,j,i), 0.0_wp, & 1314 BTEST( wall_flags_ 0(k,j,i), 2 ) )1317 BTEST( wall_flags_static_0(k,j,i), 2 ) ) 1315 1318 w(k,j,i) = MERGE( w(k,j,i), 0.0_wp, & 1316 BTEST( wall_flags_ 0(k,j,i), 3 ) )1319 BTEST( wall_flags_static_0(k,j,i), 3 ) ) 1317 1320 ENDDO 1318 1321 ENDDO … … 1395 1398 u_init(k) * dzw(k) & 1396 1399 * MERGE( 1.0_wp, 0.0_wp, & 1397 BTEST( wall_flags_ 0(k,j,nxr), 1 )&1400 BTEST( wall_flags_static_0(k,j,nxr), 1 )& 1398 1401 ) 1399 1402 1400 1403 volume_flow_area_l(1) = volume_flow_area_l(1) + dzw(k) & 1401 1404 * MERGE( 1.0_wp, 0.0_wp, & 1402 BTEST( wall_flags_ 0(k,j,nxr), 1 )&1405 BTEST( wall_flags_static_0(k,j,nxr), 1 )& 1403 1406 ) 1404 1407 ENDDO … … 1412 1415 v_init(k) * dzw(k) & 1413 1416 * MERGE( 1.0_wp, 0.0_wp, & 1414 BTEST( wall_flags_ 0(k,nyn,i), 2 )&1417 BTEST( wall_flags_static_0(k,nyn,i), 2 )& 1415 1418 ) 1416 1419 volume_flow_area_l(2) = volume_flow_area_l(2) + dzw(k) & 1417 1420 * MERGE( 1.0_wp, 0.0_wp, & 1418 BTEST( wall_flags_ 0(k,nyn,i), 2 )&1421 BTEST( wall_flags_static_0(k,nyn,i), 2 )& 1419 1422 ) 1420 1423 ENDDO … … 1444 1447 hom_sum(k,1,0) * dzw(k) & 1445 1448 * MERGE( 1.0_wp, 0.0_wp, & 1446 BTEST( wall_flags_ 0(k,j,nx), 1 ) &1449 BTEST( wall_flags_static_0(k,j,nx), 1 ) & 1447 1450 ) 1448 1451 volume_flow_area_l(1) = volume_flow_area_l(1) + dzw(k) & 1449 1452 * MERGE( 1.0_wp, 0.0_wp, & 1450 BTEST( wall_flags_ 0(k,j,nx), 1 ) &1453 BTEST( wall_flags_static_0(k,j,nx), 1 ) & 1451 1454 ) 1452 1455 ENDDO … … 1460 1463 hom_sum(k,2,0) * dzw(k) & 1461 1464 * MERGE( 1.0_wp, 0.0_wp, & 1462 BTEST( wall_flags_ 0(k,ny,i), 2 ) &1465 BTEST( wall_flags_static_0(k,ny,i), 2 ) & 1463 1466 ) 1464 1467 volume_flow_area_l(2) = volume_flow_area_l(2) + dzw(k) & 1465 1468 * MERGE( 1.0_wp, 0.0_wp, & 1466 BTEST( wall_flags_ 0(k,ny,i), 2 ) &1469 BTEST( wall_flags_static_0(k,ny,i), 2 ) & 1467 1470 ) 1468 1471 ENDDO … … 1492 1495 u(k,j,nx) * dzw(k) & 1493 1496 * MERGE( 1.0_wp, 0.0_wp, & 1494 BTEST( wall_flags_ 0(k,j,nx), 1 ) &1497 BTEST( wall_flags_static_0(k,j,nx), 1 ) & 1495 1498 ) 1496 1499 volume_flow_area_l(1) = volume_flow_area_l(1) + dzw(k) & 1497 1500 * MERGE( 1.0_wp, 0.0_wp, & 1498 BTEST( wall_flags_ 0(k,j,nx), 1 ) &1501 BTEST( wall_flags_static_0(k,j,nx), 1 ) & 1499 1502 ) 1500 1503 ENDDO … … 1508 1511 v(k,ny,i) * dzw(k) & 1509 1512 * MERGE( 1.0_wp, 0.0_wp, & 1510 BTEST( wall_flags_ 0(k,ny,i), 2 ) &1513 BTEST( wall_flags_static_0(k,ny,i), 2 ) & 1511 1514 ) 1512 1515 volume_flow_area_l(2) = volume_flow_area_l(2) + dzw(k) & 1513 1516 * MERGE( 1.0_wp, 0.0_wp, & 1514 BTEST( wall_flags_ 0(k,ny,i), 2 ) &1517 BTEST( wall_flags_static_0(k,ny,i), 2 ) & 1515 1518 ) 1516 1519 ENDDO … … 1693 1696 !-- xy-grid points above topography 1694 1697 ngp_2dh_outer_l(k,sr) = ngp_2dh_outer_l(k,sr) + & 1695 MERGE( 1, 0, BTEST( wall_flags_ 0(k,j,i), 24 ) )1698 MERGE( 1, 0, BTEST( wall_flags_static_0(k,j,i), 24 ) ) 1696 1699 1697 1700 ngp_2dh_s_inner_l(k,sr) = ngp_2dh_s_inner_l(k,sr) + & 1698 MERGE( 1, 0, BTEST( wall_flags_ 0(k,j,i), 22 ) )1701 MERGE( 1, 0, BTEST( wall_flags_static_0(k,j,i), 22 ) ) 1699 1702 1700 1703 ENDDO … … 1792 1795 !$ACC COPYIN(rho_air(nzb:nzt+1), rho_air_zw(nzb:nzt+1)) & 1793 1796 !$ACC COPYIN(ddzu(1:nzt+1), ddzw(1:nzt+1)) & 1794 !$ACC COPYIN(wall_flags_ 0(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &1797 !$ACC COPYIN(wall_flags_static_0(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) & 1795 1798 !$ACC COPYIN(bc_h(0:1)) & 1796 1799 !$ACC COPYIN(bc_h(0)%i(1:bc_h(0)%ns)) &
Note: See TracChangeset
for help on using the changeset viewer.
