Changeset 2317 for palm/trunk/SOURCE/plant_canopy_model_mod.f90
- Timestamp:
- Jul 20, 2017 5:27:19 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/plant_canopy_model_mod.f90
r2274 r2317 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! Get topography top index via Function call 23 23 ! 24 24 ! Former revisions: … … 128 128 USE indices, & 129 129 ONLY: nbgp, nxl, nxlg, nxlu, nxr, nxrg, nyn, nyng, nys, nysg, nysv, & 130 nz, nzb, nz b_max, nzt, wall_flags_0130 nz, nzb, nzt 131 131 132 132 USE kinds 133 134 USE surface_mod, & 135 ONLY: get_topography_top_index 133 136 134 137 … … 1016 1019 ! 1017 1020 !-- Determine topography-top index on u-grid 1018 k_wall = MAXLOC( & 1019 MERGE( 1, 0, & 1020 BTEST( wall_flags_0(nzb:nzb_max,j,i), 14 ) & 1021 ), DIM = 1 & 1022 ) - 1 1021 k_wall = get_topography_top_index( j, i, 'u' ) 1023 1022 DO k = k_wall+1, k_wall+pch_index 1024 1023 … … 1084 1083 ! 1085 1084 !-- Determine topography-top index on v-grid 1086 k_wall = MAXLOC( & 1087 MERGE( 1, 0, & 1088 BTEST( wall_flags_0(nzb:nzb_max,j,i), 16 ) & 1089 ), DIM = 1 & 1090 ) - 1 1085 k_wall = get_topography_top_index( j, i, 'v' ) 1086 1091 1087 DO k = k_wall+1, k_wall+pch_index 1092 1088 … … 1152 1148 ! 1153 1149 !-- Determine topography-top index on w-grid 1154 k_wall = MAXLOC( & 1155 MERGE( 1, 0, & 1156 BTEST( wall_flags_0(nzb:nzb_max,j,i), 18 ) & 1157 ), DIM = 1 & 1158 ) - 1 1150 k_wall = get_topography_top_index( j, i, 'w' ) 1151 1159 1152 DO k = k_wall+1, k_wall+pch_index-1 1160 1153 … … 1207 1200 ! 1208 1201 !-- Determine topography-top index on scalar-grid 1209 k_wall = MAXLOC( & 1210 MERGE( 1, 0, & 1211 BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 ) & 1212 ), DIM = 1 & 1213 ) - 1 1202 k_wall = get_topography_top_index( j, i, 's' ) 1203 1214 1204 DO k = k_wall+1, k_wall+pch_index 1215 1205 … … 1227 1217 ! 1228 1218 !-- Determine topography-top index on scalar-grid 1229 k_wall = MAXLOC( & 1230 MERGE( 1, 0, & 1231 BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 ) & 1232 ), DIM = 1 & 1233 ) - 1 1219 k_wall = get_topography_top_index( j, i, 's' ) 1220 1234 1221 DO k = k_wall+1, k_wall+pch_index 1235 1222 … … 1260 1247 ! 1261 1248 !-- Determine topography-top index on scalar-grid 1262 k_wall = MAXLOC( & 1263 MERGE( 1, 0, & 1264 BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 ) & 1265 ), DIM = 1 & 1266 ) - 1 1249 k_wall = get_topography_top_index( j, i, 's' ) 1250 1267 1251 DO k = k_wall+1, k_wall+pch_index 1268 1252 … … 1292 1276 ! 1293 1277 !-- Determine topography-top index on scalar-grid 1294 k_wall = MAXLOC( & 1295 MERGE( 1, 0, & 1296 BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 ) & 1297 ), DIM = 1 & 1298 ) - 1 1278 k_wall = get_topography_top_index( j, i, 's' ) 1279 1299 1280 DO k = k_wall+1, k_wall+pch_index 1300 1281 … … 1390 1371 ! 1391 1372 !-- Determine topography-top index on u-grid 1392 k_wall = MAXLOC( & 1393 MERGE( 1, 0, & 1394 BTEST( wall_flags_0(nzb:nzb_max,j,i), 14 ) & 1395 ), DIM = 1 & 1396 ) - 1 1373 k_wall = get_topography_top_index( j, i, 'u' ) 1374 1397 1375 DO k = k_wall+1, k_wall+pch_index 1398 1376 … … 1453 1431 ! 1454 1432 !-- Determine topography-top index on v-grid 1455 k_wall = MAXLOC( & 1456 MERGE( 1, 0, & 1457 BTEST( wall_flags_0(nzb:nzb_max,j,i), 16 ) & 1458 ), DIM = 1 & 1459 ) - 1 1433 k_wall = get_topography_top_index( j, i, 'v' ) 1434 1460 1435 DO k = k_wall+1, k_wall+pch_index 1461 1436 … … 1516 1491 ! 1517 1492 !-- Determine topography-top index on w-grid 1518 k_wall = MAXLOC( & 1519 MERGE( 1, 0, & 1520 BTEST( wall_flags_0(nzb:nzb_max,j,i), 18 ) & 1521 ), DIM = 1 & 1522 ) - 1 1493 k_wall = get_topography_top_index( j, i, 'w' ) 1494 1523 1495 DO k = k_wall+1, k_wall+pch_index-1 1524 1496 … … 1566 1538 ! 1567 1539 !-- Determine topography-top index on scalar grid 1568 k_wall = MAXLOC( & 1569 MERGE( 1, 0, & 1570 BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 ) & 1571 ), DIM = 1 & 1572 ) - 1 1540 k_wall = get_topography_top_index( j, i, 's' ) 1541 1573 1542 DO k = k_wall+1, k_wall+pch_index 1574 1543 kk = k - k_wall !- lad arrays are defined flat … … 1582 1551 ! 1583 1552 !-- Determine topography-top index on scalar grid 1584 k_wall = MAXLOC( & 1585 MERGE( 1, 0, & 1586 BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 ) & 1587 ), DIM = 1 & 1588 ) - 1 1553 k_wall = get_topography_top_index( j, i, 's' ) 1554 1589 1555 DO k = k_wall+1, k_wall+pch_index 1590 1556 … … 1611 1577 ! 1612 1578 !-- Determine topography-top index on scalar grid 1613 k_wall = MAXLOC( & 1614 MERGE( 1, 0, & 1615 BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 ) & 1616 ), DIM = 1 & 1617 ) - 1 1579 k_wall = get_topography_top_index( j, i, 's' ) 1580 1618 1581 DO k = k_wall+1, k_wall+pch_index 1619 1582 … … 1640 1603 ! 1641 1604 !-- Determine topography-top index on scalar grid 1642 k_wall = MAXLOC( & 1643 MERGE( 1, 0, & 1644 BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 ) & 1645 ), DIM = 1 & 1646 ) - 1 1605 k_wall = get_topography_top_index( j, i, 's' ) 1606 1647 1607 DO k = k_wall+1, k_wall+pch_index 1648 1608
Note: See TracChangeset
for help on using the changeset viewer.